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Chapter 1 

INTRODUCTION 



Fullerenes were first discovered in 1985 in an apparatus designed by 
Rick Smalley to produce atomic clusters of non-volatile elements [1]. The 
original aim of the experiment was to look for long chain HC„N molecules 
that Harry Kroto thought could be formed under the laser vapourisation 
conditions of the experiment - conditions not unlike those in carbon stars [2]. 
The discovery of this new form of pure carbon, and the insight that led to the 
interpretation of the magic number features in the observed mass spectra in 
terms of hollow carbon molecules, led to the award of the Nobel Prize in 
Chemistry to Curl, Kroto and Smalley in 1996. Many interesting review 
articles and books have been written about the events surrounding the 
discovery [2-8]. However, the fullerenes would have remained a curiosity of 
interest only to a small community of chemical physicists studying the 
properties of atomic clusters if the breakthrough achieved by Kratschmer 
and Huffman in 1990 had not occurred [9]. While studying the optical 
properties of various kinds of graphitic soots, these scientists discovered a 
remarkably simple way of producing macroscopic amounts of fullerenes. In 
addition to the production, they also showed that it was possible to separate 
the stable fullerenes (in particular Ceo and C70) from the rest of the soot, to 
purify them by standard chemical methods and to produce sufficient material 
to confirm the original hollow "soccerball" suggestion for the structure of 
Ceo- 

Since then there have been very rapid developments in the field and 
fullerenes have made an impact in many different scientific disciplines. This 
monograph is just one in a series of books devoted to the properties of 
fullerenes. Fullerenes have perhaps made their most impressive impact, 
relative to the overall size of the research activities in a particular field 
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worldwide, in the area of their birth: chemical physics and the closely related 
areas of atomic and molecular physics. The attraction of having a mass- 
selected cluster with the highest possible symmetry (Ih) that can easily be 
sublimed into the gas phase has inspired many groups with different 
experimental backgrounds to investigate the properties of these molecules. 
For many researchers from the atomic and molecular physics community 
this has opened up a whole new research area by providing a convenient and 
attractive model system for studying the behaviour of complex molecular 
systems with a large number of degrees of freedom. One nice example of 
this is a fullerene beam diffraction experiment showing quantum interference 
behaviour from Ceo [10]. In many cases the experience gained from 
experiments with fullerenes is now beginning to be applied to even more 
complex biomolecular systems. From the theoretical point of view, the 
fascinating structure and high symmetry of buckminsterfullerene has also 
inspired many works on the structural, electronic and dynamical properties 
of the molecule. A subset of the fullerene activities within atomic, chemical 
and molecular physics: collision experiments with fullerenes, is the topic of 
this book. 

The first two chapters following this introduction give a general 
introduction to the most common experimental techniques and theoretical 
models applied to the study of fullerene collisions. They supply background 
information that will be useful for understanding the discussions in later 
chapters. Chapters 4-10 then give overviews of the various sub-topics of the 
field. The aim is not to give a full comprehensive coverage of all work that 
has been carried out in any given topic but to pick out the most informative 
and illustrative examples and to point the reader towards more detailed 
presentations. On the whole, emphasis will be placed on collisions carried 
out under well-defined single-collision conditions. The topics range from 
reactive and thermal energy collisions (Chap 4) through collision induced 
dissociation (Chap 5), charge transfer (Chap 6) and electron impact studies 
(Chap 7). Chapters 8 and 9 are devoted to the more exotic members of the 
fullerene collision family: collisional production of endohedral fullerenes 
and fullerene-fullerene collisions. Finally, we finish in Chap. 10 with a 
discussion of the dynamical effects observed when fullerenes are collided 
from surfaces. 
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Chapter 2 

EXPERIMENTAL TECHNIQUES 



Fullerenes have been very fruitful systems for studying collision 
dynamics. There are a number of reasons for this. One is the sheer beauty of 
the structure and the attraction of a spherical molecule. This has enticed a 
number of atomic and nuclear physics research groups into the area of 
molecular collisions and has opened up many new possibilities for applying 
the well-developed collision techniques of these well-established disciplines 
to study the dynamical behaviour of complex molecular systems. Fullerenes 
have provided an attractive and convenient stepping stone from the world of 
atomic and simple molecular collisions to large (bio)molecules. From a more 
pragmatic point of view, fullerenes are very easy to handle experimentally 
and this has of course also played a major role in the development of our 
understanding of these materials. Since molecular beams of fullerenes can be 
obtained with relative ease, as described below, this leaves more room for 
experimental complexity and creativity in the design of the collision 
experiments that can be carried out. The following chapter discusses the way 
fullerenes are used in collision experiments and explains the basic 
techniques that are employed. Examples are given of typical experimental 
set-ups covering a wide range in collision energies from thermal reactivity 
studies to MeV surface impact. Collision experiments require ultra-high 
vacuum (< 10 '° mbar) or good high vacuum (< 10’’ mbar) conditions since 
one wants to minimise the effects of collisions with background gas and is 
generally interested in the outcome of single collisions. 
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1. FULLERENE SOURCES 

Fullerenes are prepared using the arc discharge method discovered by 
Kratschmer and Huffman in 1990 [1]. The soluble soot extract, which 
contains Ceo, C70 and some higher fullerenes [2] is further purified by liquid 
chromatography (HPLC). Fullerenes are commercially available in varying 
degrees of purity from raw soluble soot material to "Gold grade" Ceo with a 
nominal purity of 99.4 % (by weight). Many of the experiments that will be 
considered in this chapter and later in the book involve ionisation and mass 
selection of fullerene ions. For these purposes one does not require a highly 
purified sample but can use a mixture of extractable fullerenes since they 
will, in any case, be mass analysed before the collision. If neutral fullerenes 
are required - either as a molecular beam or as an effusive target, then the 
purity is more critical. Even highly purified commercial fullerene powder (> 
99.4 % Ceo) still contains a significant amount of impurities consisting of 
higher fullerenes and, particularly, solvent rest and possibly even compounds 
formed on passage through the HPLC purification column. The presence of 
these impurities also affects the vapour pressure and reproducibility of 
experimental results as will be discussed below. It is always advisable to 
carry out one or more sublimation cycles on commercial material, however 
this has only rarely been done in the experiments reported in the literature. 
Normally, the fullerene material will just be “out-gassed” at ca. 300 °C 
overnight before any experiments are carried out. The vapour pressure rises 
most steeply on the first increase in temperature (mainly solvent out-gassing) 
but on multiple heating and cooling cycles the achievable vapour pressure 
decreases. This effect cannot be solely due to solvent out-gassing. Many 
groups have also reported the formation of insoluble rest material in the oven 
after heating. The insoluble material still shows the presence of fullerenes if 
studied by e.g. laser desorption mass spectrometry and presumably has 
undergone a polymerisation reaction induced either by impurities in the 
powder and/or catalytic reactions involving trace elements from the oven 
walls. 

1.1 Neutral Beams 

The simplest way of producing a molecular beam of neutral fullerenes in 
a vacuum apparatus is to place the purified powder in a small oven with an 
aperture of a few mm and heat to a temperature of 350-500 °C. If a nozzle 
with a high aspect ratio is used with the oven then multiple collisions within 
the nozzle lead to a reasonably well-directed beam with the divergence 
decreasing with increasing nozzle length. The Kolodney group has 
developed an effusive source using ceramic material for a two-stage nozzle 
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that can be heated to temperatures of around 2100 K [3]. This allows the 
temperature of the molecular beam to be varied over a relatively wide range 
that can be an important tool for determining the influence of thermal 
effects. 

Due to the large number of vibrational degrees of freedom in Ceo (1V4), 
the molecules in beams produced by simply heating the fullerene powder to 
a temperature of around 500 °C or higher have a relatively high amount of 
vibrational excitation which can affect the outcome of collision experiments. 
Fig. 2.1 plots the vibrational energy as a function of Ceo temperature, 
calculated using the expression 



i 



eXp(-/jV; / kgT) 

[l-exTpi-hVi/kgT)] 



( 2 . 1 . 1 ) 



where g, is the degeneracy of the vibrational mode, v, is the mode 
frequency, h is Planck's constant, ks Boltzmann's constant and T the 
temperature. To a good approximation, the individual mode frequencies can 
be replaced by the average mode frequency (2.79x10^^ Hz). The vibrational 
frequencies have been taken from Schettino et al. [4] where a complete 
vibrational assignment was performed using IR, Raman and neutron 
scattering data together with density functional calculations. The average of 
the vibrational mode frequencies reported in Schettino et al., as used to 
obtain Fig. 2.1, are tabulated in Table 2.1. 




Figure 2-1. Relationship between vibrational energy and temperature for Qo, calculated using 
the vibrational frequencies in Table 2.1. The dashed line shows the calculated values using 
just the average vibrational frequency of 2.79x10*^ Hz instead of the individual mode 

frequencies. 
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Table 2-1. Frequencies of normal modes of vibration of in cm"'. Average values taken 
from Schetttino et al.[4]. Degeneracies of the modes are given in brackets. 



Symmetry 


Frequency [cm''] 


Symmetry 


Frequency [cm'*] 


Ag(l) 


491 


Audi 


973 




1458 


Tiu(3) 


527 


Tig (3) 


565 




578 




850 




1184 




1289 




1439 


Tzg (3) 


535 


Tzu (3) 


355 




764 




714 




805 




1040 




1334 




1187 


Gg(4) 


484 




1565 




570 


G„(4) 


344 




738 




760 




1135 




771 




1314 




967 




1515 




1321 


Hg(5) 


266 




1415 




429 


H„(5) 


403 




710 




486 




768 




669 




1095 




740 




1254 




1211 




1437 




1342 




1570 




1530 



It is rather difficult to reduce the vibrational excitation energy in a beam 
of neutral fullerenes below the temperature of evaporation, however, this has 
been achieved by using a source based on a cluster aggregation source where 
the fullerenes are evaporated into a stream of He gas cooled to liquid 
nitrogen temperatures [5]. The fullerene vibrational temperature is then also 
liquid nitrogen temperature (or colder, depending on the expansion 
conditions). Such a source is presently being used in collisions between 
fullerenes and highly charged ions to investigate the role of vibrational 
excitation in determining the stability of highly-charged fullerenes produced 
in the collisions [6] (see Chapter 6). 

1.2 Charged Beams 

In the early nineties beams of charged (positive or negative) fullerenes 
used for collision experiments were often produced by laser desorption. This 
method was then mainly superseded (due to stability and cost arguments) by 
electron impact ionisation of an effusive beam of fullerenes produced from a 
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small oven. Often, just a small tungsten filament biased at about 100 eV is 
enough to produce a relatively intense ion beam. The simplicity combined 
with the advantage of having either a continuous or pulsed ion source gives 
this method some clear advantages over the laser sources. Positive or 
negative ions can then be extracted from the source using standard 
techniques (some examples are given below). Many experiments that will be 
reported later in the book have used positively charged beams of fullerenes 
produced from electron cyclotron resonance (ECR) sources. Such sources 
were originally developed for producing highly charged ion beams and can 
be easily adapted for fullerenes [7]. An example of such a source is shown 
in Fig. 2.2. Beams of C6o‘*^ (q = 1-5) with intensities sufficiently high to 
carry out collision experiments can be routinely produced from ECR 
sources. 



Microwave 




Figure 2-2. Schematic diagram of ECR source of multiply charged fullerenes, adapted from 

18] 

The big disadvantage of all the above-mentioned ion sources is the 
uncertainty in the internal excitation energy of the fullerene ions produced. 
An indication of the amount of excitation can be obtained from studying 
metastable decay of the fullerene ions on the microsecond time scale. This 
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can be done, either in a conventional double sector mass spectrometer [9] or 
e.g. in areflectron time-of-flight mass spectrometer [10]. The internal energy 
of the fullerene ions as they enter the field-free region of the mass 
spectrometer can be estimated by determining the amount of unimolecular 
decay (generally, evaporation of C 2 molecules from the cage) that occurs 
within a time window set by the parameters of the experiment. In this way 
the rate constant for unimolecular decay can be directly determined and, 
assuming the mechanism and energy for the unimolecular decay process are 
known, this can be used to estimate a vibrational temperature [11]. 
Unfortunately, the decay mechanism and dissociation energy for the 
fragmentation of Ceo are still rather uncertain [12] (see Chap. 3.2.1), 
however, in any case it is safe to state that the fullerene ions emerging from 
electron impact ionisation (and laser desorption) sources are "hot" with 
internal energies on the order of tens of eV. Not only the magnitude but also 
the energy distribution is uncertain and is probably very broad (greater than 
thermal) for most of the ion sources used. 

The situation is more promising for negatively charged fullerene ions. 
Recently, an electrospray ion source, shown in Fig. 2.3, has been developed 
that can produce fairly intense beams of Ceo as long as a suitable solvent 
mixture is used [13]. The advantage of this kind of source is that the internal 
energy of the ions is much less than with conventional sources and is also 
better defined since the ionisation method is much more gentle. Electrospray 
fullerene sources are rather new in the collision community and have not yet 
been extensively used for collision experiments. 



Lenses 




Figure 2-3. Schematic diagram of an electrospray ion source, adapted from [13]. 
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For many of the experiments reported later in the book, it is necessary to 
have a well-defined target of neutral fullerenes. In order to determine cross 
sections for reactions it is also necessary to know the density of the neutral 
fullerene target. The fullerene vapour pressure has been determined as a 
function of temperature by a number of groups. Unfortunately, the values are 
not in agreement and are in some cases out by orders of magnitude. Some of 
the available data is plotted in Fig. 2.4 with the fitted parameters and details 
of the measurements given in Table 2.2. The corresponding target density is 
plotted as a function of temperature in Fig. 2.5. 




Figure 2-4. Log of vapour pressure as a function of inverse temperature. Straight lines: Data 
from Abrefah et at. [14] Their values for C 70 are also shown; Dashed lines: Data from 
Piacente et al. [15], the upper line is obtained on the first sublimation, the lower line for a 
subsequent sublimation; Dotted line: from Mathews et al. [16]; Dash-dot from Jaensch and 
Kamke [17]; Dash-dot-dot from Popovic et al. [18]. 
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Table 2-2. Comparison of Vapour Pressure Measurements. 







log(P/kPa) = A 


- B/(T/K) 




Method 


Temperature 
Range /K 


A 


B 


Reference 


Knudsen Effusion 


673-873 


6.6 + 0.1 


8281 +873 


[14] 


Knudsen Effusion (C70) 


703-873 


8.2 + 0.3 


9848 + 228 


[14] 


Torsion Effusion (1st) 


730-990 


8.28 + 0.20 


9154+150 


[15] 


Torsion Effusion (2nd) 


730-990 


8.73 + 0.25 


9668 + 200 


[15] 


Mass Spec. 


600-800 


8.18 + 0.03 


9475 + 80 


[16] 


Mass Spec. 


680-790 


7.49 + 0.14 


8267 + 120 


[17] 


Knudsen Evaporation 


663-823 


8.02 + 0.19 


8801 + 145 


[18] 




Figure 2.5. Target density as a function of temperature for (with one set of data for C70). 

Same sources as Fig. 2.4. 

Although there is still some disagreement concerning the correct vapour 
pressure data, the collision community generally uses the values obtained by 
Abrefah et al. [14]. Normally, purification by sublimation is not carried out 
prior to collision experiments and the data obtained with the presence of 
residual solvents is most relevant for the determination of the target 
densities. A detailed discussion of the reasons for using the Abrefah data for 
analysing collision cross sections has been given by Schwartz et al. [19]. 
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These authors base their choice of this data partly on their theoretical 
understanding of single-electron capture from Ceo to highly charged ions (see 
Chapter 7). Their results are very consistent with the Abrefah data and not at 
all consistent with the Piacente, Popovic or Kamke values. The Piacente 
values have, however been shown to be consistent with UV-vis absorption 
measurements on pre-sublimed samples [20]. The Popovic and Kamke 
results may lie too high due to the complexity of the data analysis and the 
assumptions that have to be made concerning either ionisation cross sections 
and electron gain on the one hand [18], or details of the Ceo - Ceo interaction 
on the other [17]. All absolute cross section values discussed in later 
chapters have been obtained by using the Abrefah et al. vapour pressure 
results, unless otherwise stated. 

1.3.1 Determination of Cross Sections 

In all of the fullerene experiments which determine absolute reaction 
cross sections, a projectile beam is sent through a neutral target (normally a 
gas cell rather than a second beam) and the reaction products are determined 
as a function of collision energy (or scattering angle if differential cross 
sections are being determined). The intensity of the reaction product, Ipmduct, 
is related to the cross-section for that reaction, a, as follows: 

hroducr = hn(Tl ( 2 . 1 . 2 ) 

where Ig is the initial intensity of the projectile beam, n is the number density 
of the target gas and / is the length of the target gas cell. This expression 
holds in the target density range where single-collision events dominate and 
the target is not depleted. The absolute cross section for a given collision 
energy is determined by measuring the intensity of the reaction product as a 
function of the pressure (number density) in the target cell. For most 
experiments with rare gas targets this pressure can be determined directly 
with rather good accuracy. For fullerene targets, the temperature of the target 
cell, containing some fullerene powder, is varied and n determined with the 
help of the vapour pressure measurements discussed in the previous section. 
An example of typical data is shown in Fig. 2.6. The gradient in the linear, 
low-pressure range of the data, where single-collisions dominate, yields the 
cross-section. Deviation from straight-line behaviour occurs due to the 
increasing importance of multiple collisions in the target cell. Such 
pressure/temperature dependent measurements can be very time consuming 
and often this will be carried out at just two or three values of the collision 
energy. Relative determinations will then be carried out for the other 
collision energies at constant pressure/temperature. 
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Figure 2.6. Example of the determination of collision cross-sections. The reaction product ion 
intensity (in this case CnC from + C^o collisions) is determined as a function of target 
gas number density. The cross section is obtained from a least squares fit to the data in the 
low target density range where single-collision conditions are operating. Adapted from [21]. 

The detection efficiency of the projectile beam (normally positive ions) 
should also be considered, in particular if a wide range of masses is being 
detected and/or collisions over a wide collision energy range are being 
investigated. MicroChannel plates are often the ion detectors of choice. The 
detection efficiency in this case depends strongly on the velocity of the 
species being detected. Heavy fullerene ions will be detected with much less 
efficiency than faster, smaller fragment ions. This has to be investigated for 
each individual experimental set-up (different microchannel plates can show 
quite different behaviour) and the raw data corrected accordingly. An 
example of such a procedure is given by Campbell et al. [22] and applied to 
fullerene projectile collision experiments in [23]. An example of the effect 
on a fragmentation mass spectrum of correcting for the mass (or velocity) 
dependence of the detection efficiency is shown in Fig. 2.7. The problem 
becomes much less severe when the fullerene ions have velocities beyond 
ca. 5x10“* ms ' (or collision energies beyond ca. 10 keV for Ceo^) where the 
detection efficiencies converge. 





Relative Int. In /leo 
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Fig. 2. 7. Relative fragment intensities as measured directly in Cgo* + Ar experiments at a 
centre of mass collision energy of 85 eV (open bars) compared to the intensity corrected for 
different detection efficiences of the carbon clusters (filled bars). From [23]. 



2. EXAMPLES OE COLLISION EXPERIMENTS 



2.1 Thermal Collision Energies 

The majority of reactive collision experiments carried out at thermal 
collision energies have used the technique known as SIFT: Selected Ion 
Flow Tube. This technique is a fast flow tube / ion swarm method for the 
study of the reactions of ions (positive or negative) with atoms and 
molecules under thermalised conditions over a wide range of temperatures. It 
has been extensively used to study ion-molecule kinetics. In contrast to most 
of the techniques used to investigate fullerene collisions, this does not study 
the products of individual collisions but is used to determine reaction rate 
constants. In the SIFT apparatus the ions are created in an ion source which 
is external to the flow tube. In order to study the reactions of fullerene ions, 
C6o‘*^, the ions were produced by electron impact ionisation of vaporised 
material [24]. The ions are then extracted from the ion source, selected 
according to their mass-to-charge ratio using a quadrupole mass fdter, see 
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Turbo Rotary Lobe 




Diffusion Interfeice Diflusion 

Pump Pump 



Figure 2-8. Schematic diagram of a SIFT (selected ion flow tube) apparatus used to determine 
reaction rate constants involving fullerene cations. Q1 and Q2 are quadrupole mass filters. 

Figure 2.8, and injected into a flowing carrier gas (usually helium) via a 
small orifice. The ions undergo many collisions with the carrier gas 
(typically 10^) and are thermalised. To study a particular ion-molecule 
reaction, a reactant gas is introduced at a measured flow rate into the carrier 
gas at a position along the flow tube where thermalisation has occurred. By 
measuring the count rates of the reactant and product ions with the second 
quadrupole mass spectrometer and relating them to the number density of the 
reactant gas it is possible to determine the rate constant for that particular 
reaction. 

The reaction rate coefficient determined in such flow tube experiments 
(in units of cm^mof^s ') can be related to the cross section for the given 
reaction as k(v) = vct(v) with v the relative collision velocity. Assuming 
that the velocity distribution in the experiments is thermal with a 
temperature, T then the normalised distribution for the relative velocity in 
the centre-of-mass frame of reference is 

f(v)dv = 42rv^(/i/22Tk^Tf"exp(-/iv^/2kj,T)dv ( 2 . 2 . 1 ) 

where p is the reduced mass. The rate constant at temperature T is then 
given by 

k(T) = j f(v)a(v)dv. (2.2.2) 

If we assume that the cross section depends only on the collision velocity or 
relative collision energy, Ecm (and not on the details of the initial and final 
states) and make the further rather extreme assumption that the collision 
cross section can be given by the line-of-centres model, discussed in Chapter 
9 then 
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a = 7rd\\-EJE^^)foxEc^>E^ 

= Ofor£„<£„ 

with Eg the threshold energy for reaction. Eq. (2.2.2) then becomes [25] 

k{T) = Kd'^ exp(-£'o (2.2.4) 

This is just the standard "collision theory" rate constant with an Arrhenius- 
like temperature dependence. Note that a number of simplifying assumptions 
have gone into deriving this expression and normally the pre-exponential 
factor determined in experiments is not in good agreement with the simple 
derivation. Also, it is quite common in SIFT experiments of ion-atom or ion- 
molecule collisions to apply a small electric field to "pull" the ions through 
the neutral target. In this case the relative velocity distribution is different 
from eq. (2.2.1). 



2.2 Hyperthermal collision experiments 

Hyperthermal collision experiments can be divided into two different 
energy regimes: low and high. "Low" in this case refers to collisions where 
the collision velocity is lower than the orbital velocity of the electrons. In 
this case, the main excitation channel is expected to be vibrational excitation. 
At "high" collision energies the excitation is dominated by electronic 
transitions. This is described in more detail in the following Chapter. Most 
of the studies in the literature fall into the "low" energy category with 
collision energies on the order of a few keV. The same kind of experimental 
set-up is used for both categories of hyperthermal collision. For these 
collisions fullerenes can either be the projectile ions (positively or negatively 
charged) colliding with a neutral gas target or the neutral targets in collisions 
with singly or multiply charged atomic ions or both projectile and target. 
Figure 2.9 shows an example of an experimental set-up used to study low 
energy collisions of fullerene ions with neutral targets, including neutral 
fullerenes [21]. This is one of the few set-ups used to study fullerene 
collisions where it is also possible to determine the angular distributions of 
the collision products. The fullerene ions are produced by electron impact 
ionisation and are then accelerated with a pulsed electric field into a small 
time-of-flight mass spectrometer. The ions are focused in time according to 
their mass to charge ratio [26] and can be mass-selected by briefly switching 
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off an applied electric field when the ions of choice are present at the "mass 
selector". The mass-selected ions are then steered and focused into a small 
target chamber where the collisions with the neutral target gas take place 
under single-collision conditions. The collision energy is controlled either by 
changing the acceleration voltage in the source region or by changing the 
potential of the target chamber and can be varied from 50 eV to 5 keV in the 
laboratory frame of reference. The energy resolution is typically on the order 
of 5% of the laboratory collision energy. Charged products from the 
collision are detected in a reflectron time-of-flight mass spectrometer [27] 
with an additional retarding field energy detector. The reflectron detector 
assembly can also be rotated around the target using a specially constructed 
target oven with a horizontal exit slit. This has been used to determine 
differential cross sections in fullerene cation - neutral fullerene collisions 
[28]. 

The apparatus has also been used to study positively charged products of 
hyperthermal collisions of neutral fullerenes with neutral gas targets [29]. In 
this case, a laser desorption source was used to produce an intense beam 
consisting of both positively and negatively charged fullerenes. The 
negatively charged fullerenes were accelerated to the desired laboratory 
collision energy and mass-selected. After mass selection a significant 
proportion of the negatively charged ions undergoes thermal electron 
emission due to the high internal vibrational excitation they receive from the 
laser excitation. Just before entering the collision cell, the remaining 
negative ions are removed from the beam with an electric field and only the 
accelerated neutral fullerenes enter the target chamber. 

Instead of time-of-flight mass spectrometers, it is quite common to use 
electric and magnetic field sector mass spectrometers for ion selection before 
and after collision. This gives the advantage of higher mass and energy 
resolution but they tend to be less flexible and it is more difficult to scan 
parameters such as the collision energy over a large range. One nice example 
of a set-up based on electric and magnetic sector mass spectrometers is that 
developed by the Schwarz group [30] and used mainly to study collision 
induced dissociation of molecular cations (including fullerenes). This is the 
apparatus that was used in experiments that first showed the possibility of 
collisionally capturing an atom inside the fullerene cage to produce 
endohedral fullerenes (see Chapter 8). A schematic diagram of the apparatus 
is shown in Fig. 2.10. It is optimised to work at a laboratory collision energy 
of 5 keV. The big advantage of the apparatus is the possibility of carrying 
out successive collisions followed by mass/energy selection. With the 
constellation shown in Fig. 2.10 it is e.g. possible to carry out the following 
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Figure 2.9. Schematic diagram of SOCCER apparatus (Scattering of Carbon Clusters Energy 
Resolved) used to study collisions between positively charged fullerene ions and neutral 
fullerene targets. From [21], 

experiment. Fullerene ions are produced by electron impact ionisation in the 
source region and then selected according to their mass/charge ratio. They 
then collide with a neutral target gas, e.g. He and a small fraction of the 
fullerene ions capture the He atom inside their cage. These ions are then 
selected by sector 2, collide in a second target chamber where they can 
undergo a neutralisation reaction. After this collision chamber any residual 
positive ions are removed from the beam by applying an electric field and 
the fast neutral molecules enter the third collision chamber where they can 
again undergo a charge transfer reaction to emerge as cations. These undergo 
a final mass-selection before being detected. Such a complicated and 
challenging series of collisions was carried out as described [31] and was 
one of the early strong indications of the endohedral nature of the ions 
produced in the first collisions. 

Many hyperthermal energy collisions have been carried out using neutral 
fullerenes as the target. Often, highly charged ions are used as the projectiles 
and the cross sections for electron transfer from the fullerene to the projectile 
are studied (see Chapter 7). Examples of these kinds of collision experiments 
are shown in Figs. 2.11 and 2.12. In Fig. 2.11 the charge state of the 
projectile ion after collision as well as any kinetic energy gain or loss is 
determined as a function of collision energy and initial charge state. After 
collision the projectile ions are selected according to their energy and charge 
state by means of a 180° cylindrical energy analyser. Alternatively, the final 
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charge state and energy has been detected by using deflector plates. The 
latter method also allows the possibility of measuring the cross sections for 
projectile neutralisation since the projectile particles have sufficient kinetic 
energy to be detected without additional acceleration at the detector. By 
using a position sensitive detector it is also possible to determine angular 
differential cross sections [32]. 




Figure 2-10. Schematic diagram of Schwarz experiment. Cl-7, collision cells; Dl-4, ion 
detectors. Adapted from [33]. 







Energy 

Analyser 



Figure 2.11. Apparatus used to measure total electron capture cross-sections in collisions 
between multiply charged ions, A'*'^, and Cgo- Adapted from [19] 



In the second set-up (Fig. 2.12), the ionised target fullerenes are detected 
using a linear time-of-flight mass spectrometer. This technique has the 
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advantage of simultaneously detecting the parent target ions as well as their 
fragment ions thus giving an indication of the amount of energy transferred 




Ceo jet 



to internal degrees of freedom in the collisions. 

Figure 2-12. Apparatus for studying charge transfer cross sections in collisions between 
multiply charged ions (including 060“*^) and neutral Ceo- Product ions and their fragments are 
detected with a time-of-flight mass spectrometer. Adapted from [34]. 



2.3 Coincidence experiments 

In all of the experiments discussed so far it is only the positively charged 
products of the collisions that are detected. This is of course a disadvantage 
when it comes to determining collision mechanisms and obtaining 
information on energy transfer etc. since much information is "hidden" in the 
neutral fragment distributions or emitted electron kinetic energies, for 
example. Much more detailed information can be obtained from coincidence 
experiments and a few coincidence measurements have been carried out in 
recent years to study the products of fullerene collisions. 

An early example is the set-up used by Vandenbosch and co-workers, 
shown in Fig. 2.13 [35]. In this experiment, an energetic beam of mass- 
selected negatively charged Ceo is produced and collides with a light target 
atom or molecule (e.g. H 2 ). Due to the kinematics of the collision, the 
fullerene fragments continue in a narrow cone about the original projectile 
beam axis. The positively charged fragment ions produced in the collision 
are then detected in coincidence. This allows the size distribution of the 
partner fragment ions to a particular mass fragment to be determined. The 
more energetic heavier fragments pass through a hole in the first detector 
and enter the second detector. The heavy fragments detected in this way are 
then measured in coincidence with the different light fragments detected 
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with detector 1. This experiment provided evidence for a binary break-up of 
the fullerene cage followed by sequential evaporation of C 2 molecules from 
the hot fragment ions produced in this way (see Chapter 5). 
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Figure 2.13. Apparatus used to measure coincidences between fragment ions produced in 
energetic collisions between C^o anions and light target gases. Adapted from [35]. 

The apparatus shown in Fig. 2.14 is quite different from that of Fig. 2.13 
but is used to obtain similar information [36]. In this case, fast protons 
collide with a neutral fullerene target. Electrons produced in the energetic 
collisions are used as the start signal for a time and position sensitive 
detector. In this way all positively charged fragments are detected in 
coincidence with one of the emitted electrons. 

The kind of set-ups in Figs. 2.11 and 2.12 have been combined to 
determine the target charge and fragmentation distributions in coincidence 
with a given charge-exchanged projectile ion, Fig. 2.15. Information on the 
kinetic energy distributions of the electrons emitted in collisions between 
highly charged ions and neutral fullerenes has been obtained by introducing 
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an electron energy analyser and measuring the energy distributions in 
coincidence with time-of-flight mass spectra (similar to the set-up in Fig. 
2.15) [37]. 




Figure 2.14. Apparatus for determining coincidences between fragment ions produced in 
collisions between energetic protons and neutral Ceo- Adapted from [36]. 
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Figure 2.15. Apparatus to measure the energy distributions of the charge-exchanged projectile 
ions in coincidence with the fullerene fragmentation spectrum. Adapted from [37]. 



2.4 Electron impact 



The collision energy dependence of electron impact ionisation and 
fragmentation of fullerenes has been studied extensively. Many of the results 
involving energetic electron collisions with neutral fullerenes have been 
obtained in the Mark group using the apparatus shown in Fig. 2.16 [9]. This 
is a double-focusing two-sector field mass spectrometer. An effusive 
fullerene beam interacts at 90° with an electron beam. For cross section 
measurements the electron current is reduced to ca. 10 pA to ensure single- 
collision conditions. Fligher currents can be used to produce significant 
amounts of multiply-charged ions in higher order reactions. The electrons 
are guided by a weak magnetic field and have kinetic energies in the range 
1-1000 eV with an energy spread of about 0.5 eV. The positive ions 
produced in the collision are extracted at right angles and enter the two- 
sector mass spectrometer through the slit SI. Metastable decay reactions 
occurring within the second field-free region of the mass spectrometer can 
be studied by using the so-called MIKE scan technique [38]. Measurement 
of the metastable fragment peak shape, obtained by scanning the electric 
field in the second sector can be directly related to the kinetic energy 
released during the metastable fragmentation. This, in turn, can provide 
information on dissociation energies although it is dependent on the model 
assumed for the fragmentation mechanism. This has led to problems 
concerning the determination of the dissociation energies for fullerenes (see 
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Figure 2.16. Two-sector mass spectrometer used extensively to study the metastable 
fragmentation of fullerene ions produced by energetic electron impact. Fullerene ions are 
mass-selected by the magnetic sector and the fragmentation that occurs in the second field- 
free region (ff2) is studied by scanning the electric field in the second sector of the machine, 

Adapted from [39]. 

Chap. 3.2.1). Recently the apparatus has been upgraded to a three-sector 
instrument, which has led to a substantial improvement in the achievable 
resolution [39]. 

The same apparatus has also been used to study electron collisions with 
positively charged fullerene ions by inserting an intense electron beam 
between the two sectors [40]. However, most of the literature results 
concerning cross sections for electron impact ionisation and fragmentation of 
fullerene ions have been obtained with the apparatus developed by the 
Salzborn group and shown in Fig. 2.17 [41]. The fullerene ions were 
produced in an ECR ion source (see 2.1.2) and then crossed an intense 
electron beam with energy in the range 10-1000 eV. After the interaction the 
product ions were separated by a 90° magnet and detected with a single- 
particle detector. Absolute cross sections were obtained by using the 
dynamic crossed beam method where the electron beam is moved through 
the ion beam with simultaneous registration of the parent and product ion 
beam intensities [42]. 




Figure 2-17. Experimental set-up for determining absolute cross-sections in ion-electron 
collisions. Fullerene ions are produced in an ECR source, mass and energy selected and then 
collide with a high power electron gun. Both the fragment ion distribution and the parent ion 
beam intensity is recorded. Adapted from [41]. 
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2.5 Surface Collisions 




Figure 2.18. Apparatus used for studying the collisions of fullerene ions on surfaces. Adapted 

from [43]. 



Fullerenes colliding with surfaces at relatively high impact energies of a 
few hundred eV and higher were shown to be very resilient with a tendency 
to "bounce" away from the surface. This led to a number of studies of the 
dynamics of fullerene surface collisions, discussed in detail in Chapter 10. A 
typical surface scattering apparatus is shown in Fig. 2.18 [44]. The set-up is 
quite similar to the gas phase scattering experiment shown in Fig. 2.9. In the 
surface scattering experiment, laser desorption is used as the fullerene ion 
source. The fullerenes are directly desorbed into a pulsed expansion of a rare 
gas. Fullerene cations, anions and neutrals are cooled and carried by the rare 
gas pulse through a skimmer and into an acceleration region. The ions are 
accelerated perpendicular to the direction of expansion of the gas pulse and 
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mass selected with a pulsed electric field. A pulsed deceleration lens is used 
to decelerate the fullerene ions to the desired collision energy thus allowing 
the surface region to remain field-free. After scattering from the surface the 
ions are pulsed accelerated and detected with a rotatable time-of-flight mass 
spectrometer. The velocity distribution of the scattered ions can be mapped 
out by changing the time delay between ion impact on the surface and 
pulsing on of the extraction acceleration field. The scattering dynamics can 
be strongly influenced by impurities on the target surface. For this reason it 
is essential to clean the surface by sputtering and to control the surface 
quality by FEED and Auger spectroscopy. 

Fullerenes have also been used to study the effects of implanting high 
energy (MeV) clusters into solids [45]. For these studies the fullerene ions 
were produced in a tandem accelerator where a Cs ion gun is used to 
bombard a Ceo target in order to produce negatively charged ions with a few 
keV energy. These ions are then accelerated to MeV energies and collide 
with N 2 gas to produce positive ions that are then further accelerated in the 
second part of the tandem. 
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Chapter 3 

THEORETICAL MODELS 



In this chapter we will review some basic collision dynamics and 
convenient frames of reference. More detailed discussions can be found in 
text books dealing with atomic and molecular collisions e.g. [1-4]. We will 
also discuss some of the theoretical models used to describe the general 
dynamical behaviour of fullerenes. Theoretical considerations that are 
specific to a certain sub-set of collision experiments will be covered in the 
relevant chapter. 



1. FRAMES OF REFERENCE AND CLASSICAL 
MECHANICS OF COLLISIONS 

Two frames of reference are normally used to describe collision 
experiments: the laboratory and the centre-of-mass frames. The laboratory 
frame of reference is the one most convenient for experimentalists. It 
describes the experimental results in terms of the kinetic energy or velocity 
of the projectile ion beam and, if appropriate, the scattering angle of the 
projectile beam as directly measured in the laboratory. This is a convenient 
frame of reference to use at high collision energies where the velocity of the 
projectile with respect to the electron orbital velocities in the target atom or 
molecule becomes a relevant factor (see below). However, at low collision 
energies it is normally more convenient to think in terms of the centre-of- 
mass reference frame. The force between the particles depends only on the 
relative distance between them. For this reason it is of advantage to describe 
the collision in terms of the relative separation of the colliding partners and 
the motion of their centre-of-mass rather than the motion of each individual 
partner. If no external forces are acting then the motion of the centre-of-mass 
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is unchanged during the collision. It can therefore be extracted from the total 
energy and we only need to consider the remainder, which is the kinetic 
energy of the relative motion. The relative energy of the collision partners is 
important since this determines the maximum energy available for 
conversion to internal excitation in the collision. In a fully inelastic collision 
the entire centre-of-mass collision energy is converted to internal excitation 
and the relative velocity of the collision partners after collision is zero (e.g. 
in the experiments discussed in Chapters 8 and 9). The velocity of collision 
partner i in the laboratory frame can be written as 

v,=v,„ + w, (3.1.1) 

where is the velocity of the centre of mass 

and w, is the velocity of particle i with respect to the centre-of-mass. 
Conservation of linear momentum implies that 

m,Wi + m2W2=0 (3.1.3) 

where ntj and m 2 are the masses of the collision partners. The relative 
velocity is given by 



(3.1.2) 




V = Wj -W2 



(3.1.4) 



and the relative kinetic energy in the centre-of-mass frame is 



1 2 1 2 

2 + 2 ^2^2 




(3.1.5) 



where p is the reduced mass 



M = 



m/iij 



(3.1.6) 
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For hyperthermal collision experiments, where the thermal velocity of the 
target atom/molecule can normally be neglected compared to the high 
velocity of the projectile, it is convenient to give an expression for Ecm in 
terms of the laboratory collision energy, Eiai,: 



1 , ,2 

2 (m, + m^) 

1 m.m^ 9 

2 (m^ + ^2 ) (mj + ^2 ) 



where the subscript 1 refers to the fast projectile and 2 to the slow target. 

Equation (3.1.7) implies that the amount of energy available for internal 
excitation in a hyperthermal atom -fullerene collision will depend very much 
on whether the fullerene is the projectile or the target. If a fast fullerene ion 
collides with a slow He target the centre-of-mass collision energy is only 
0.55% of the laboratory kinetic energy whereas if a fast He collides with a 
thermal fullerene target then the centre of mass collision energy is 99.4% of 
the laboratory kinetic energy. 

The conversion between scattering angles and scattering intensities in the 
lab and centre-of-mass frames of references is most conveniently made with 
the help of Newton diagrams. Fig. 3.1. In the laboratory system, the centre- 
of-mass moves with a constant velocity Vcm, which is unchanged by the 
collision. In the centre-of-mass system, the centre-of-mass is at rest and the 
deflection dcm is the deflection of the relative velocity vector v'. In elastic 
collisions, the magnitude of the relative velocity does not change. The final 
velocity of a particle in the laboratory frame of reference is given by a vector 
addition of its final relative velocity with respect to the centre-of-mass and 
the velocity of the centre-of-mass itself. The particle speeds with respect to 
the centre-of-m a ss also do not change 



= m^v'l (mj + m 2 ) 
w'2 = my/ (mj + m2 ) 



However, the final laboratory velocities are different from their initial 
values in both magnitude and direction 

< = Ve„ + w'i 

v'2 = + w'2 



(3.1.9) 
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Figure 3.1. Newton diagrams illustrating the conversion between laboratory and centre-of- 
mass coordinates for collision processes, (a) elastic collision between two particles with initial 
velocities Vi and V 2 impinging at 90° and with equal masses, (b) as (a) but with particle 2 
initially at rest, (c) Inelastic collision between a fast Qo projectile and a stationary Xe atom 
where half of the centre-of-mass energy is converted to internal excitation, (d) the same 
centre-of-mass collision energy, energy loss and centre-of-mass projectile scattering angle but 
this time for a fast Xe projectile colliding with a stationary C^o. 



Fig. 3.1a gives an example of an elastic collision between two particles 
intersecting at 90° with initial velocities Vi, V 2 and masses ni] = m 2 . Particle 1 
is scattered through an angle in the laboratory frame of reference. This 
corresponds to a centre-of-mass scattering angle of 0cm- We now need to 
convert the centre-of-mass scattering intensity (which is what is normally 
calculated from model potentials) Icm( 0cm) to the corresponding laboratory 
intensity habi &iah)- For each 0cm, conservation of intensity gives 

hab {Qlab)d^lab = hm {^cm)dO) (3.1.10) 

where di2iab is the solid angle element in the laboratory system 
corresponding to dcOcm in the centre-of-mass system. The "Jacobian", or solid 
angle ratio of the transformation dcOcm/dC2iah can be read-off from Fig. 3.1a to 
be 
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dCOcn, _ vf 

d^iab w'l'lcos^ 



(3.1.11) 



If more than one centre-of-mass angle corresponds to a given laboratory 
scattering angle this procedure is repeated for each dcm and the centre-of- 
mass intensities are summed to yield the total observable laboratory intensity 

3fab( ^lab)- 

Further examples of Newton diagrams for elastic and inelastic scattering 
are given in Fig. 3.1. 

1.1 Collisional Energy Transfer 

The simplest possible situation is a collision between two structureless 
bodies that interact via a central conservative force. In this case, no energy 
can be exchanged between the relative motion of the particles and any 
internal modes. Collisions of this type are called elastic. The only outcome 
of the collision is that the particles deflect one another from their original 
straight-line trajectories, as discussed above. Inelastic collisions, on the other 
hand, involve the transfer of energy from or to internal modes while in 
reactive collisions chemical bonds are broken or formed. At low energies, 
collisions between rare gas atoms are elastic but for molecules the collisions 
may induce changes in internal energy or cause chemical rearrangement. In 
order to fully describe the excitation mechanisms and collision outcome one 
needs to have accurate potential energy surfaces and coupling constants, 
which is unrealistic at present for a complex molecule as large as Ceo- 
Instead, much simpler considerations are generally made to model the 
energy transferred in the collisions. 

One extreme simplification which seems to work rather well for low 
energy atom-fullerene collisions is that the atomic collision partner only 
"sees" a small part of the fullerene cage and undergoes an elastic collision 
with this'. The kinetic energy transferred to the part of the cage involved in 
the collision is then subsequently transferred to kinetic energy of the other 
atoms in the cage corresponding to excitation of many vibrational modes 
(normally a fast equilibration of energy among the vibrational degrees of 
freedom is assumed). This is known as a "spectator" model [2]. The idea is 
illustrated in Fig. 3.2 for an atom-fullerene collision where the fullerene is 
initially at rest. The atom. A, collides with just a small part of the cage, B - 



* For convenience in the discussion that follows we will normally regard the atom (or atomic 
ion) collision partner as being the projectile and the fullerene molecule as the target. Of 
course, exactly the same arguments hold for the opposite case: fast fullerene ions colliding 
with stationary atomic targets. 
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for Ceo this could reasonably be assumed to be e.g. a carbon atom or a C 2 
molecular unit - and the rest of the carbon cage, C, is simply a spectator to 
the collision. Any change in the velocity of B due to an elastic collision with 
A is thus a change in the relative velocity with respect to C and hence a 
change in the vibrational energy of the fullerene. The initial momentum Pa is 
partitioned after the elastic collision between A and B so as to conserve 
energy and momentum. The additional assumption is made that the collision 
is collinear i.e. we consider only a head-on collision between A and B. The 
velocity of B after the collision is then 

v/ = 2pA/(«tA + «tB) (3.1.12) 



The vibrational excitation is just given by the change in the relative 
velocity of B and C 




1 muffle 

2 (m^ -h -h trig f 



(3.1.13) 



The ratio of vibrational energy transfer to initial centre-of-mass collision 
energy is thus 

^vib _ 1 4p>gfflc 2m^ {m^ + ) 

^cm ^ {m^ + mgf {mg+m^) pl{mg + m^) (3.1.14) 

4m^mgm^ (^m^ + mg+ ) 

{mA + rngf {nig +m^f 



This expression is plotted in Fig. 3.3 as a function of mg for three different 
atoms colliding with Ceo'. He, Ar and Xe. The energy transfer efficiency is 
always highest when the unit of the carbon cage taking part in the collision 
has the same mass as the collision partner. 




3/ 3V 
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Figure 3.2. Spectator Model 




Figure 3.3. Energy loss for Ceo-rare gas collisions, from the spectator model. 



The spectator model regards parts B and C of the molecule to be 
connected by a very loose spring. So, when the colliding atom A interacts 
with B, eq. (3.1.10) gives an upper limit for the fractional energy transfer. In 
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physical terms, the limit of a weak oscillator spring is when the duration of 
the collision, to is short compared to the oscillation period, t^ [2]. We can 
estimate the duration of the collision as 

t^=a/Y ( 3 - 1 - 15 ) 

where a is a typical interaction distance and v is the relative velocity during 
the collision. The spectator model can be expected to be valid for relative 
collision velocities 



Y>a/ty=av (3.1.16) 

where v is the oscillation frequency. 

If we take the average Ceo vibration frequency, v = 2.79x10^^ s ', and the 
diameter of the fullerene cage to estimate the interaction distance, a = 10 A, 
then we obtain v > 27900 mi ' . For collisions in which energetic Ceo ions 
collide with stationary targets this corresponds to a laboratory kinetic energy 
£i«,s2900eV. 

At lower collision velocities the duration of the collision is long 
compared to the oscillation period. In this case the oscillation can "adjust" to 
the perturbation caused by the collision or, in other words, the collision 
partner "sees" the entire molecule. In this limit, the energy transfer efficiency 
is very low. 

Evidence that this simple model provides a rather good picture of the low 
energy hyperthermal collisions between fullerenes and small atomic targets 
was obtained in experiments in which a Ceo^ projectile beam was collided 
with a stationary He target [5]. The kinetic energy lost by the projectile ion 
was measured in the experiments and related to the effective mass on the 
fullerene cage, m^, according to eq. (3.1.13). The results are plotted in Fig. 
3.4. At low centre-of-mass energies, corresponding to relative velocities 
V < flv , the atom "sees" the entire fullerene cage whereas for velocities the 
atom only interacts with a small part of the fullerene cage, in this case a C 2 
unit. The dotted line on Fig. 3.4 marks the centre-of-mass energy for which 
the relative velocity is 27900 ms * and thus marks the boundary between the 
high and low velocity limits, according to eq. (3.1.16). The agreement with 
the experimental results is surprisingly good. 
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Figure 3.4. Experimentally determined effective mass on the fullerene cage that interacts with 
a colliding He atom, plotted as a function of centre-of-mass collision energy [5] The dotted 
line marks the separation between the low velocity regime and the high velocity regime where 
the spectator model is expected to give a good approximation. 

This simple model is fairly good for velocities beyond those given by eq. 
(3.1.16) but not indefinitely. As the interaction time during the collision 
becomes shorter and shorter the energy transfer efficiency due to elastic 
collisions with individual atoms or small molecular units decreases again. 
For the high energy limit, where the collision duration is comparable to the 
period of the electron orbits in the target then the energy is not efficiently 
transferred to the nuclear motion of the target (vibrational excitation) but 
instead there is predominantly a direct electronic excitation. This is often 
referred to as the crossover from the "nuclear stopping" regime to that of 
"electronic stopping". This terminology is taken from the ion implantation 
community and, to a reasonably good approximation, one can treat fullerene 
target molecules as a very small double layer of graphite and use the models 
developed for describing energetic ion implantation of materials to estimate 
the extent of energy transfer. 

The collisions in the "low energy" limit, i.e. where nuclear stopping 
dominates, are treated in a fashion similar to that discussed above for the 
simple spectator model. The collisions are treated as binary elastic collisions 
between the atomic projectile and individual atoms in the target. The main 
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difference is that the collisions are not restricted to collinear, zero impact 
parameter collisions but the scattering due to the interaction potential of the 
projectile and target is also taken into consideration and multiple binary 
collisions can occur during passage of the projectile atom through the target 
molecule. The nuclear stopping is most simply and generally described by a 
suitable scaling of energy and cross-section as described in detail by 
Lindhard et al. [6]. Here, the dimensionless quantities 



£ = E 



aM^ 



(3.1.17) 



and 



p = RNM^An:a^ 






(3.1.18) 



are measures of collision energy and range. Zj, Mj and Z 2 , M 2 are the atomic 
numbers and mass numbers of the projectile ion and target atom, 
respectively, E is the projectile ion kinetic energy. R is the range of the 
interaction (dependent upon the interaction potential), N is the number of 
atoms per unit volume and a = 0.8853a„Z”'" , a characteristic screening length 
where = Zj^" -i-zj” and ao is the Bohr radius^. The nuclear stopping power 
is defined (in atomic units) as 



^ _ d£ AjtaZ^Z^e ^ 

dp (Mj-I-Mj 

The derivative d£/dpi& a dimensionless measure of the reduced stopping 
power that allows all nuclear stopping cross sections to be described by one 
curve. Fig. 3.5 shows the computed reduced stopping power for two 
common forms of the interaction potential; a Thomas-Fermi type potential 
(dashed line) and a screened Bohr potential (dotted line). Details of the 
calculation can be found in Lindhard et al.[6]. 

The ion atom potential must have the form 



^ Atomic units are used throughout this section unless otherwise stated. 
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V(R) = S^j, 

R (3.1.20) 

where the function u must tend to unity when /? — > 0 , and must vanish when 
R ^oo ^ because the Coulomb potential is screened at large distances. The 
screened Bohr potential assumes a power law potential, V(R) R~" or 

(3.1.21) 



with s = 2 and k 2 = 2/(0. 885 3e). The slightly more realistic simplified 
Thomas-Fermi potential describes the static interaction between two atoms 
as u-u(Rla). The classical deflection angle, and thus the elastic 
scattering cross section, can be shown in this case to depend only on the 
reduced parameters e and p [6]. The reduced stopping power calculated 
from the Thomas-Fermi potential. Fig. 3.5, clearly shows a maximum in the 
nuclear stopping with a rapid fall in energy transfer efficiency due to elastic 
scattering as the reduced collision energy increases beyond this maximum. 
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Figure 3.5. Reduced stopping power as a function of the square root of the reduced collision 
energy (proportional to projectile velocity). Dashed: Thomas-Fermi interaction potential, 
Dotted: screened Bohr potential, Full: electronic stopping. After Lindhard et al. [6]. 
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The electronic stopping is approximately proportional to the projectile 
velocity for velocities that are low compared to the orbital velocities of 
electrons which could be carried by the ion, i.e. according to the Thomas- 
Fermi estimate [ 6 ] 



v<Vi=Zf'\ (3.1.22) 

where v^=e^/h. In this range, the electronic stopping power can be 
approximated by 



de 

dp 



= ke 



1/2 



(3.1.23) 



with ^ = 0.15 [7]. This has also been plotted on Fig. 3.5 (full line). Electronic 
and nuclear stopping are of the same order of magnitude when Vj » 0 . 1 v„. or 
f = 2 . For higher projectile velocities the electronic stopping rises sharply 
to a maximum and then falls off again as 1 /v^ [ 8 ]. 

Let us now consider some absolute values for the energy transfer in 
collisions between projectile ions and Ceo. The stopping can be converted to 
the energy transfer to internal modes of the fullerene in eV by making some 
additional assumptions. We regard only central collisions and regard the Ceo 
as a thin carbon foil with a density 

^ = 2 60 ■ 12m ^0 15 (3.1.24) 

Ajtr 

where r = 3.55 A is the radius of the fullerene cage and the factor 2 takes 
account of collisions at both the front and the back of the cage. The energy 
lost by the projectile ion in eV (assuming interaction occurs with individual 
C atoms in the cage) is then 

AE" = D^pg ! cm^^S^eV /(//g/cm^)] . (3.1.25) 



This is plotted as a function of projectile kinetic energy in Fig. 3.6 for H"^, 
He"^, Ar"^ and Xe"^ projectile ions colliding with Ceo, using a simplified 
Thomas-Fermi potential to describe the nuclear stopping. 

This simple picture is not successful at very low collision energies where 
the Thomas-Fermi approximations are not expected to work very well (e.g. 
predicting energy transfer that is significantly higher than the collision 
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energy!). In this region, the spectator model is a better approximation, as 
discussed above. For comparison, we show the predictions of the spectator 
model, assuming collisions with a single carbon atom on the cage, for the 
He"^ collisions in Fig. 3.6. However, at energies around or beyond the 
maximum in the nuclear stopping power the results using the Thomas-Fermi 
approximation should give reasonable values. They are also in very good 
agreement with somewhat more sophisticated although similar calculations 
using the TRIM programme (Transmission and Ranges of Ions in Matter [9]) 
[10]. Note that the energy transfer efficiency increases strongly with the 
mass of the projectile ion. Electronic stopping dominates already at 
relatively low laboratory energies for collisions with protons and helium 
ions. This also has consequences for the fragmentation patterns observed in 
experiments as will be discussed in Chap. 5. 

Some more sophisticated calculations have been carried out e.g. taking 
into account the impact parameter dependence [11], using a Monte Carlo 
method to estimate the nuclear stopping [12] or using an improved treatment 
of the electronic stopping [13] but the overall picture remains close to that 
estimated from the Lindhard model, discussed above. 




Projectile Kinetic Energy, [eV] Projectile Kinetic Energy. E^, [eV] 



Figure 3.6. Energy transfer to the fullerene cage in collisions with energetic projectile ions, as 
a function of ion kinetic energy. Dashed line: nuclear stopping calculated using a simplified 
Thomas-Fermi potential [6], full line: electronic stopping (eq. 3.1.23). The estimate using the 
spectator model for He+ colliding with an individual C atom on the fullerene cage, (eq. 
3.1.14) is shown for comparison (dotted line). Note the different scales on the left and right 

sides of the figure 
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1.2 Scattering Dynamics 

There have been relatively few studies of the scattering dynamics of 
fullerenes so far. For hyperthermal collision energies it is possible to make 
use of simple classical models of the interaction behaviour for heavy particle 
collisions (i.e. atom-fullerene, ion-fullerene, molecule-fullerene or fullerene- 
fullerene). 

Elastic collisions can be described by the following relations. For a fixed 
centre-of-mass collision energy E, interaction potential U(r) and impact 
parameter b, the centre-of-mass deflection angle in an elastic collision is 
given by the classical deflection function [1,2] 



Z{E,b) = 7T-2bj 

^0 



drlr^ 



(3.1.26) 



where Uejj is the effective interaction potential given by the sum of the 
interaction potential and the centrifugal energy and r is the separation of the 
colliding particles 







(3.1.27) 



The expression for the differential cross section can be obtained from the 
deflection function and is written as 






b 

{dz/db)smz 



(3.1.28) 



The expression for the total elastic scattering cross section is given by the 
integration of eq.(3.1.28) over all scattering angles 



^ total {E) = ^(j{Z,E)27t?,mzdz. 

0 



(3.1.29) 
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Thus by using a suitable interaction potential U(r) it is possible to 
calculate the cross sections and to compare them with those obtained from 
the experiments. Alternatively, if the experiments have sufficient resolution 
one would use the experimentally determined differential cross sections to 
determine the interaction potential. Unfortunately the experiments are not 
yet sufficiently advanced to allow us to do this. Also, there are many 
competing reaction channels in fullerene collisions that make the situation 
rather complicated [14]. 

Some estimates of interaction potentials have been made in the literature 
for fullerene-atom collisions, in particular with reference to Penning 
ionisation experiments [15] in order to estimate the expected magnitude of 
the cross sections and probability of energy transfer to compare with the 
experimental results. 

Presently the most reliable interaction potential was calculated by 
Girifalco for fullerene-fullerene collisions using a Lennard-Jones type 
potential [16, 17]. In this potential, the C-C interactions within the Ceo 
molecule are averaged over the surface of a sphere. This reduces each Ceo 
molecule to a pseudo-atom with just 3 centre-of-mass coordinates and 
introduces rotational isotropy. The potential has the form 



U (r) = -a 



- + - 



5(5-1) 5 ( 5 -tl) 






+A 



-+- 



-+- 



r( 5 -l)^ 5(5 + 1)^ 

(3.1.30) 



where the parameters are a= 0.04678 eV, fd = 8.486 x 10’^ eV and 5 = r/7.1 
A. The simplicity of the Girifalco potential is one of its advantages. 
Although it was shown that the Girifalco potential underestimates the 
stability of the liquid phase of Ceo [18, 19] it is relevant for the solid phase. It 
is also consistent with elastic scattering measurements for Ceo^ + Ceo gas 
phase experiments [14] when an additional polarisability term is added 



V^,{r) = 



T]e^ 

2 ? 



(3.1.31) 



with the polarisability of Ceo, p = 80 ± 4 A^ [20, 21]. 

The Girifalco interaction potential (eq. 3.1.30) is plotted in Fig. 3.7. The 
potential rises very rapidly as the fullerenes approach each other and is quite 
close to a hard sphere potential. Due to the rather unusual properties of the 
intermolecular potential it is not possible to use the conventional unique 
small angle scattering relationship between and b [22] even for relatively 
high collision energies of a few hundred eV. Instead, the deflection function 
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has to be calculated for each collision energy [14]. At the relatively high 
collision energies that have been studied so far (180-1500 eV) there is very 
little influence of the shallow potential well and the interaction probed in the 
collisions is purely repulsive. Lower energy collisions, although difficult to 
carry out experimentally, should be feasible and would provide a good and 
needed test of this simple but useful interaction potential. 




Figure 3.7. "Girifalco" potential [16, 17] for the Qq-Csq interaction 



2. FRAGMENTATION: STATISTICAL MODELS 

Generally the fragmentation behaviour of collisionally excited fullerenes 
can be reasonably well understood in terms of simple statistical models. In 
these models it is invariably assumed that the excitation energy is rapidly 
statistically distributed among all internal degrees of freedom of the 
molecule. This seems to be a good assumption for relatively low collision 
energies where the energy transfer is predominantly via vibrational 
excitation (nuclear stopping) and where there is sufficient time for energy 
redistribution to occur (> ps) before significant fragmentation takes place. 
However, the assumptions are more questionable for high energy collisions 
dominated by electronic stopping or where there is simply not sufficient time 
for equilibration to occur before the molecule disintegrates. Some reports of 





3. Theoretical Models 



45 



non-statistical fragmentation have been made and will be discussed in 
chapter 5, however, the evidence for non-statistical behaviour is still rather 
limited and in some cases the results can also be consistent with a statistical 
model. 

Here we present the simplest and most commonly used statistical models 
applied to the fragmentation of fullerenes. 



2.1 Arrhenius Fragmentation Rates and the Finite Heat 
Bath Correction 

A highly excited fullerene can release its energy in different ways; it can 
fragment, it can emit electrons or it can emit photons. One of the most 
interesting aspects of excited neutral fullerenes is that all three decay 
channels compete with each other for excitation energies and time scales 
typical of many experiments [23]. This is unusual for large molecules and is 
a consequence of the high dissociation energy combined with a relatively 
low ionisation potential as well as the ability of the molecule to "store" large 
amounts of energy for relatively long times due to the large number of 
vibrational degrees of freedom. For positively charged fullerenes the 
situation is simpler since statistical electron emission ("thermionic 
emission") does not compete significantly with fragmentation and radiative 
cooling due to the high value of the second ionisation potential (11.4 eV). 
The importance of competing statistical decay channels has not always been 
recognised and this has contributed to the difficulty of obtaining a reliable 
value for the dissociation energy for Ceo i-e. the energy required for the 
reaction Ceo — > C 58 - 1 - C 2 . Values for this energy reported in the literature 
vary over a large range [24] and there is still a certain amount of controversy 
and uncertainty concerning the correct value. A further problem with 
determining this value is that it cannot be measured directly in an 
experiment. Most experimental determinations involve invoking a statistical 
model and making some assumptions about the mechanism such as the 
looseness of the transition state [25, 26]. Until recently the experimental 
values were significantly below estimates obtained from high-level density 
functional calculations [27], however, the experimental values are creeping 
upwards and the values being reported are now compatible with the most 
reliable theoretical estimates [24, 25, 28-31]. In the following, where 
necessary, we will use values that are presently considered to be the most 
appropriate to describe the metastable fragmentation (or C 2 evaporation) of 
fullerenes, however, it should be noted that this may change as our 
knowledge of the complex mechanisms involved increases. The present 
status of the dissociation energies of the ions is summarised in Fig. 3.8. 
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Although the absolute values are still under investigation the relative values 
of the dissociation energies should be reliable. 




Figure 3.8. Dissociation energies for fullerene cations. Experimental: Circles from [30]. Up- 
triangles from [29]. Down triangle from [32]. Squares are relative values from [33] that have 
been normalised to the n=54 value of Tomita et al. [30]. Theory: Crosses from [31]. 



The dissociation energies of the neutral fullerenes are related to the 
fullerene cation dissociation energies according to the scheme shown in Fig. 
3.9. The dissociation energy for C 2 loss from Ceo is thus 0.54 eV higher than 
that from singly ionised Cgo- The ionisation potentials of the fullerenes are 
tabulated in Table 3.1. 
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Figure 3.9. Diagram illustrating the reasoning behind the deduction of the cation C 2 
dissociation energy, D^at, from that of the corresponding neutral fullerene value, D„eatrai- 





3. Theoretical Models 



47 



Table 3-1. Ionisation Potentials of the Fullerenes 



Fullerene 


IP [eV] 


Reference 


Fullerene 


IP [eV] 


Reference 


C48 


7.t7+0.t 


134] 


Cso" 


11.4+0.04 


135] 


C50 


7.61+0. tl 


134] 


c 


16.6+1.0 


136] 


C52 


7.07+0.1 


134] 




15.6+0.5 


137] 


C58 


7.07+0.1 


134] 




17.0+0.75 


138] 


Ceo 


7.59+0.02 


139] 


p 3 -h 
'-60 


20+1.0 


140] 


Cei 


7.07+0.1 


134] 


Cso"" 


22.3 


[41] 


Cee 


7.27+0.14 


134] 


Cro" 


11.4+0.06 


[35] 


Ce8 


7.27+0.14 


134] 


p 2+ 

'-70 


16.5+2.0 


[42] 


C70 


7.48+0.05 


139] 


p 3+ 
'-70 


19.2+5.0 


[42] 


C72 


6.97+0.1 


134] 


p 4+ 
'-70 


29.2+8.3 


[42] 


C74 


6.84+0.1 


134] 








C76 


7. 1+0.1 


139] 








C78 


7.05+0.1 


139] 








Cso 


6.84+0.1 


134] 








C84 


7.15+0.1 


139] 








CloO-200 


6.2-6.45 


134] 









The simplest description of the dissociation rates is given by the 
Arrhenius-like equation 






f-D^ 



k,T 



(3.2.1) 



where T is the internal temperature of the fullerene and D is the dissociation 
energy. A similar expression is used to describe the delayed "thermionic" 
emission of electrons with D given by the ionisation potential and Ab„s 
replaced by the appropriate pre-exponential factor for ionisation [43]. 

The appropriate parameters to use in order to describe the C 2 emission 
from Ceo neutrals, ions and other fullerene species are still a matter of 
investigation, as mentioned above. The most recent experimental 
determination of the dissociation energy for C 2 loss from neutral Ceo places it 
at 10.8±0.3 eV [32], in good agreement with theoretical values (ca. 11-11.5 
eV [27, 31]). This corresponds to a dissociation energy for the singly 
charged cation of ca. 10.3 eV. The present status of our knowledge 
concerning the dissociation energies of the singly charged fullerene cations 
is plotted in Fig. 3.8. Values for the dissociation energies of doubly charged 
cations can be found in [31] and references given there. Ceo has an 
exceptionally high dissociation energy, compared to the other fullerenes but. 
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even more unusual is the extremely high value of the pre-exponential factor 
Aoiss that is needed for Ceo and the other fullerenes in order to obtain 
dissociation rates that are compatible with experiment [25, 30]. The most 
recent evaluation and modelling of experimental data provides a pre- 
exponential factor of 3.4x10^' s’^ for C 2 loss from Ceo with somewhat lower 
values for the other fullerenes [32]. Theoretical estimates can obtain values 
within the range 10^°- 10^^ [25, 44] but for many researchers it is intuitively 
difficult to accept such high values and the fragmentation mechanism is still 
under debate. A simplified expression for the temperature dependence of the 
Ceo dissociation pre-exponential factor that is often used is [25] 

=3.0xl0'^r'(l-exp(-1441/r))%-' (3.2.2) 

The temperature commonly used in eq.(3.2.1) and (3.2.2) is not simply 
the vibrational temperature of the molecule but is a "heat bath" temperature. 
This is the temperature that a heat bath must be set to so that the canonical 
rate constant k(T) is equal to the micro-canonical rate constant k(E)\ The 
finite heat bath correction [45] relates the heat bath temperature to the 
vibrational temperature 



T,,{E) = T{E) 



, 1 D 1 

1 H 1 

2TC 12 




+ .... 



(3.2.3) 



where En, is the vibrational temperature of the fullerene, C is one less than 
the heat capacity of the emitting species at temperature T (in units of ks) and 
D is the activation energy of the reaction under consideration. 

It is generally accepted that the most significant fragmentation channel of 
an excited fullerene at modest excitation energies is the successive loss of C 2 
units (with the exception of fission reactions of highly charged fullerenes as 
discussed in Chapter 6). This is also the channel that is observed in 
experiments that study the metastable decay of fullerenes on the 
microsecond time scale. However, for high levels of excitation it is possible 
that the fullerene cage rapidly loses its structure (on the sub-ps time scale) 
[46] and it is then energetically more favourable for it to break up into ring 
and chain fragments. The cross-over from predominantly C 2 evaporation to 
rapid break-up into two or more larger fragments shows some signatures of a 

^ A micro-canonical ensemble is an ensemble of sub-systems where each sub-system has a 
fixed energy. A canonical ensemble is an ensemble of sub-systems where energy can flow 
between the sub-system and a heat reservoir so that the temperature is fixed but the energy 
of each sub-system can fluctuate. 
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phase transition [46, 47]. This is also seen in maximum entropy calculations, 
discussed in the next section. 



2.2 Maximum Entropy Calculations 

The maximum entropy model calculates the most probable distribution of 
a collection of atoms using the heats of formation A//^. , the ionisation 
potentials, the number of isomers gj and the mean vibrational frequencies of 
all possible species that can be formed from the collection of atoms. It was 
developed to describe the fragmentation patterns of organic molecules [48] 
and has also been applied to atomic cluster fragmentation [49] and to nuclear 
collisions. The most probable distribution is given by finding the maximum 
entropy of the system. The only constraints present are the conservation of 
energy, mass and charge. No assumptions whatsoever are made concerning 
the fragmentation mechanism, only that the energy is statistically distributed 
among all degrees of freedom. The model has been applied to the 
fragmentation of Ceo, Ceo^ [50, 51], La@C 82 ^ [51] and Ci 2 o^ [52]. The simple 
model gives remarkably good agreement with experiment for fullerene 
fragmentation - much better than for organic molecules. This may be partly a 
consequence of the relative simplicity of the fullerenes, consisting of only 
one atomic element, but is also a strong sign of the highly statistical nature 
of the fullerene fragmentation. Starting with a given temperature of the 
fragments, the model iterates through the translational, vibrational and 
rotational partition functions to find the maximum total entropy. It can be 
shown [49] that the abundance of species j is given by: 

= gjQ y exp [-pMl] - Yji. - r,Cj ) (3.2.4) 

where gj is the number of isomers of species j, Qj is the partition function per 
unit volume, V is the interaction volume, rij and cj are the number of atoms 
and the charge of species j, respectively. A//* is the heat of formation 
relative to that of the parent species AH* = AH. - AH^^^^^ '‘.The Lagrange 
multipliers and and the initial mean internal energy of the parent 
molecule (£) have to be determined in order to find the maximum entropy 
starting with a given fragment temperature. The mean energy is given by 

= (3.2.5) 



* Note that this was not done correctly in the original application of this model to Ceo 
fragmentation [50] but was corrected in later work [51, 52]. 
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The mean energy of species j is 



(£^ = ^£yXy 
i 



(3.2.6) 



where e.j is the energy of quantum state i calculated from a common zero of 
all isomers j adding the vibrational, rotational and translational energies to 
the heats of formation. Note that the heats of formation AHj and are 

negative. AH* will therefore be a positive number. 



IJ I I I J 



(3.2.7) 



and the fraction of molecules of species j in the quantum state i is 

^„ = e;‘exp(-/;(£,-AH*)) (3.2.8) 



The partition function can be calculated from standard partition functions for 
the vibrational, rotational and translational degrees of freedom, given e.g. in 
the JANAF tables [53]. The conservation of mass and charge provides the 
constraints 



N = '^n,X. (3.2.9) 

j 

where TV = 60 for Cao and 

^ = (3.2.10) 

j 

with C=0 for a neutral parent molecules and C=1 for a singly charged 
cation. The interaction volume V (eq.3.2.4) is a parameter that influences 
the relationship between the initial internal excitation energy of the parent 
molecule and the resulting fragment temperature but does not critically 
affect the qualitative excitation energy dependence of the fragment 
distributions. For all applications to fullerene fragmentation, V was set to a 
value equal to approximately three times the volume of the parent molecule. 

The heats of formation, ionisation potentials, vibrational frequencies and 
number of isomers are rather well known for the carbon species from C to 
Ceo and, where unknown, it is possible to extrapolate from known values to 
obtain a reasonable "guestimate". For the parameters used in the calculations 
see [50] and references given therein. Extrapolations were also made to 
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obtain the input parameters for calculations involving Ci 2 o^ [52] and 
La@C82M51J. 

The abundance of a particular species predicted by the model depends, to 
a large extent, on the population density in the available translational and 
vibrational states. To create smaller fragments, the difference in binding 
energies between the parent fullerene and the particular fragment has to be 
overcome. At low energies not enough energy is left for higher lying, more 
degenerate states to be populated. The larger fullerenes with a lower barrier 
to overcome are therefore favoured at lower excitation energies. With 
increasing excitation energy, entropy starts to play a role. The smaller 
fragments become dominant as they are favoured due to their high number 
of translational states. The model distributions are therefore obtained as a 
consequence of competition between these two tendencies. Figure 3.10 
shows the results of maximum entropy calculations of the fragmentation of 
Cgo^- The upper part shows the calculated fragment distribution for an initial 
parent excitation energy <E> = 110 eV. The typical fullerene bimodal 
fragmentation pattern is clearly reproduced. The relative intensities of the 
various fragments are critically dependent on the heats of formation and 
ionisation potentials used in the calculation [51]. No attempt has been made 
so far to optimise this to obtain perfect agreement with experiment (a change 
of 0. 1 eV in binding energy can make a significant difference to the relative 
intensity of a given species). Instead the literature values and simple 
extrapolations from these have been used. The main features observed in low 
energy collision experiments and ns laser excitation work are reproduced. 
The fragmentation pattern moves from that of sequential C 2 evaporation 
from fullerene ions at low excitation energies to a bimodal distribution 
containing fullerenes ions as well as small fragment ions with predominantly 
ring and chain geometries as the excitation energy increases beyond ca. 70 
eV. The magic numbers observed in the experiments for e.g. n = 7, 11, 15 
etc. are reproduced. For very high excitation energy only small fragment 
ions, predominantly and € 3 "^ remain. The lower part of the figure plots the 
relationship between the mean initial excitation energy <E> and the resultant 
fragment temperature. Flere one can clearly see the signature of a phase 
transition that occurs at a temperature of slightly less than 4000 K. Evidence 
for a phase transition at such elevated temperatures is also found in 
molecular dynamics simulations [46] (section 3.3). The excitation energy at 
which the phase transition onsets corresponds to the threshold energy for 
production of the small ring and chain fragment ions Cn"^ with n < 30. 
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Number of Carbon Atoms in Fragment Ion 




Figure 3.10. Maximum entropy calculations for Ceo^. (a) Calculated fragment distribution 
for an initial internal excitation energy of 1 10 eV (b) relationship between average internal 
energy <E> and fragment temperature showing the signature of a phase transition. 



3. MOLECULAR DYNAMICS SIMULATIONS 

Molecular dynamics (MD) simulations are a very powerful tool to help 
visualise and intuitively understand the dynamics of collisions and half- 
collisions. Many studies involving fullerenes have been carried out. It is, 
however, important to remember that the results of such calculations can be 
critically dependent on the parameters used in the particular study and one 
should always be careful not to over-interpret the data however tempting it 
may sometimes be. It is also the case that MD simulations are generally only 
valid as long as there is no influence of electronically excited states. It is 
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only very recently that methods are being developed that are able to treat 
electronic excitation in dynamical studies of fullerenes [54, 55]. Generally, 
MD simulations involve the calculation of trajectories of an ensemble of 
particles in space-time. The particles move due to the forces acting on them 
according to Newton's second law. The equation of motion of a particle of 
mass m moving in the presence of (N-1 ) interacting particles is thus: 

F, = (3.3.1) 

at 

j*i 

where r,y is the distance of the particle i to another particle j and Ffry) is the 
force on the particle i exerted by particle j. If one knows the interaction 
potential at any given time for any given set of coordinates then one can 
calculate the trajectory of the particle using eq.(3.3.1) and the relation 
F(r) = -VV(r). To realise this on a computer, one makes use of time- 
discretisation techniques using finite difference methods [56]. One of the 
most common methods is to make use of the Verlet algorithm [57]: 

r, (t + St)^ 2r. (f ) -r, (f - St) + —F. (t) (3.3.2) 

m 

where is a suitable time step, normally on the order of 10'*^ to 10’*® s. The 
velocity of the atoms can be obtained from 

\i{t)-\i(t-St) + —F.(t) (3.3.3) 

m 

The major problem of molecular dynamics simulations is to define the 
correct interaction potential. It is often a very non-trivial task and most of the 
computation time is allocated to the calculation of the potential. There are 
different ways of doing this. In the following sections the main methods 
applied to fullerene collisions will be outlined. General introductions to MD 
simulation techniques can be found in e.g. [58]. 

3.1 Classical MD Simulations 

Due to the large number of atoms involved in fullerene collisions and the 
need to calculate many trajectories in order to obtain a statistically relevant 
picture of the collision dynamics, classical MD simulations have been 
extensively used. Classical MD simulations are fast, computationally 
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efficient and based on the construction of parameterised, empirical 
potentials. Carbon has the ability to form many kinds of structures and this 
must be accounted for in the construction of the empirical potentials, 
especially if fragmentation processes are to be investigated. The treatment of 
the covalent character of the carbon-carbon bonds requires the incorporation 
of three-body forces. One of the most common and successful 
parameterisations for carbon is due to Tersoff [59]. This works very well for 
bulk diamond-like structures but is less suitable for treating graphite and is 
not able to correctly reproduce the structure of the small carbon molecules. 
This deficiency was overcome by Brenner who generalised the Tersoff 
potential to be able to treat hydrocarbons and small carbon molecules at the 
expense of adding more parameters [60] (with corrections to the original 
paper given in [61]). 

The Brenner parameterisation has been found to be very suitable for 
modelling fullerenes [62]. The carbon-carbon interaction potential has the 
form 



‘''" = jii[»'.(':,)-s/4':>)] (3.3.4) 

where V^(ry) and VA(ry) are the repulsive and attractive parts of the potential, 
respectively. The parameters for the carbon-.carbon interaction can be 
written as 



('[,■) = /c {>'y)UR,A exp {-K,Ahj 

k=\ 

k^ij 



G{0)^b, 



l + ."o 



^0 (^ 0 ^ -t (l-tcos(6’)) 



(3.3.5) 



where 
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(3.3.6) 



Bij is the bond-order function of three-body bond-angle dependent 
interaction terms, 9 is the angle between two bonds, G( 9ijk) is a function of 
the angle between bonds i-j and i-k, fdrn) restricts the pair potential to 
nearest neighbours, Rc is the equilibrium distance, Dc is the well depth and 
are Morse parameters. The function G( 0) is parameterised with respect 
to the binding energies of the C 2 molecule [60]. N is the number of carbon 
atoms in the system. The values of the various parameters found to be most 
appropriate for fullerene simulations are given in Table 3-2. 



Table 3-2. Brenner parameters used for classical MD simulations involving fullerenes. Taken 
from [60]. All values in atomic units. 



Ur 


19.0496 


Ua 


12.0544 


Ar 


1.27498 


Aa 


0.98835 


bo 


1.1304x10'^ 


Co 


19 


do 


2.5 


Rc 


3.496 


Dc 


0.283459 



The above parameterisation and other similar approaches (e.g. [46]) have 
been used to simulate "half-collisions" in which the dynamical behaviour of 
a highly excited fullerene is studied as a function of time. They have also 
been used to study the behaviour of fullerene-fullerene collisions [63] (see 
Chapter 9). 

In many cases one wishes to simulate the dynamical behaviour in 
collisions between a fullerene molecule and an atom. This is important for 
understanding the collision energy dependence of fragmentation or for 
understanding the mechanisms by which atoms are captured in gas phase 
collisions to produce endohedral fullerenes (see Chapter 8). In this case the 
interaction between the foreign atom and the carbon atoms is most simply 
described by a screened Coulomb potential [62, 64] 

V 7 7 ^2 

M Gj 



(3.3.7) 
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where Z; and Z 2 are the charge numbers of the atoms, N is the number of 
carbon atoms in the fullerene, 2" is the Moliere screening function [65] and a 
is an appropriate screening radius. An example of an MD simulation 
showing capture of a Ne atom by Ceo at a collision energy of 80 eV is shown 
in Fig. 3.11. These calculations can provide very useful qualitative 
information concerning the dynamics of fullerene collisions but are likely to 
be less reliable when describing the details of the fragmentation behaviour 
[66]. It should also be remembered that such MD simulations can only 
reliably and realistically calculate the time dependence on a scale of ps. 
Experiments, on the other hand, normally detect the products of collisions on 
a time scale of at least many microseconds after the collision event. Some 
attempts have been made to combine MD studies with statistical models to 
reproduce experimental fragmentation cross sections [67]. 



3.2 MD Simulations using Density Functional Theory 

A more reliable but considerably more time consuming treatment uses ab 
initio MD concepts based on density functional theory (DFT). DFT was 
created in 1964 with the work of Hohenberg and Kohn [68] who proved that 
the electronic ground state energy and all related properties can be defined 
from the knowledge of the electronic ground state density alone. Slightly 
later, Kohn and Sham showed [69] that the exact electronic ground state 
energy Eg for a system of N electrons and M nuclei with the ground state 
electron density p can be written as 



^ i=l A=1 V 



+ - 
2 



1 jj/' (>')/' (O 



r-r 



)^drdr'+E^^^[p{r)] 



(3.3.8) 



where the first term is the kinetic energy of the electrons, the second term is 
the electron-ion interaction energy and the third term is the electron-electron 
interaction energy. E^xc is the exchange correlation energy that depends on 
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Figure 3.11. Classical MD simulations of the collision of a Ne atom with Ceo at a collision 
energy of 80 eV. The Ne atom penetrates the cage and is captured to form a highly excited 
endohedral molecule. From [70]. 



the density p. The functions xj/lr) are the so-called "Kohn-Sham" orbitals. 
They have no physical meaning but are used to find the ground state density 
p which is given by 

= ( 3 - 3 - 9 ) 

t=l 

The Kohn-Sham orbitals are obtained by solving the one-electron equations 

(3.3.10) 

where H^s is the one-electron Hamiltonian 



KS 2 






+n,,(r) (3.3.11) 



and Vexc is the exchange-correlation potential. According to [69], this 
potential is a functional derivative of E^xc with respect to p 
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y 

Sp ■ (3.3.12) 

In the ideal case the exchange-correlation potential handles all the 
electronic contributions not accounted for by the third term in eq.(3.3.8) and 
is only a functional of the electron density, p. In general, DF methods differ 
from each other in the form of the exchange correlation energy Eexdp]- The 
solution for Eg is found interatively: 

• One starts with an initial guess for p and then calculates the initial 
Yexc from eq.(3.3. 12) 

• The Kohn-Sham orbitals, y/j, are numerically calculated using 
eqs.(3.3.10 and 3.3.11). 

• The thus determined y/-?, are used in eq.(3.3.9) to obtain a new value 
of p, 

• One starts again from the beginning of the cycle with the new p until 
the calculation converges. The ground state energy, E„, is then 
obtained from the converged p using eq.(3.3.8). 

The main error in the calculation results from the approximation used for 

Eexc- The local density approximation (LDA) is widely used, whereby 

Eexc\^P{^)\ depends only on the electron density 



= (3.3.13) 

Calculations in which the electronic states were treated within the LDA 
with a Fermi-Dirac distribution of the excited levels were carried out for 
LC + C^ collisions [71]. These calculations showed the possibility of the 
direct insertion of the Li ion to form an endohedral fullerene (see Chap. 8). 
However, there have been very few truly ab initio molecular dynamics 
simulation carried out for fullerene collisions. The high demands on 
computing capacity that are needed for simulating such large systems have 
stimulated the development of semi-empirical techniques which can have 
almost the accuracy of self-consistent LDA techniques but are much less 
demanding on computer time. 

Tight-binding schemes have generally been very successful and can be 
seen as simplified ab initio methods which still use quantum mechanics to 
determine the electronic properties of the system. In these calculations the 
elements of the Hamiltonian matrix are fitted to reproduce an arbitrary set of 
input data. Normally, only two-centre contributions to the Hamiltonian 
matrix are considered. However, the fitting can be a rather complicated and 
difficult procedure. Zhang et al. have applied such calculations to simulate 
fullerene collisions [72]. 
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Another semi-empirical approach that avoids the problem of fitting the 
matrix elements has been applied extensively to the study of fullerene 
collisions [73]. This is based on the work of Seifert, Eschrig and co-workers 
[74, 75] and applies the formalism of optimised linear combination of atomic 
orbitals (LCAO). In this approximation, the Kohn-Sham orbitals are 
expanded in terms of atom-centred localised basis functions (r - ) 

(centred at nucleus J). 



(3.3.14) 

ft 



The Kohn-Sham equations (3.3.10) can be transformed into a set of 
algebraic secular equations 



where the matrix elements are defined as 



(3.3.15) 



(3.3.16) 



and the Kohn-Sham one-electron Hamiltonian is given by eq.(3.3. 1 1). 

The interatomic forces for MD applications can be derived from a 
calculation of the gradients of the total energy at the atom sites, divided into 
an electronic and a nuclear part 

F =-V£ = F, -yv^^ (3.3.17) 

J J ,e U U 

J^I 



where the last term is the repulsion due to the other nuclei. Fy,,, can be 
written as a sum of orbital contributions 

F.,. = (-VV + V + (3.3.18) 

i fi y 

where n, is the occupation number of the ith orbital, are the orbital energy 
eigenvalues and V ^ are the matrix elements of the last two terms in 
eq.(3.3.11). 
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It has been shown [74, 75] that the nuclear repulsion force from 
eq.(3.3.17) and the contribution arising from the last term in eq.(3.3.18) 
largely compensate each other and these two terms can be described together 
as an empirical repulsive force parameterised by a potential of the form 

(3.3.19) 

[0, R> R, 

where the parameter a is the vibrational frequency and Ro the equilibrium 
bond distance for the C 2 molecule. 

This model is thus a hybrid of the DF methodology together with 
empirical potentials. The Newton equations of motion are solved 
simultaneously with the approximate Kohn-Sham equations that are solved 
self-consistently. 

3.3 Non-adiabatic MD Simulations 

The above MD treatments can only describe systems in the electronic 
ground state. This is not such a problem for low energy collisions where the 
excitation is predominantly via vibrational excitation but it is a problem for 
high energy collisions and electron impact collisions (as well as laser 
excitation) where one can expect a substantial amount of electronic 
excitation. Time-dependent density functional methods have been developed 
in recent years that can describe photoabsorption spectra and excitation 
energies of atoms and small molecules [76], however, it is still 
computationally very challenging to apply such methods to collisions 
involving molecules as large as fullerenes. Kunert and Schmidt have applied 
a model to ion fullerene collisions where the electronic and vibrational 
degrees of freedom are treated simultaneously and self-consistently by 
combining time-dependent density functional theory with classical MD [54, 
77]. They studied the energy loss for various ion - Ceo collisions (H"^, C^, 
Ar"^) as a function of impact velocity in the range 0.01-0.5 a.u. The relative 
contributions (electronic and vibrational) to the total excitation energy were 
shown to depend strongly on ion mass and impact velocity up to a value of 
0.25 a.u. At higher velocities the excitation energies were predicted to be 
velocity-independent, in contradiction with the "stopping power" predictions 
discussed in section 3.1.1. Results of Ar^Ceo collisions are shown in Fig. 
3.12 for two different impact velocities (0.02 and 0.45 a.u.). At the lower 
velocity, the excitation is predominantly vibrational and leads to a very rapid 
(< 100 fs) break-up of the carbon cage. At the higher velocity, no momentum 
is transferred to the cage on impact but the fullerene becomes electronically 
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excited (with approximately the same total excitation energy as for the low 
velocity impact). The cage starts to visibly vibrate after about 45 fs but 
fragments on a much longer time scale than in the first example. This type of 
non-adiabatic MD simulation is potentially a very powerful tool to predict 
and demonstrate the dynamics of collisions involving large molecules and 
clusters and it is to be hoped and expected that further development and 
application of such models will occur in the near future. 



t=0 






77 fe 



45 fs 



.-N 










Figure 3.12: Snapshots of the reaction dynamics as a function of time for Ar'^-Qo collisions 
calculated using a combined time-dependent DFT, classical MD approach, (a) b = 2 a.u., v = 
0.02 a.u. Excitation is predominantly vibrational (b) b = 2 a.u., v = 0.45 a.u. Excitation is 
predominantly electronic. Although the total excitation energy is approximately the same, the 
nature of the excitation strongly influences the time scales and mechanisms of fragmentation. 

Adapted from [54]. 
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4. CHARGE TRANSFER 

Many experimental studies have been made of charge transfer collisions 
involving fullerenes (see Chaps. 4,6 and 9). In order to interpret the results 
of experiments, relatively simple classical or semi-classical models are 
invoked that were originally developed to describe atom-ion collisions. 

4.1 Classical "Over-the-Barrier" Model 

The most extensive charge transfer studies in the literature involve slow 
collisions between highly charged ions and neutral fullerene targets (Chap. 
6). In this case the classical "over-the-barrier" model is often used as a basis 
for models used to describe and interpret the experimental results [12, 78]. 
This approach considers the classical situation of a slow ion with charge q 
(treated as a point charge) at a distance R outside an infinitely conducting 
sphere with radius r$. Fig. 3.13. The potential outside the sphere consists of 
two parts: one Coulombic term from the ion itself and one part from the 
image charge induced by the ion. If an electron tries to leave the sphere it 
will feel forces from the ion and its induced dipole, from the dipole that it 
induces itself and from the electrically charged sphere. The sphere surface 
must remain an equipotential surface for all distances R between the ion and 
the centre of the sphere. This criterion gives a unique solution for the 
magnitude q ' = -qr^ / R and position R' = r^ I R of the image charge. In order 
to conserve the total charge, an equally large positive charge -q' = +qr^ / R 
must be placed at the centre of the sphere. The energy required to ionise the 
sphere in the presence of the projectile ion is shifted due to the shift in the 
electric potential, (j). For the point on the surface of the sphere closest to the 
projectile ion this is given as 



= — i i 

R-rs r^-R' h 

_ q q ^ ‘It's _q 

R — r^ R — r^ r^R R 



(3.4.1) 



so that the required ionisation energy, IP* = IP-i-q/R. 

The critical distance between target and projectile at which the first 
electron may leave the sphere is given by the classical over-the-barrier 
criterion V[,' (max) = //]* where V[,‘(max) is the maximum of the potential 
seen by the electron when it moves towards the projectile ion. It can be 
written in terms of the electron's interactions with the projectile charge, the 
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q 



Figure 3.13. An ion of charge at a distance R outside an infinitely conducting neutral sphere 
of radius rg. The position of the image charge q' = -r^q/R is given hy R'=rg/R. An equally 
large positive charge is positioned at the centre of the sphere. Adapted from [12] 



induced charges and the charge left behind on the sphere (represented by a 
+1 charge at the sphere centre) 
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(3.4.2) 



The last two terms in eq. (3.4.2) are due to the electron's interaction with its 
self-induced image charges. One can easily extend this to the sequential 
transfer of electrons. In this case the potential energy for transfer of the ith 
electron from an ( i-1 )+ charged target to a (q^-i-i-i)-i- charged projectile is 



{q-i + l) ^ r^{q-i + l) 
R — x Rx — r^ 



.[q-i + l) + lR 1 



Rx 



+ - 
2 




(3.4.3) 



For electron transfer to occur this should be equal to the shifted ionisation 
energy for the removal of an electron from the (i-l)+ target ion, 
IP* = IP. + {q-i + l)/R . 

The potentials that fulfill the over-the-barrier criterion for transfer of 
electrons between a Ceo and an Ar*"^ projectile ion are illustrated in Fig. 3.14. 
The full line illustrates the situation for transfer of one electron from the 
neutral fullerene to the projectile ion. The critical distance for this is 
approximately 26 ao- The dashed and dotted lines show the situation for 
transfer of one electron from the singly charged target to the 1+ projectile 
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and from the 7+ charged target to the 1+ charged projectile ion, respectively. 
As the projectile ion moves past the fullerene target it may be fully 
neutralised in a single collision following the sequential transfer of electrons. 
Note that even for electron transfer from to Ar"^ the critical radius at 
which transfer can occur is still significantly larger than the fullerene cage 
radius. For these calculations the radius of the sphere is taken to be 8.2 ao, 
derived from the experimental polarisability of free Ceo molecules [20]. The 
total electron transfer cross section in collisions with a projectile of charge q 
can then be estimated from 



K, = ^RMf 



(3.4.4) 




Figure 3.14. Potential energy seen by the active electron in Ar®"^ + Qo collisions at the critical 
distances for charge transfer where the height of the potential barrier matches the shifted 
ionisation potential I*. Full line: transfer of the first electron from the neutral sphere to Ar*"^. 
Dashed line: transfer of the second electron from C^o* to Ar’"^. Dotted line: transfer of the 
eighth electron from to Ar"^. 

The above version of the over-the-barrier model assumes that the 
fullerene behaves like an infinitely conducting sphere. In other words, the 
charge mobility is assumed to be so high that all effects of localisation of 
individual charge carriers are averaged out. Another approach has also been 
applied to fullerene collisions [79-81]. In this second approach the positive 
holes left behind on the fullerene cage are assumed to be movable on the 
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molecular surface in a localised fashion. During the time it takes for the 
electron to move from the sphere to the ion, the hole is assumed to be 
localised closest to the projectile. When the electron has reached the ion, the 
hole moves towards the position that minimises its potential energy. In this 
case the energy needed to ionise a neutral sphere is IP* = IP^ + ql{R-a) . It is 
more complicated to calculate the electronic potential with this model, 
especially if more than 2 electrons are involved. The results are not too 
different from the infinitely conducting sphere model. The discrepancies are 
larger for the transfer of a small number of electrons but as the number of 
holes on the sphere surface increases, the electric field outside the sphere 
approaches that of the infinitely conducting sphere model. For typical 
projectile velocities, the times between the successive critical distances are 
in the range 0.1-1 fs which implies that the hole re-arrangement times must 
be on the order of 0. 1 fs or faster. 

4.2 Rapp and Francis Model 

The Rapp and Francis model [82] provides a semi-quantitative 
calculation of the collision energy dependence of charge transfer in atom-ion 
collisions 



A^ + B^A + B^+^E (3.4.5) 

where AE is the difference in ionisation potentials of A and B. This has been 
applied to charge transfer collisions between neutral fullerenes and fullerene 
ions, discussed in Chap. 9. The method is applicable in the collision velocity 
regime where the semi-classical impact parameter method is applicable 
(roughly (10^/p'^^)ms ' < v < 10® ms ^ where p is the reduced mass in atomic 
units). If one first considers symmetric charge transfer, -i- A — > A -i- AA, 
i.e. AE = 0, the probability of charge transfer in a collision with velocity v 
and impact parameter b is 



P{b,v)- sin^ 






2hv 



dx 



(3.4.6) 



This expression is obtained from considering the A 2 "^ collision complex 
formed during the collision as a one-electron problem. The nonstationary 
state (representing the collision) can be expressed in terms of symmetric and 
antisymmetric stationary states with energies £« and E^, respectively. The 
cross section for charge transfer is then calculated from 
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c7(v) = 2;r£ P{b,\)bdb 



(3.4.7) 



For ion-atom collisions this can be calculated by using the wavefunctions 
of the molecular ion. For more complex reactants, including fullerenes, this 
is not feasible and Rapp and Francis adopted an approach using semi- 
empirical orbitals. For R/uo >1, with R the internuclear separation, this 
becomes, in units of eV [83] 



(£,-£j = 4(/F)(R/a,)exp(-^(/F)/13.6(R/ao)) (3-4.8) 



where IP is the ionisation potential of atom A. Combining this with eq. 
(3.4.6) one can obtain an approximation to the charge transfer probability 
[82] 



Po(6,v) sin^ 



I 2ti /P ^3/2 



ya^ hv 



^ a ^ 
1 + ^ 

V 



exp 



yb'' 
V 



(3.4.9) 



with y = V/P/13.6 . Over a fairly large velocity range this leads to a charge 
transfer cross section with the simple form 



(7 =-^jlnv-l-^2 



(3.4.10) 



where 






(3.4.11) 



and ^2 is a rather complicated expression, depending only on IP and a 
velocity-averaged impact parameter for which P(6,v)=0.25 [82]. 

For asymmetric charge transfer (eq. (3.4.5)) we consider that angular 
momentum and spin are conserved during the charge transfer reaction. This 
is justifiable since charge transfer occurs at fairly large impact parameters 
for collisions in the "intermediate" velocity range where the model is 
applicable. In this case we only need to consider charge transfer between the 
two electronic states of AB"^ with identical symmetry which go to (A"^ -i- B) 
and (A -i- B"^) as asymptotic forms. This is known as a "two-state 
approximation", it implies that charge transfer involving other asymptotic 
states with a high AE is disregarded. The fraction of states that allow charge 
transfer (i.e. have identical symmetry to the initial collision partners) is 
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denoted as /, the statistical weight factor. The calculation of the charge 
transfer cross section is quantum mechanical. The LCAO ansatz is used to 
obtain the total wave function for the active electron in the combined field of 
A-" and 

V^c^{t)(/)^ {r^)exp{-i(Oj) + Cg {t)(p^ {r^)exp{-icOgt) (3.4.12) 



where the c's are time-dependent coefficients of the expansion, (pA and (j)B the 
atomic orbitals for the electron on on nuclei A’^ and B"^, and are the 
distances of the electron from A"^ and B"^, co^=e/h and Ea and Eb are the 
ionisation potentials of A and B in states and respectively. These 
states are solutions of the atomic Schrddinger equations 



2m^ 












(3.4.13) 



in which j can be either A or B and V) are the effective potentials binding the 
electron to the nuclei. The time-dependent coefficients Cj are then obtained 
by solving the time-dependent Schrddinger equation 



2m^ 




+ V,{r,) + V,{r,) 



Iff -hi 



diff 



(3.4.14) 



Eqs. 3.4.13 and 3.4.14 thus provide a system of equations for the coefficients 
Cj and their time derivatives. The reaction probability, P(b,v) is then given by 

T’(Zi,v) = |c^ (t — > °o)|/ . (3.4.15) 



After a number of rearrangements and approximations that will not be 
detailed here [82], the charge transfer cross section can be obtained from 



f sech^ 


CO a^Tcb 


J 

0 


< 



27Tbdb 



(3.4.16) 



The impact parameter b / is obtained from the relationship 



sech^ 



0) a^TTb 

v]l 2r 



4fo(bi,v) 



(3.4.17) 
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where Po is obtained from eq. 3.4.9 and a> = AE/h. 

Examples of the application of this model to charge transfer collisions 
involving fullerenes are discussed in Chap. 9. 
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Chapter 4 

CHEMICALLY REACTIVE AND THERMAL 
ENERGY COLLISIONS 



In this chapter we will review collision experiments at thermal or quasi- 
thermal energies that have been used to study mainly charge transfer and 
collision reactions. The first part of the chapter will consider the wide range 
of reactions that have been studied using the SIFT method (selected ion flow 
tube). We will also touch on some complementary techniques applied to 
chemical reaction studies involving fullerenes, also at higher collision 
energies. In the second part of the chapter we will discuss Penning ionisation 
studies of collisions between neutral fullerenes and electronically excited 
rare gas atoms. 



1. LOW ENERGY ADDITION AND CHARGE 

TRANSEER REACTIONS WITH EULLERENES 



1.1 Addition Reactions 

Low energy addition reactions and charge transfer involving fullerene 
ions have been extensively studied by Bohme and co-workers [1-4] using the 
SIFT technique (see Fig. 2.8). The reactions of singly, doubly and triply 
charged Ceo with a wide range of collision partners were studied at room 
temperature (294 ± 2 K) and at a helium bath pressure of 0.35 ± 0.01 Torr 
[ 2 ]. 
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The positive charge on the fullerene can be regarded as being delocalised 
over the carbon surface. The interaction with an incoming molecule will then 
localise the charge as a consequence of the electrostatic interaction between 
the charge and the induced or permanent dipole of the molecule. Chemical 
addition can then occur at the site of the charge. However, the carbon atom 
on the cage that undergoes the chemical bond has to experience a change in 
hybridisation from sp^ to sp^ that implies that exothermic covalent bond 
formation may involve an activation barrier. The singly-charged fullerene 
cation is in fact rather unreactive towards many molecules [2]. The reactions 
that have been observed in SIFT experiments are summarised in Fig. 4.1. As 
the charge on the cation increases, it becomes easier to derivatise the cage 
and considerably more reaction products could be observed in collisions 
with than with [2]. This is presumably because the stronger 

electrostatic interaction between and an incoming molecule can help to 
overcome the activation barrier associated with the hybridisation of the 
bonding site. As the positive charge on the fullerene increases further, the 
ionisation energy of the fullerene increases and electron transfer becomes 
more exothermic and kinetically more favourable, in spite of the activation 
barrier that arises from the production of two charged products. Electron 
transfer therefore becomes more likely as a competing or dominant reaction 
channel. For this reason the number of observed addition reactions with a 
triply charged fullerene cation is lower than that for a doubly charged cation 
[ 2 ]. 

Interesting trends are seen in the reactions of fullerene cations with 
ammonia and amines [5]. Ammonia adds only very slowly to the singly 
charged cation but the rate of addition approaches the collision rate with 
increasing alkyl substitution. Table 4.1. The reactions occur exclusively by 
adduct formation 



C,; + M^C,oM^ (4.1.1) 

where the dot indicates a free radical species. The adduct formation is 
attributed to the formation of a covalent C-N bond with the fullerene cation 
through donation of an electron pair from the nitrogen atom to the site of 
charge concentration on the surface of the fullerene, as illustrated in Fig. 4.2. 
This requires a transition from a reactant ion with extreme charge 
delocalisation to a product ion where the charge is extremely localised on the 
nitrogen atom. This formation mechanism is quite different to the proposed 
mechanism for amine addition to neutral Ceo in solution which involves two 
steps: single electron transfer to a 7i bond followed by covalent bond 
formation and proton transfer [6]. 
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Figure 4.1. Overview of addition reactions of observed with a SIFT apparatus at room 
temperature and at a helium bath gas pressure of 0.35 Torr. From [2]. 

Table 4-1. Products and rate coefficients for reactions of Cso* with ammonia and amines at 
294 K in helium, determined using the SIFT technique. From [5]. The rate coefficients are 
given in units of 10 "^cm^molecule‘*s''. kc is the collision rate coefficient. 



neutral reactant 


product 


^bs 


kc 


NHj 


CsoNH3^ 


< 0.001 


1.69 


CH 3 NH 2 


C6oCH3NH2-^ 


0.015 


1.44 


CH 3 CH 2 NH 2 


C6oCH3CH2NH2-^ 


0.050 


1.34 


(CH 3 ) 2 NH 


C6o(CH3)2NH+ 


0.85 


1.22 


(CH3)3N 


C6o(CH3)3N" 


1.0 


1.05 



The observed increasing reaction rate with increasing degree of aliphatic 
substitution of the amine can be explained by considering the association to 
occur in two steps. Firstly, a collision complex is formed which is stabilised 
by colllision with the He atoms of the buffer gas. The collisional stabilisation 



(CgoNRj^ )* +He ^ CgoNRj^ + He’ 



( 4 . 1 . 2 ) 
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( 2 ) 



Figure 4.2. Reaction of ammonia and amines with the Cso* cage. From [5]. 
competes with the unimolecular dissociation of the collision complex 
(C,oNR3^)*^C,; + NR3 (4.1.3) 



Application of the steady-state assumption to the concentration of the 
collision complex leads to a termolecular rate constant given by 

where kc and kg are the rate constants for the formation of the collision 
complex and its collisional stabilisation and Td is the lifetime for 
unimolecular dissociation [5]. The large increase in the rate constant is due 

k(3) = kcksT^ (4-1.4) 



primarily to To since kc and kg are not expected to be very sensitive to 
aliphatic substitution. Various statistical theories can be used to provide an 
estimate of To of which the simplest is [7] 



F rt Ff. 



(D + 3RT) 
3RT 






(4.1.5) 



where Fq is a collision lifetime corresponding to the frequency of vibration 
along the reaction coordinate, D is the dissociation energy of the collision 
complex, S is the number of effective degrees of freedom available in the 
complex for energy transfer and R is the gas constant. Although Ceo has 60 
atoms, it has no internal modes of rotation and, given its highly rigid 
structure, the vibrational modes are widely space and will not allow for a 
great deal of energy dispersion within the collision complex under the 
experimental conditions. Trimethylamine contains much fewer degrees of 
freedom than Ceo however, the methyl groups have rotational modes that 
have a much higher energy density than the vibrational modes of Ceo- The 
bending vibrational modes of trimethylamine will also have a relatively high 
energy density, suggesting that they may also aid in energy dispersal within 
the collision complex and thus prolongation of the complex lifetime. 
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The doubly charged fullerene cations show quite different behaviour, 
Table 4.2. In this case adduct formation competes with charge transfer. 
Charge transfer is endothermic for the reactions with the singly charged 
cation but exothermic for reactions with the doubly charged cation. 
However, Coulomb repulsion between the singly charged cations produced 
in a charge transfer reaction introduces a barrier to charge transfer. Therefore 
charge transfer should only be observable at room temperature for reactions 
that are sufficiently exothermic to overcome this barrier [8]. The ionisation 
potential of Ceo^ is 1 1 .4 eV which implies that charge transfer with all of the 
amines in Table 4.2 is exothermic by more than 2 eV. Charge transfer is 

Table 4-2. Products and rate coefficients for reactions of with ammonia and amines 
at 294 K, determined using the SIFT technique. Rate coefficients in units of 10'^ 
cm^molecule'*s‘*. From [5]. 



neutral reactant 


product 


branching ratio 


kobs 


K 


NHj 


C6oNH3^-^ 


1.00 


1.2 


3.37 


CH3NH2 


C6oCH3NH2"-^ 


0.90 


2.6 


2.88 




C6o" + CH3NH2-^ 


0.10 






CH3CH2NH2 


C6oCH3CH2NH2^-^ 


0.70 


2.6 


2.68 




C6o" + CH3CH2NH2-^ 


0.30 






(CH3)2NH 


C6o(CH3)2NH"+ 


0.10 


2.9 


2.44 




Ceo^ + (CH3)2NH+ 


0.90 






(CH3)3N 


Cfio" + (CH3)3N" 


1.00 


2.2 


2.10 



indeed observed to compete with association for all the amines and, 
according to the measured branching ratios, does so increasingly as the 
charge transfer reaction becomes more endothermic until it is the only 
observed channel in the reaction with trimethylamine which has the lowest 
ionisation potential (7.82 eV) of all the amines in the study. No charge 
transfer is observed with ammonia (IP = 10.16 eV) even although this 
channel is also exothermic. However, the exothermicity of the charge 
transfer reaction (1.23 eV) is less than the contribution to the barrier height 
which arises from the Coulomb repulsion [5]. The situation is illustrated 
schematically in Fig. 4.3. 

No association reaction products were found in reactions between the amines 
and Ceo^^- The triply charged cation reacted exclusively by charge transfer, 
in accordance with the much higher ionisation energy of Ceo^^- 

Triple quadrupole mass spectrometry experiments have also been applied 
to the study of the association reactions of fullerene ions [9]. Several factors 
have been shown to influence the addition reactions: the charge state and 
size of the fullerene ion and the size, structure and degree of unsaturation of 
the target. 




76 



Chapter 4 



The reactivity of fullerene anions has been studied with the flowing 
afterglow technique [10]. They were shown to be rather unreactive with a 
variety of acids but did react with NO 2 . 



> 

o> 

l_ 

(13 

C 

LU 

(D 

> 

0 
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Figure 4.3. Hypothetical potential energy / reaction coordinate profile which accounts for the 
occurrence and nonoccurrence of charge transfer in reactions between C(,a* and amines. The 
interaction between the product ions is dominated by Coulombic repulsion. 

Time of flight mass spectrometry experiments, at much higher collision 
energies than discussed above, have also been applied to the study of 
fullerene reactions. Callahan et al. carried out detailed studies of the reaction 
products in collisions between fullerene ions and O 2 [11]. Their results shed 
some light on the photostability of fullerenes in air and indicated that the 
reaction with O 2 and O 3 (but not O) led to the disruption of the cage structure 
and fragmentation to produce CO, CO 2 and C 2 O 3 . The paper also discussed 
an empirical link between the extractability of macroscopic amounts of 
endohedral metallofullerenes and their oxidative stability. 

Guided ion beam experiments were applied to the study of reactive 
collisions between neutral fullerenes and atomic cations. Some of these 
results involve the endohedral capture of the ion inside the carbon cage [12- 
15] and are discussed at length in Chapter 8. Three different types of 
bonding were studied in the experiments: the formation of weakly bound 
substituents attached to the carbon cage, the replacement of carbon atoms in 
the fullerene cage by heteroatoms and the above-mentioned formation of 
endohedral complexes. A summary of the results obtained for the different 
collision partners investigated is given in Table 4.3. One of the most 
interesting systems was - 1 - Ceo where evidence was found for both a 
weakly bound C attached to the cage and for a substitution reaction in hich 
the projectile ion replaced one of the carbon atoms [16]. The cross sections 
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measured for the range of products is shown as a function of laboratory 
collision energy in Fig. 4.4. The product ion appears over a narrow 
energy window, peaking at an energy of approximately 4 eV. The most 
dominant process over the entire collision energy range investigated is 
charge transfer to produce Ceo^- The nature of the adduct was addressed 

Table 4-3. Summary of results from reactive ion-fullerene collisions obtained by the 
Anderson group [17,18]. 



Projectile ion low energy 
adduct 
formed? 


Approx, 
disappearance 
energy of 
adduct [eV] 


Substitutional Endohedral 
; compound complex 

formed? formed? 


Efficient 

charge 

transfer? 


Li-", Na-", K-" no 


- 


no 


yes 


no 




yes 


10 


no 


no 


yes 




yes 


15 


yes 


no 


yes 


N-" 


no 


- 


yes 


no 


yes 


0* 


yes 


28 


yes? 


yes? 


yes 


F+ 


no 


- 


no 


no 


yes 


Ne-" 


no 


- 


no 


yes 


yes 


Mn-" 


yes 


10 


no 


yes 


yes 


Fe-" 


yes 


15,40 


no 


7 


yes 


C 2 " 


yes 


36 


yes? 


no 


yes 


C0-" 


no 


- 


no 


no 


yes 


Mni"" 


only MnCeo' 


' 13 


no 


only Mn@C6o^ yes 



by using a projectile ion and measuring the isotopomer distribution 
resulting from the process [Cei^]* — > + C. There are two possible 

limiting cases. If the extra carbon atom is only weakly bound to the cage 
then the dominant decay route will be the loss of the added C atom with no 
isotope scrambling. At the other limit, the could incorporate the extra 
atom, producing a distorted fullerene. In this case the atom is more likely 
to be retained in the subsequent decay process. Since the is internally 
hot, bond rearrangements could lead to a completely random distribution in 
the cage and therefore a small probability for losing the '^C. Figure 4.5 
shows the isotope distributions for Ceo^ and produced in the reactions 
between and Ceo [16]. Fig. 4.5(a) shows the peak under conditions 
where the contribution to the signal contains both charge transfer products as 
well as the decay product from Cei"^. The pure charge transfer product will 
simply show the natural isotope distribution. The distribution in Fig 4.5(a) 
shows some evidence for retention of '^C. Parts (c) and (d) of the figure 
show the expected isotope distributions assuming the natural abundance 
and also assuming that there is complete isotope scrambling. In Fig. 4.5(b) 
only forward scattered ions are detected and therefore the pure charge 
transfer signal is discriminated against. Under these conditions the isotope 
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distribution clearly shifts towards a higher content indicating the 
presence of isotopic scrambling in the cage [16]. Tight binding MD 
simulations were able to provide good overall agreement with the 
experimental results [19]. However, their predictions of endohedral 
formation at high collision energies could not be confirmed by experiment. 
Indications for the possibility of endohedral formation in + Ceo collisions 
at high collision energies (velocity = 0.02 a.u.) were also found in 
nonadiabatic MD simulations [20] but such collision complexes would not 
survive intact for times long enough to be detected in the experiment. 

o 




Figure 4 . 4 . Relative cross sections for products of C*' + C^o collisions as a function of 
collision energy. Filled diamonds: CT (Ceo^), note the right hand scale; open circles: C^*', 
filled circles: Csg"^; open down triangle: C56"^; filled down triangle: C5T; open square: €52^; 
filled square: €50"^; open up triangle: C48+; filled up triangle: C4T. Adapted from [ 16 ]. 

Another interesting system that was studied is Fe"^ + Ceo [IV]. Two 
FeCeo^ components are formed at rather different collision energies that 
appear to have different stabilities although the fragmentation patterns argue 
against an endohedral product. The CID dissociation of a relatively strongly 
exohedrally bound Fe-Ceo complex formed by ligand exchange reactions 
with fullerenes had been investigated previously [21]. 

1.2 Charge Transfer 

Low energy charge transfer reactions are closely related to the addition 
reactions discussed in the previous section and often provide a competing 
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reaction channel. The most extensive low energy studies have been carried 
out by Bohme and collaborators [1]. Various experiments have shown that 
Ceo can be chemically ionised at thermal and higher energies in single- 




Number ofi3c 



Figure 4.5. Isotopomer distributions resulting from reaction of with Ceo- (a) at low 
collision energy and high octapole trapping voltage (b) at low collision energy with 
octapole set to discriminate against charge transfer products, (c) isotopomer distribution 
expected for no isotope scrambling (d) isotopomer distributions expected for full isotopic 
scrambling. Adapted from [16]. 

electron transfer reactions with cations that have an electron recombination 
energy that is greater than or equal to the ionisation potential of Ceo- Early 
observations involved charge transfer between atomic metal ions and Ceo 
[21,22]. Rate coefficients and branching fractions for charge transfer were 
not reported but Fe"^ was said to react with Ceo predominantly via charge 
transfer. Reactions of Ceo with other atomic ions having higher electron 
recombination energies than 7.6 eV e.g. (13.6 eV), Si"^ (8.151 eV) and the 
rare gas ions have been investigated in SIFT experiments at 294 K [23,24]. 
All of these ions were found to extract an electron from Ceo under the 
conditions of the experiment. 

Reactions of He"^ and Ne"^ with Ceo, studied in SIFT experiments, 
produced substantial amounts of doubly charged in electron transfer, 
electron detachment reactions of the type 

X^ + C,o^C+X + e- 



(4.1.6) 
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that competed with single electron transfer [23]. No dication formation was 
found for collisions with Ar"^. Several mechanisms have been proposed to 
describe removal of two electrons from Ceo in its reaction with He"^ and Ne"^ 
[23]. One possibility is that the high exothermicity of the transfer of one 
electron from Ceo to the rare gas cation leaves the fullerene cation 
sufficiently electronically excited to autoionise (or to undergo fast thermal 
ionisation before energy equilibration with the vibrational degrees of 
freedom, see next section). A second possibility that has been proposed 
involves an intial electron transfer to produce a metastable rare gas atom 
which then Penning ionises the Ceo^ at smaller separations. Another 
alternative mechanism views Ceo with its 240 valence electrons to have a 
metalloid character [1]. An Auger emission mechanism has long been used 
to explain the ejection of electrons from metals by slow rare gas ions. With 
Ceo this mechanism would involve the initial transfer of a deeply embedded 
electron followed by Auger emission of an electron. 

Extensive measurements have been carried out using SIFT to determine 
the rate constants for electron transfer in collisions with multiply charged 
fullerene ions [1], some of which are discussed in section 4.1.1 above. The 
absence of electron transfer to from molecules having ionisation 

energies in the range from 9.58 to 1 1.36 eV, still below the ionisation energy 
of Ceo^, is accounted for by the model in which a reaction barrier arises from 
the Coulomb repulsion between the singly charged product ions of the 
electron transfer process, see 4.1.1 above and Fig. 4.3. This assumes that 
electron transfer is most likely to occur when the reactants are in close 
proximity and one of the charges has been localised by the electrostatic 
interaction between the reactants with the other charge diametrically 
opposed on the fullerene surface. The requirement for electron transfer from 
a dication to a neutral molecule, producing two monocations initially in 
close proximity, is that the reaction exothermicity must exceed the 
electrostatic potential which denotes the Coulomb repulsion between the 
monocations at the moment of electron transfer (see Chap. 3.4.1). The 
Coulomb repulsion depends on the distance between the charges. If it is 
assumed that the charge separation at the moment of electron transfer is 
equal to or greater than the diameter of the fullerene cage then the 
electrostatic potential is < 2.06 eV [1]. This agrees well with the 
experimental observations. 

Single electron transfer reactions with triply charged Ceo^^ have also been 
studied for a wide range of molecular targets [25,26]. The ionisation energy 
threshold for the occurrence of electron transfer to Ceo^^ lies well below the 
thermodynamics threshold by about 4.4. eV. This is higher than the 
analogous differnce observed for dication reactions. Again, the explanation 
lies in the Coulombic repulsion since the products of single-electron transfer 
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are both charged. In this case we have one singly and one doubly charged 
ion initially in close proximity to one another. The threshold for production 
of such an ion pair will be significantly above the threshold or production of 
ions at infinite separation. A barrier of 4.5 eV has been estimated from a 
geometrical model for the charge distribution in and the effects of 
Coulomb repulsion between the products [27]. 

Double-electron transfer has been observed between and polycyclic 
aromatic hydrocarbons (PAHs) at room temperature [28]. Both sequential 
and concomitant electron transfer have been proposed as classical 
mechanisms of double-electron transfer to [28]. In the sequential 

mechanism, the transfer of a single electron is proposed to occur first, 
followed by reorganisation of the remaining two charges on Ceo and then the 
transfer of the second electron (see also 3.4.1 and Fig. 3.15). This 
mechanism assumes that the intial three charge sites on the cage are too far 
apart for two electrons to be transferred at the same time to two different 
charge sites. The second classical mechanism proposed involves the 
concomitant transfer of two electrons to form an intermediate zwitterion 



2. PENNING IONISATION 

Penning ionisation of Cm 

RG* -tCgo ^ RG -t C*,; -t e“ (4.2.1) 

has been observed in flow tube studies [29] and studied in more detail by 
two other groups using molecular beams [30-32]. The ionisation cross 
section was found to rise exponentially with increasing excitation energy. 
The experimental results from [30,31] are shown on Fig. 4.6. The values 
reported by Brunetti et al. [31] were only relative but they provide a point at 
low excitation energy (8.32 eV) for collisions with Xe*. The two data sets 
agree very well where identical collision partners were used so the data of 
Weber et al [30] were used to normalise the data from Brunetti et al. [31]. 
The data has been fitted by a model in which the energy transferred in the 
collision was assumed to be statistically distributed among the electronic 
degrees of freedom [33]. The ionisation was regarded as being thermal and 
the total ion yield was quenched on a time scale of the order of 240 fs by 
coupling to the thermal sink of the nuclear degrees of freedom of the 
fullerene. The coupling is characterised by a time constant which, in the 
original study [30] was fitted to be about a picosecond but with later, 
improved calculations of the electronic density of states this value was 
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Figure 4.6. Yield of ions in Penning ionisation experiments on C^o- The x-axis gives the 
energy transferred from each excited rare gas. Data from [30,31]- The full line is a fit to the 
data assuming that electrons are emitted thermally before the energy can be equilibrated with 

the degrees of freedom [33]. 



revised downwards [33]. The fit of the time constant requires an expression 
for the electron emission rate and, since the model is statistical, also an 
expression for the electronic level density of Ceo and Ceo^- The rate constant 
is given as [34,35] 



k{E,e)de- 



ttV p^(E) 



(4.2.2) 



where e is the kinetic energy of the electron, m,, the mass of the electron, p 
the level density of the daughter ion and parent molecule, E the internal 
energy of the molecule and a(e) the capture cross section for an electron in 
the Coulomb potential of the daughter ion. The total rate is then 

E 

k(E)^jk(E,e)de. (4.2.3) 

0 

The rate constant of eq. (1) is shown in Fig. 4.7 together with the rate 
constant expected when the same excitation energy is distributed over the 
vibrational degrees of freedom (see 3.2.1). 
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The model was later applied to understand the photoelectron spectra 
produced from the interaction of Ceo with ultrashort laser pulses [33,36]. 




E(eV) 

Figure 4.7. Rates for thermal electron emission from Qo as a function of internal energy. 
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Chapter 5 

COLLISION INDUCED DISSOCIATION 



1. SINGLY-CHARGED FULLERENE PROJECTILE 
IONS 

Collision-induced dissociation (CID) is a powerful technique to 
investigate the structure and stability of molecular ions. The fragmentation 
of fullerene ions induced in collisions with atomic and molecular targets has 
been widely studied since macroscopic amounts of fullerenes were first 
made available. In the first report [1] a reversed-geometry, double-focusing 
mass spectrometer, with a collision cell situated between the magnetic and 
electric sectors, was used to study the CID reactions of mass-selected Ceo^ 
and Cvo"^ at laboratory collision energies in the range 5-8 keV. Fragmentation 
mass spectra were reported down to C 4 o^ with mass peaks separated by 24 u 
(C 2 ). In the same paper, similar mass spectra for the CID of doubly charged 
fullerenes and also from charge inversion of negatively charged fullerenes in 
collisions with O 2 were reported. The mass spectra were very similar to 
those reported some years earlier by the Smalley group for the 
photodissociation of mass-selected fullerenes using ns laser pulses [2]. A 
number of other groups subsequently reported related results using a variety 
of target gases [3-9]. In these experiments, the laboratory energy was mainly 
fixed and the centre-of-mass energy was varied by changing the target gas. 
For sufficiently high excitation energies, a bimodal fragment pattern was 
observed [3,4,6,7]. The development of such a collisionally induced bimodal 
fragment pattern in shown in Fig. 5.1(a). These spectra were measured with 
a time-of-flight mass spectrometer and are for collisions between Ceo^ and 

85 




86 



Chapter 5 



Ar atoms [10]. The ion intensities have not been corrected for detection 
efficiency so that the relative intensity of the small fragment ions is 
exaggerated in this plot [11]. An extract of the fragment ion spectrum 
obtained for collisions with Ar at a centre-of-mass collision energy of 100 
eV is shown in Fig. 5.1(b). This very clearly shows the cross-over at n = 32 
from larger fullerene-like fragments with the characteristic C 2 spacing 
between the mass peaks, to the small, predominantly ring isomers, separated 
in mass by a carbon atom. The geometry of the various carbon cluster ions 
has been studied in beautiful gas ion chromatography experiments. The 
geometry can be determined from a measurement of the drift time of ions 
through a buffer gas. Fig. 5.2 shows the mobility measured for different 
carbon cluster ions [12]. The results can be divided into different families 
depending on the geometry. Fullerenes have a high mobility due to their 
compact structures, and can be observed down to n = 30. Many different 
isomers can co-exist for cluster sizes of around n =30 but a simple ring 
structure is the most likely isomer for 10 < n < 20. Below this, linear 
structures dominate. 




Figure 5.1. Fragmentation in € 50 ^ + Ar collisions, from [10]. (a) Development of the bimodal 
mass distribution with increasing collision energy, (b) extract of the fragment distribution at a 

collision energy of 100 eV. 

The bimodal fragment distribution has been explained in terms of 
successive statistical C 2 evaporation down to an unstable € 30 "^ that then 
breaks apart, or by a cleavage of the fullerene cage in which large cluster 
fragments are emitted. For low excitation energies it is now clear that the 
dominant fragmentation mechanism is successive evaporation of C 2 . As the 
excitation energy increases there is a small but increasing probability for the 
emission of larger neutral fragments such as C4 or Ce- The initial 
fragmentation process on the picosecond time scale can, however, if the 
excitation energy is high enough, lead to a catastrophic break-up of the 
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Number of carbon atoms 

Figure 5.2. Ion mobility versus cluster size for positively charged carbon cluster ions from n 

= 3 to n = 60. The lines through the data are to guide the eye for the four labeled families of 

isomers. Adapted from [12]. 

fullerene cage. Multiple numbers of large fragments result which can 
undergo a subsequent metastable fragmentation. For fullerene-like fragments 
this will be C 2 loss whereas for the smaller fragment ions the dominant 
metastable fragmentation pathway is via the loss of C 3 or a carbon atom. 
This scenario is supported by some molecular dynamics simulations [13-16] 
and has been shown explicitly for photodissociation experiments using 
nanosecond laser pulses [13,17]. Experimentally, the detection of large 
neutral fragments from collision experiments, under conditions where the 
positively charged fragments show a bimodal distribution, was reported by 
McHale et al. [18]. In these experiments, the neutral fragments were post- 
ionised in subsequent collisions with O 2 . A mass distribution closely 
mirroring that of the positively charged fragments was obtained and 
provided conclusive evidence for a cleavage mechanism producing large 
neutral fragments at these collision energies (340 eV in the centre of mass) 
rather than the dominance of C 2 evaporation. Fig. 5.3. The spectra of neutral 
fragments, determined in this way, as well as the positive ion fragment 
spectra from CID of are detected on a time scale of microseconds after 
the initial impact and the ions have time to undergo further metastable 
cooling evaporations. The detected mass distributions are thus strongly 
influenced by the thermodynamic stability of the neutral and charged carbon 
clusters as well as the ionisation potentials. This is the reason for the typical 
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appearance of the "magic numbers" in the mass spectra, in particular for n 
=7,11,15,19,23,50 [3,4,10,18,19]. The detected masses therefore do not 
necessarily provide direct information on the dynamics of the initial 
fragmentation step. 




m/q 



Figure 5.3. a), b) CID mass spectra produced in collisions between Qo"^ and He or O 2 
collisions at a laboratory collision energy of 8 keV. c) collisional dissociation /reionisation 
experiment where the fullerene is first dissociated in collisions with He. The charged 
fragment ions are removed from the beam and the neutral fragments are post-ionised in 
collisions with O 2 . Adapted from [18]. 

Similar fragmentation patterns are observed for the positively charged 
fragments produced in collisions with negatively charged ions [1,9,20,21], 
but the mass distribution of negatively charged fragments looks quite 
different. Fig. 5.4. Very few negatively charged fullerene fragments are 
detected. This is attributed to the very high probability of thermal electron 
emission leading to only neutral fragments that generally go undetected. 
Thermal electron emission from highly vibrationally excited fullerenes and 
metal clusters has been addressed in review articles [22,23]. 

Detailed collision energy dependencies of the fragmentation cross 
sections in low energy (5-300 eV centre of mass) collisions of positively 
charged fullerene ions with a variety of targets have been reported [10,24]. 
The form of the cross sections for production of large fullerene fragments 
was shown to have a predominant maximum that shifted to larger collision 
energies as the size of the fragment ions decreased and as the complexity of 
the target molecules increased. A combination of classical molecular 
dynamics simulations with statistical RRKM calculations to take into 
account the metastable fragmentation processes on the microsecond time 
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Figure 5.4. Time-of-flight mass spectrum of the negatively charged products of collisions 
between Cgg' and N 2 at a laboratory collision energy of 30 keV. Adapted from [20]. 

scale of the experiments was able to give reasonable qualitative agreement 
with the measured cross sections for collisions with Ar [10]. The comparison 
is shown in Fig. 5.5. The calculations took the total internal excitation 
energy into account. An average initial internal energy of the projectile ion 
was estimated (based on the amount of metastable fragmentation of the 
projectile beam, determined in the absence of collisions) and an estimate of 
the energy transferred in the collisions was obtained from classical MD 
simulations. The agreement is quite satisfactory considering that there is a 
large uncertainty in the absolute value of the experimental cross sections due 
to uncertainties in determining the target vapour pressure (see Chap. 2.1.3). 
The relative cross section values as a function of collision energy are much 
more accurate. Also there are uncertainties associated with the modelling. In 
particular, the width and form of the initial internal energy distribution in the 
projectile beam is now known to play a crucial role for determining the 
appearance of fragmentation spectra, as shown in laser ionisation / 
fragmentation studies [25]. We do not know the initial internal energy 
distribution in the fullerene ion beams and can only make a rough estimate 
of the average internal energy. Also, the parameters used for the RRKM 
modelling are different from the parameters believed today to provide the 
best description of the metastable decay of fullerenes (Chap. 3.2.1), 
however, under the conditions of the laser experiments, this should not make 
a major difference. The pronounced maxima in the cross sections were seen 
to be due predominantly to the metastable evaporation of C 2 with the actual 
shape of the maximum depending on the experimental time window for 
detection of the fragment ions. The initial internal energy of the projectile 
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Figure 5.5. Collision energy dependence of the absolute fragmentation cross section in C^o* + 
Ar collisions, from [10]. Top: experimental data. Bottom: results from cacluations in which 
MD simulations are combined with statistical RRKM calculations. 

ion does indeed strongly affect the measured form of the cross sections for 
production of the large fullerene fragment ions. Fig. 5.6 shows a comparison 
of the total CID cross sections determined for collisions between Ceo^ and 
different targets. Two sets of collisions with Ar are reported where the 
fullerene projectile ion beam was produced via laser desorption/ionisation 
with different laser fluence [26]. The cross section labelled Ar(a) was 
measured using a laser fluence of approximately 20 mJ/cm^ (as used for the 
data in Fig. 5.5) whereas Ar(b) was obtained with a laser fluence of 
approximately 5 mJ/cm^ which is just slightly above the threshold for 
production of Ceo^- From measurements of black-body radiation from laser 
desorbed fullerenes [27] it is known that the internal energy of the desorbed 
molecules scales with the laser fluence. By applying the same simple 
MD/RRKM modeling to the Ar data of Fig. 5.6 as used to fit the data in Fig. 
5.5, it could be estimated that the initial average internal energy of was 
ca. 30 eV in the case of the Ar(a) cross section and ca. 15-20 eV in the case 
of Ar(b). The lower absolute value of the measured total fragmentation cross 
section for Ar(b) compared to Ar(a), as well as its shift to higher collision 
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Figure 5.6. (a) Total fragmentation cross sections measured for collisions of Co* with Ar, Xe 
and SFs as a function of the centre of mass collision energy. The laser fluence used for 
production of the fullerene projectile ion beam was 20 mj/cm^ for all data except for Ar(b) 
which was 5 mj/cm^. (b) Sum of the individual fragmentation cross sections for CC 
production for n < 52. From [26]. 

energies, is due to a lower total excitation energy of the Ceo^ and the fact that 
metastable C 2 evaporation which, according to the model, is assumed to be 
the major contributor to the large peak, can occur on a time scale of more 
than a few microseconds which was the limit of the experiments. In other 
words, the lower total excitation energy (initial internal energy plus energy 
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transferred in the collision) for Ar(b) leads to a smaller rate constant for 
metastable C 2 loss and thus less metastable fragmentation will be detected on 
the microsecond time scale of the experiment [26]. The classical MD 
simulations indicate that high impact parameter, grazing collisions lead to a 
small energy transfer to the cage [28]. This is apparently sufficient to induce 
substantial metastable C 2 loss from the more highly excited projectile ion but 
not for the colder ones. The Xe and SFe cross sections, also shown on Fig. 
5.6, were measured under laser desorption conditions corresponding to the 
internally hot projectile beam. The total cross sections for these targets do 
not show a pronounced maximum below a centre of mass collision energy of 
100 eV but are much broader. The large difference in magnitude of the SFe 
and Xe cross sections is mainly due to the intensity of the large fullerene-like 
fragments. This can be seen in Fig. 5.6(b) where the sum of the individual 
fragment ions up to and including € 52 ^ has been plotted. There is still a 
residual effect of metastable C 2 loss in the Ar(a) cross section around 50 eV 
but it has been considerably reduced compared to the total cross section 
results. The summed cross sections for the other targets are all very similar. 
Although the Xe and SFe results were obtained with projectile ions produced 
under the same conditions as used in the Ar(a) experiments i.e. with an 
initial internal energy of around 30 eV, the total cross section behaviour is 
much more similar to the Ar(b) measurement. More information can be 
obtained from plotting the sum of fullerene-like fragments separately from 
the small fragment ions. This has been done in Fig. 5.7. Although the 
fullerene-like fragment intensities are smaller for collisions with Xe, the 
small fragment ions are higher in intensity. The reasons for this are not yet 
fully understood. It should also be noted that for the laboratory collision 
energy range investigated there will be significant scattering of the fragment 
ions in the laboratory frame for collisions with the heavy partners (Xe and 
SFe) and signal intensity will be lost due to this. The effect will be most 
pronounced for the smaller fullerene-like fragments [29]. This is the reason 
for the exceptionally low cross section for fragment ions in the range 32 < n 
< 52 measured for Xe collisions (Fig. 5.7(a)). The slight increase observed 
for these cross sections at centre of mass collision energies beyond 200 eV is 
a consequence of the smaller laboratory scattering angles at these higher 
energies. It is more difficult to understand the higher cross sections 
measured for the small fragment ion production. It may be that the stronger 
interaction potential simply leads to an overall higher collision cross section. 
Alternatively, it may have to do with the details of the elastic collisions with 
the carbon cage that lead to an overall higher energy transfer in this collision 
energy range for the larger Xe atom (see e.g. Fig. 3.3) and therefore a higher 
tendency for prompt cage break-up to form bimodal fragment distributions 
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Figure 5.7. Sum of the individual fragmentation cross sections for the data plotted in Fig. 5.6 
for (a) 32 < n < 52 (b) n < 26. From [26] . 

Some interesting details can be extracted from the cross sections in Fig. 
5.7. The fullerene-like fragments from collisions with both Ar(a) and Xe, 
which were obtained with hot (ca. 30 eV) projectile ions have maxima at the 
same collision energy. However, the onset for the Xe collisions is slightly 
shifted to higher collision energies. This was also remarked on in connection 
with earlier experiments [10] where the appearance energies of the 
fragments were discussed. These are plotted in Fig. 5.8. The appearance 
energies measured in collisions with Ar and Ne are very similar over the 
entire fragment ion mass range. In contrast, the appearance energies in Xe 
collisions are shifted to higher energies with the SFe thresholds lying even 
higher for loss of more than 10 carbon atoms. The appearance energies for 
Ar and Ne are in very good agreement with maximal entropy calculations 
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Figure 5.8. Appearance energies of the various fragment ions in collisions between Ceo^ and 
rare gas atoms or SFs molecules. Filled symbols are experimental data from [10]. Open 
triangles are from maximal entropy calculations [25]. 

assuming an initial internal energy of 30 eV and total transfer of the centre 
of mass energy to internal degrees of freedom. These results are also shown 
on Fig. 5.8. Note that the calculated appearance energies differ from the 
original publication [30] and have been obtained using a corrected energy 
scale [25]. Originally, with the incorrect energy scale in the calculations, 
better agreement was obtained between the maximal entropy model and Xe 
collisions [30]. The lower appearance energies measured for Ar and Ne were 
therefore attributed to impulsive "knock-out" of fragments from the cage 
without the transfer of large amounts of energy to vibrational degrees of 
freedom. Some evidence for such processes is discussed below but we now 
believe that this is not the main cause of the discrepancy between the Ar, Ne 
and Xe results in Fig. 5.8. Instead, we relate it to the greater tendency for 
transient endohedral fullerene formation in the collisions with the smaller 
rare gas atoms (see Chap. 8). This provides a very efficient mechanism for 
transfer of the entire centre-of-mass collision energy to internal degrees of 
freedom of the cage for those collision trajectories that lead to capture. 
Endohedral fullerene formation is much less efficient for collisions with Xe 
and is expected to occur only for collision energies higher than most of the 
fragment appearance energies. Therefore there is a much greater tendency 
for elastic scattering with cage components leading to transfer of energy less 
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than the centre of mass value. At the centre-of-mass collision energies of 
interest, the laboratory energy of the fullerene projectile in Xe collisions 
only covers the range up to slightly less than 2000 eV. For these energies, 
according to the simple model discussed in Chap. 3, we would expect the Xe 
to interact with a large region of the cage surface leading to a rather low 
transfer of energy due to elastic scattering with the cage surface (Fig. 3.3). 
The even higher appearance energies for collisions with SFe can be 
rationalised by considering the partition of energy between the internal 
degrees of freedom of the fullerene projectile and the internal degrees of 
freedom of the molecular target. The drop of the threshold energies at n = 30 
for SFe is due to the formation of doubly charged fullerene ions in these 
collisions and the appearance energies plotted in Fig. 5.8 are actually for 
Cjn'^ in the region between n = 20 and 30. 

The onset of the fullerene-like fragments and the position of the total 
cross section maximum are shifted by at least 10 eV to higher collision 
energies for Ar(b) compared to Ar(a) (Fig. 5.5a). This corresponds well to 
the difference in internal energy of the projectile ions in the two series of 
measurements. Th SFe cross section for fullerene-like fragments is very 
similar to the Ar(b) cross section. This is mainly a consequence of the 
statistical distribution of excitation energy among all vibrational degrees of 
freedom in the collision system leading to less vibrational excitation of 
for a given collisional energy loss (remember that the fullerene ions in the 
collision experiments with Xe and SFs have a high initial internal energy of 
ca. 30 eV). At collision energies around the maxima of the cross sections this 
energy difference is roughly the same as the difference in internal excitation 
between the fullerenes used in collisions Ar(a) and Ar(b). 

There is a significant difference in the small carbon fragment ion mass 
distributions measured for Xe and SFe collisions. A comparison is shown in 
Fig. 5.9. The distribution is clearly shifted towards the low mass range for 
SFs. There are also additional mass peaks in the SFg spectrum, doubly 
charged fullerene ions and CnF"^ that are reaction products of the collision 
[26]. There is no evidence for corresponding reaction products with 
fullerene-like fragment ions. The enhanced formation of singly-charged 
small fragments may be explained by the ionisation of small neutral carbon 
cluster fragments in the course of the collision. Fragments produced on 
collision induced dissociation of the SFe target have high electron affinities 
e.g. EA (F) = 3.4 eV and EA(SF5)=3.7 eV [31], much higher than the parent 
molecules, EA(SE6) = 1.05 eV. These fragments can strip an electron from 
the neutral carbon fragments during the course of the collision. The 
formation of CnE"^ competes with the charge transfer reaction [26]. This 
implies that the carbon cage must break-up on a very short time scale (a few 
hundred fs) in order for the small neutral carbon fragments to react with the 
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fragmentation products from the target. This time scale for cage break-up is 
confirmed by MD simulations and also by surface collision experiments 
where the neutral fragments can pick-up an electron from the surface (see 
Chap. 10). 




Figure 5.9. Time-of-flight mass spectra of products from collisions between Ceo^ and Xe (top) 
and SFj (bottom) at a centre of mass collision energy of 236 eV. Only the small fragment 
range of the bimodal fragment distribution is shown. The arrows indicate CnF"^, n = 1-11. 

From [26]. 

The results of the low energy collision experiments with fullerene 
projectile ions can be well explained within the context of statistical models. 
However, for very small target atoms there is some evidence for an 
impulsive "knock-out" mechanism. Classical MD simulations have shown 
trajectories where individual atoms are "knocked out" of the cage in elastic 
collisions [32]. Experimental evidence for this kind of mechanism has been 
found for collisions at 50 keV (lab) with He atoms but not for collisions with 
Ne or Ar targets [6]. The detection of a small signal due to a fragment 
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ion in collisions with He, shown in Fig. 5.10, was taken as providing 
evidence for the impulsive "knock out" of a cage atom without the transfer 
of large amounts of energy to vibrational degrees of freedom. Non-adiabatic 
MD simulations also showed trajectories corresponding to such a process but 
for collisions with fast Ar ions [33]. However, due to the large amount of 
electronic energy transferred in these collisions, it would not be possible to 
directly observe the remaining € 59 "^ since it would very rapidly undergo 
unimolecular decay to smaller and more stable ions. 
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Figure 5.10. Positive ion fragmentation spectra around a mass to charge ratio of 59 for C^o 
collisions with He, Ne, Ar at a laboratory collision energy of 50 keV. 059 "^ is detected for 
collisions with He but not with Ne or Ar. Adapted from [ 6 ]. 

Large differences were found in the fragmentation spectra measured for 
different target gases (all at a laboratory collision energy of 50 keV). The 
differences could not be explained simply in terms of the different total 
energy available for excitation (centre of mass energy). The energy transfer 
to the cage was modeled using the Lindhard model, described in Chap. 3.1.1. 
The experiments fall into the nuclear stopping regime and the model was 
used to estimate the energy transferred in each collision as a function of 
impact parameter. They then considered the probability for a "prompt" 
fragmentation where single carbon atoms can experience an energy transfer 
in a single collision that is larger than the energy required for bond breaking 
all the bonds of the atom in the cage. The rest of the collisionally transferred 
energy was considered to be statistically distributed among the vibrational 
degrees of freedom of the remaining cage atoms leading to the normal 
statistical fragmentation mechanisms discussed above. Using this simple 
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model they were able to obtain very good qualitative agreement with the 
experimentally measured fragment ion distributions [6]. 



2. FULLERENE TARGETS 

Collision experiments between atomic singly charged ions and neutral 
fullerenes were carried out using a guided ion beam, apparatus [34-37]. The 
results of these studies will be discussed in more detail in Chap. 8, with 
reference to the collisional form a tion of endohedral fullerenes. The centre of 
mass collision energies were in the range 20-100 eV, comparable to those 
discussed above and shown e.g. in Fig. 5.5. Unfortunately, for experimental 
reasons in guided ion beam experiments it was only possible to determine 
the cross sections for the production of the large fullerene-like fragment 
ions. The experiment was blind to the small (n < 32) fragment ions. The 
cross section behaviour was very similar to that discussed above with a clear 
maximum that shifts to higher collision energies as the size of the fragment 
ion decreases. Although the time scale for metastable fragmentation os 
orders of magnitude longer in the guided ion beam experiments compared to 
time-of-flight mass spectrometry, the relative behaviour of the cross sections 
is rather similar. This is a consequence of the influence of radiative cooling 
on the hot fragment ions [38]. The internal energy is cooled down to below 
the fragmentation threshold on a time scale of a few tens of microseconds. 

Many experiments have been carried out using highly charged projectile 
ions colliding with neutral fullerenes. The charge transfer processes in such 
collisions are described in Chap. 7. Here we will just focus on the aspects of 
the experiments and modeling that provide insight into the energy transfer 
and fragmentation behaviour of the fullerenes. 

Experiments involving high energy collisions with singly or multiply 
charged ions make it possible to scan the impact velocity range where one 
expects to see a change in the nature of the excitation mechanism. For 
relatively low energy collisions, energy is predominantly transferred via 
vibrational excitation (nuclear stopping) while for high energy collisions it is 
electronic energy transfer (electronic stopping) that is expected to dominate 
(Chap. 3.1.1). It may also be expected that direct fragmentation processes 
can be observed for high energy collisions where fragmentation occurs 
before the energy can be statistically distributed among the vibrational 
degrees of freedom. Many groups have looked for fingerprints of non- 
statistical behaviour. 

A very spectacular change in fragmentation behaviour as a function of 
impact velocity was observed by Schlathdlter et al. in collisions between 
He‘*'^ projectile beams and neutral Ceo [39,40]. Mass spectra recorded for 
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collisions with He"^ are shown in Fig. 5.11. At projectile velocities below 
0.45 a.u. (corresponding to a collision energy of 20 keV in Fig. 3.6) there is 
a significant amount of large fullerene fragments indicating the statistical 
evaporation process seen at lower collision energies and discussed above. 
These fullerene-like fragments decrease with increasing collision velocity 
and have almost vanished by 0.45 a.u. while the yield of small fragment ions 
increases strongly beyond this point. The velocity of 0.45 a.u. is close to 
where the electronic excitation is predicted to start to dominate over 
vibrational excitation, according to the Lindhard model. The experiments 
were nicely explained by invoking the concepts of the Lindhard model but 
using a more elaborate modeling procedure [40]. Relatively little energy was 
shown to be transferred in the elastic collisions between the projectile and 
target atoms at the high collision velocities of these experiments. It is 
sufficient to produce evaporation of C 2 molecules to form large fullerene- 
like fragments but is not sufficient to lead to the phase transition and 
bimodal fragment distribution discussed above. The electronic excitations, 
due to the electronic stopping of the projectile moving through the electron 
gas of the fullerene, that dominate for v > 0.45 a.u., are more efficient at 
transferring energy to the fullerene. The simulated projectile energy losses 
are on the order of 100 eV [40] and thus sufficient after energy coupling to 
vibrational degrees of freedom has occurred to cause break-up of the cage 
with emission of ring and chain fragments, leading to the mass distributions 
shown in Fig. 5.11. The observation of multiply-charged fullerene ions is 
also in agreement with the occurrence of strong electronic excitations, even 
for the lowest velocity studied. This interpretation of the experimental results 
has been qualitatively confirmed by nonadiabatic MD simulations, discussed 
in Chap. 3.3.3 [33]. In this approach, electronic and vibrational degrees of 
freedom are treated simultaneously and self-consistently by combining time- 
dependent density functional theory with classical MD. Snapshots of 
trajectories for collisions between Ar"^ and Ceo are shown in Fig. 3.12. The 
low velocity (0.02 a.u.) trajectory leads to immediate strong deformation of 
the cage corresponding to the nuclear stopping regime. For the high velocity 
(0.45 a.u.) trajectory, the ion seems to pass straight through the fullerene 
cage without causing any significant deformation but transfers a substantial 
amount of energy to electronic excitation. The electronic excitation relaxes 
on a time scale of a few hundred femtoseconds [41] thereby exciting 
vibrational degrees of freedom and ultimately causing fragmentation. There 
is no reason, from the MD simulations, to invoke exotic fragmentation 
mechanisms to explain the appearance of the mass spectra. The strongest 
experimental signature of the crossover between nuclear and electronic 
stopping should be found in angular distributions of the fragments. Such 
experiments have not yet been reported. 
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Figure 5.11. Mass spectra obtained for He"^ + C^o collisions [39,40], At slow projectile 
velocities ( < 0.45 a.u.), fragmentation is dominated by statistical evaporation of C 2 units. 
With increasing projectile velocity the yield of small fragment ions increases while that of the 
large fullerene-like fragments decreases. 

Calculations of the energy transfer versus collision velocity are shown in 
Fig. 5.12 for three different projectile ions and compared with TRIM 
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Figure 5.12. left: results of nonadiabatic MD simulations for the total kinetic energy loss 
(open circles) and its electronic contribution (open triangles) as a function of collision 
velocity for collisions between different ions (ff^, C*, Ar’^) with Qo and a fixed collision 
geometry with b = 2 a.u. Experimental data for + Cso are shown by filled circles [42]. 
Right: TRIM calculations [43] for the same collision systems. Solid lines: total kinetic energy 
loss. Dashed lines: electronic contributions. Adapted from [33]. 

calculations [33]. There are also some data points obtained for proton 
collisions [42] shown on the plot. The total energy transfer and the 
component due to electronic excitation alone are shown. Some important 
differences can be seen between the nonadiabatic MD calculations and the 
TRIM code. There are large differences in the absolute values of the energy 
transfer for the different ions. This will result in qualitatively different 
fragment distributions. More surprisingly, the calculations provide velocity- 
independent excitation energies with dominating electronic excitation in the 
high velocity range (0.25-0.5 a.u.) for all three ions. These two results 
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provide an explanation for the outcome of collisions with different ions in 
the high velocity range (v > 0.4 a.u.) where very different but velocity 
independent fragmentation patterns were seen for the different collision 
partners [44]. Some examples of mass spectra are shown in Fig. 5.13. A 
velocity-independent electronic energy transfer deviates from the predictions 
of the electronic stopping power model that shows an almost linearly 
increasing electronic stopping as the collision velocity increases. The 
saturation could be caused by the 2-dimensional electronic density 
distribution on the fullerene surface instead of a three-dimensional 
homogeneous electron gas as used for the derivation of the stopping power 
[43,45]. Finally, the transferred energy only increases smoothly with the 
velocity in the low velocity range (v < 0.3 a.u.) for proton projectiles 
whereas the heavier projectiles have maxima in the transferred energy. This 
results in an increasing and then decreasing vibrational energy transfer [33]. 
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Figure 5.13. Time-of-flight spectra of ions produced in collisions with H’^, He"^ (200 keV 
laboratory collision energy) and Ar^"^ (600 keV laboratory collision energy). Adapted from 

[44]. 
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3. COINCIDENCE EXPERIMENTS 

Coincidence experiments have provided more information on the 
fragmentation dynamics. Often these experiments have involved higher 
projectile charge states so that it is not always possible to directly compare 
the results with the low energy, singly-charged fullerenes but they can be 
compared with the energetic atomic ion - fullerene collisions discussed in 
5.2, above. An early study, carried out for collisions between doubly charged 
and He at a collision energy of 120 keV showed that Cn"^ fragments as 
large as n = 9 could be measured in coincidence with heavy fullerene-like 
fragments [46]. Fig. 5.14 shows the coincidence signal obtained in these 
experiments. Coincidence experiments have also been reported for collisions 
between neutral fullerenes and highly charged ions [47,48]. In these 
experiments, the observed fragmentation behaviour was strongly dependent 
on the charge state of the fullerene after collision as well as the amount of 
energy transferred to internal degrees of freedom. Experiments in which 
fragments from 75 keV collisions between Ceo and H 2 were measured in 
coincidence are particularly revealing [49]. The apparatus used for these 
studies is shown in Fig. 2.13. The experiments showed that the emission of 
species larger than C 2 is competitive with C 2 emission in events that leave a 
large fullerene fragment. The measurements of small fragment ions obtained 
in coincidence with Cso"^ is shown in Fig. 5.15. Small fragments with both 
odd and even numbers of carbon atoms are found in coincidence with 
fullerene-like ions, indicating subsequent unimolecular decay of the light 
fragment if the initial step is a binary fragmentation. The emission of small 
odd numbered neutral clusters (C3) is the most energetically favourable 
unimolecular decay path for the small fragment ions [13]. The sum of the 
number of carbon atoms of the two species detected in coincidence is 
generally less than 60 showing that an initial binary fragmentation (or 
higher) must be accompanied by some neutral fragments. The experimental 
results were modeled assuming an initial binary fragmentation with 
subsequent fragmentation estimated according to RRKM theory (similar to 
[10], discussed above) [50]. The important features and correlations in the 
fragmentation yields measured in the experiments could be well reproduced 
assuming an exponential distribution of excitation energies. Such an 
excitation energy distribution is in agreement with the results of MD 
simulations for H 2 [51] and He [32] collisions for similar centre of mass 
collision energies. A different model was put forward to explain the 
fragmentation pattern observed in collisions between positive fullerene ions 
and H 2 at 20-200 keV laboratory collision energy [5,52]. here the authors 
assumed that multiple collisions between the target and carbon atoms in the 
cage could occur, inducing a large rift in the cage structure. Successive C 2 
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evaporation from the damaged region could then occur befor the rift "re- 
closes". 




Figure 5.14. Number of coincidences with large fullerene-like clusters registered as a function 
of cluster size for the products of + He collisions at a laboratory energy of 120 keV. 

Adapted from [46]. 
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Figure 5.15. Yield of small carbon cluster ions detected in coincidence with Cso* for 
collisions between Ceo and H 2 at a laboratory collision energy of 75 keV. From [49]. 

A multi -coincidence set-up (Fig. 2.14) for investigating the ionisation 
and fragmentation of Ceo in collisions with various fast (100-300 keV) low- 
charged projectile ions provided evidence for both bimodal "fission" and 
multifragmentation into three or more small charged fragments accompanied 
by neutral fragments [44]. Examples of time-of -flight fragment spectra in 
which 1 to 7 particles are detected in coincidence after collisions with Ar^"^ 
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are shown in Fig. 5.16 [53]. Each /-particle spectrum contains all the 
detected particles that came from events in which exactly i fragments were 
detected in coincidence. The angular distributions of the charged fragments 
determined in the same experimental setup were isotropic [53] i.e. no 
evidence was found for non-statistical "knock-out" or other direct 
fragmentation processes. 
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Figure 5.16. Time-of-flight spectra for 1 to 7 coincidentally measured particles after 
collisions between Ar^’^ and Qo at a laboratory collision energy of 300 eV. Adapted from 

153]. 

Evidence for multi-fragmentation has also been found in measurements 
of the kinetic energies of fragment ions produced in 2 MeV collisions 
between Si"^"^ and Ceo [54]. 

Detailed modeling of -i- Ceo collisions for the collision energy range of 
50-300 keV was carried out using a combination of models but again 
assuming statistical processes [55]. The ions obtained as the products of such 
collisions are predominantly fullerene-like fragment ions i.e. the intensity of 
small fragment ions from a prompt cage break-up is small (5-10%), as can 
be expected from the calculations of energy transfer in the collisions (Eig. 
5.12). At the very high collision velocities considered in these experiments 
(1-3.5 a.u.), the electronic stopping power actually decreases again with 
increasing velocity [55]. The coincidence experiments show that these small 




106 



Chapter 5 



fragment ions are predominantly from the fission of multiply charged 
fullerenes [55]. Typical mass spectra for the proton collisions are shown in 
Fig. 5.17. Three steps were considered in the calculations: firstly the energy 
transferred to the neutral fullerene target was modeled using the Lindhard 
model (Chap 3), secondly the electron emission from the excited neutral 
target was treated using the statistical energy deposition model [56] and 
thirdly the statistical C 2 emission was considered. A comparison of the 
experimental and calculated evaporation fraction is shown in Fig. 5.18 as a 
function of collision velocity. The evaporation fraction is defined as 
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Figure 5.17. Part of the time-of-flight spectra of positively charged products of collisions 
between protons and C^o at a laboratory collision energy of 300 keV. Adapted from [55]. 

Experimental results obtained by Itoh et al. for even higher collision 
velocities [57] are shown for comparison on the same plot. The two different 
experimental data sets fit well together and both sets of results are very well 
reproduced by the calculations. In order to obtain the agreement with 
experiment, the parameters used to model the thermal electron emission 
were adjusted in such as way as to decrease the energy required for electron 
emission while keeping the exponential pre-factor at the value generally 
used to model thermionic emission from vibrationally excited fullerenes. 
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The same effect could be obtained by assuming that fast thermal electron 
emission occurs from the hot electronic subsystem of Ceo before the energy 
is equilibrated with the vibrational degrees of freedom. Such a picture has 
been used to successfully model the electron emission from fullerenes 
excited with femtosecond laser pulses or in Penning ionisation experiments 
(Chap. 4) [41]. The thermal ionisation was found to be an important 
parameter for correctly describing the relationship between evaporation 
fraction and fullerene charge state. The higher fragmentation for higher 
charge states has been observed in many experiments including fs laser 
excitation and electron impact ionisation as well as charge stripping in 
collisions with highly charged ions. 




projectile velocity [a.u.] 



Figure 5.18. Evaporation fractions for charge states q = 1 (squares), q = 2 (circles) and q = 3 
(triangles) as a function of the projectile velocity. Solid symbols and crosses are from [55], 
hollow symbols are from [57], Lines are the results of model calculations, discussed in the 

text. Adapted from [55]. 



REFERENCES 

1 A. B. Young, L. M. Cousins, and A. G. Harrison, Rap. Comm in Mass Spec. 5, (1991) 
226. 

^ S. C. O'Brien, J. R. Heath, R. F. Curl, and R. E. Smalley, J. Chem. Phys. 88, (1988) 220. 

3 K. A. Caldwell, D. E. Giblin, and M. L. Gross, J. Am. chem. Soc. 114, (1992) 3743. 

^ R. J. Doyle and M. M. Ross, J. Phys. Chem. 95, (1991) 4954. 

^ P. Hvelplund, L. H. Anderson, H. K. Haugen, J. Lindhard, D. C. Lorents, R. Malhotra, and 
R. Ruoff, Phys. Rev. Lett. 69, (1992) 1915. 





108 



Chapter 5 



^ M. C. Larsen, P. Hvelplund, M. O. Larsson, and H. Shen, Eur. Phys. J. D 5, (1999) 283. 

^ S. W. McElvany and J. H. Callahan, J. Phys. Chem. 95, (1991) 6186. 

^ T. Weiske, J. Hrusak, D. K. Bdhme, and H. Schwarz, Chem. Phys. Lett. 186, (1991) 459. 

^ D. H. Yu, L. H. Andersen, C. Brink, P. Hvelplund, D. C. Lorents, and R. Ruoff, Mol. Mat. 
4,(1994) 237. 

R. Ehlich, M. Westerburg, and E. E. B. Campbell, J. Chem. Phys. 104, (1996) 1900. 

E. E. B. Campbell, G. Ulmer, B. Hasselberger, H.-G. Busmann, and I. V. Hertel, J. Chem. 
Phys. 93, (1990) 6900. 

1^ G. V. Helden, M.-T. Hsu, P. R. Kemper, and M. T. Bowers, J. Chem. Phys. 95, (1991) 
3835. 

13 H. Hohmann, R. Ehlich, S. Furrer, O. Kittelmann, J. Ringling, and E. E. B. Campbell, Z. 
Phys. 0 33,(1995) 143. 

14 E. Kim, Y. H. Lee, and J. Y. Lee, Phys. Rev. B 48, (1993) 18230. 

13 S. G. Kim and D. Tomanek, Phys. Rev. Lett. 72, (1994) 2418. 

16 Y. Xia, Y. Xing, C. Tan, and L. Mai, Phys. Rev. B 52, (1995) 110. 

1"^ K. R. Lykke, Phys. Rev. A 52, (1995). 

1^ K. J. McHale, M. J. Polce, and C. Wesdemiotis, J. Mass Spectr. 30, (1995) 33. 

1^ E. E. B. Campbell, in Clusters of Atoms and Molecules, edited by H. Haberland (Springer, 
1992). 

30 D. Mathur, C. Brink, P. Hvelplund, N. Jensen, and D. H. Yu, Rap. Comm, in Mass Spec. 
9,(1994) 114. 

31 H. Shen, C. Brink, P. Hvelplund, and M. O. Larsson, Z. Phys. D 40, (1996) 371. 

33 E. E. B. Campbell and R. D. Levine, Annu. Rev. Phys. Chem. 51, (2000) 65. 

33 J. U. Andersen, E. C. Bonderup, and K. Hansen, J. Phys. B 35, (2001) Rl. 

34 E. E. B. Campbell, R. Ehlich, M. Westerburg, and I. V. Hertel, , Aarhus, 1993 (AIP Press), 
p. 697. 

33 A. Lassesson, A. Gromov, K. Mehlig, A. Taninaka, H. Shinohara, and E. E. B. Campbell, 
J. Chem. Phys. 119, (2003) 5591. 

36 R. Ehlich, H. Sprang, M. Westerburg, and E. E. B. Campbell, J. Chem. Phys. 108, (1998) 
9390. 

33 R. Mitzner and E. E. B. Campbell, J. Chem. Phys. 103 No. 7, (1995). 

38 R. Ehlich, E. E. B. Campbell, O. Knospe, and R. Schmidt, Z. Phys. D 28, (1993) 153. 

30 A. V. Glotov and E. E. B. Campbell, Phys. Rev. A 62, (2000) 033202. 

30 E. E. B. Campbell, T. Raz, and R. D. Levine, Chem. Phys. Lett. 253, (1996) 261. 

31 S. G. Lias, J. E. Bartmess, J. F. Liebman, J. L. Holmes, R. D. Levin, and W. G. Mallard, J. 
Phys. Chem. Ref. Data Suppl. 17, (1988) 872. 

33 R. Ehlich, O. Knospe, and R. Schmidt, J. Phys. B 30, (1997) 5429. 

33 T. Kunert and R. Schmidt, Phys. Rev. Lett. 86, (2001) 5258. 

34 J. E. Christian, Z. Wan, and S. L. Anderson, J. Chem. Phys. 99, (1993) 3468. 

33 z. Wan, J. F. Christian, and S. L. Anderson, Phys. Rev. Lett. 69, (1992) 1352. 

36 z. Wan, J. F. Christian, and S. L. Anderson, J. Chem. Phys. 96, (1992) 3344. 

33 z. Wan, J. F. Christian, Y. Basir, and S. L. Anderson, J. Chem. Phys. 99, (1993) 5858. 

38 K. Hansen and E. E. B. Campbell, J. Chem. Phys. 104 No. 13, (1996). 

30 T. Schlatholter, O. Hadjar, R. Hoekstra, and R. Morgenstern, Phys. Rev. Lett. 82, (1999) 
73. 

40 T. Schlatholter, O. Hadjar, J. Manske, R. Hoekstra, and R. Morgenstern, Int. J. Mass 
Spectr. 192,(1999) 245. 

41 K. Hansen, K. Hoffmann, and E. E. B. Campbell, J. Chem. Phys. 119, (2003) 2513. 

43 J. Opitz, H. Lebius, S. Tomita, B. A. Huber, P. Moretto-Capelle, D. Bordenave- 
Montesquieu, A. Bordenave-Montesquieu, A. Reinkoster, U. Werner, H. O. Lutz, A. Niehaus, 




5 . Collision Induced Dissociation 



109 



M. Benndorf, K. Haghighat, H. T. Schmidt, and H. Cederquist, Phys. Rev. A 62, (2000) 
022705. 

J. F. Ziegler and J. P. Biersack, TRIM-code, (1991). 

A. Reinkoster, U. Werner, and H. O. Lutz, Europhys. Lett. 43, (1998) 653. 

J. Lindhard, V. Nielsen, and M. Scharff, Mat. Fys. Medd. Dan. Vid. Selsk. 36, (1968) 1. 

C. Brink, L. FI. Andersen, P. Hvelplund, and D. H. Yu, Z. Phys. D 29, (1994) 45. 

L. Chen, J. Bernard, A. Denis, J. Desesquelles, and S. Martin, Physica ScriptaTSO, (1999) 
209. 

S. Martin, L. Chen, A. Denis, and J. Desesquelles, Phys. Rev. A 57, (1998) 4518. 

R. Vandenbosch, B. P. Flenry, C. Cooper, M. L. Gardel, J. F. Liang, and D. I. Will, Phys. 
Rev. Lett. 81,(1998) 1821. 

R. Vandenbosch, Phys. Rev. A 59, (1999) 3584. 

G. Seifert and J. Schulte, Physics Letters A 188, (1994) 365. 

D. H. Yu, L. H. Anderson, C. Brink, and P. Hvelplund, Z. Phys. D 29, (1994) 53. 

A. Reinkoster, PhD Thesis, (2001) Universtat Bielefeld. 

A. Itoh, H. Tsuchida, K. Miyabe, M. Imai, and N. Imanishi, Nucl. Instr. Meth. Phys. Res. 
B 129,(1997) 363. 

A. Reinkoster, U. Werner, N. M. Kabachnik, and H. O. Lutz, Phys. Rev. A. 64, (2001) 
023201. 

N. M. Kabachnik, V. N. Kondratyev, Z. Roller-Lutz, and H. O. Lutz, Phys. Rev. A 56, 
(1998) 2848. 

H. Tsuchida, A. Itoh, Y. Nakai, K. Miyabe, and N. Imanishi, J. Phys. B 32, (1998) 5289. 




Chapter 6 

CHARGE TRANSFER IN COLLISIONS WITH 
HIGHLY CHARGED IONS 



One of the presently most active research fields involving fullerene 
collisions is that of collisions between neutral fullerene targets and highly 
charged ions. This allows the study of processes in which a large number of 
electrons is captured from the fullerene targets under the action of the high 
potential energy of the projectile ions. The studies are fundamentally 
interesting from at least two main perspectives. One question involves the 
charge stability of the molecule. How many electrons can one remove from 
Ceo and still retain the molecular structure? High impact parameter collisions 
with highly charged ions provide a means of gently removing electrons from 
the molecule without transferring significant amounts of energy to 
vibrational degrees of freedom. This allows the production of high fullerene 
charge states and their stability can be studied using time-of-flight mass 
spectrometry or other techniques. 

The second questions concerns what happens to the projectile ion. 
Collisions between highly charged ions and surfaces have been extensively 
studied for a number of years [1]. So-called "hollow atoms" are produced 
during such collisions. Hollow atoms are short-lived, multiply excited atoms 
which carry a large proportion of their electrons in high n levels while inner 
shells remain transiently empty. This unusual population inversion arises for 
a time on the order of a few lO's of femtoseconds during the interaction of a 
slow highly charged ion with a solid surface. Despite the limited lifetime it 
is possible to obtain information concerning the formation and decay of the 
hollow atoms by studying electron and X-ray emission as well as the final 
charge and energy of the scattered projectiles [1]. Fig. 6.1 shows a schematic 
picture of the formation and decay of hollow atoms at surfaces. In many 
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ways a fullerene behaves like a small piece of bulk matter. Very similar 
phenomena can occur in collisions between slow highly charged ions and 
fullerene molecules. The big advantage of these systems is that they are gas 
phase collisions and the projectile is not immediately destroyed after 
formation. This makes it possible to probe the projectile properties after the 
collisions and to study the detailed dynamics with the help of coincidence 
experiments. 




Figure 6.1. Schematic picture of the processes that occur when a slow highly charged ion 

collides with a surface. 



1. HIGH CHARGE STATES OE C^o 

Multiply charged fullerenes can decay via fission processes due to 
Coulomb instability in which they emit charged fragments. The Rayleigh 
limit against fission for charged objects is defined via the equality between 
the dissociative Coulomb force and the attractive cohesive force [2]. It is 
difficult to probe this limit for finite highly excited systems such as clusters 
and large molecules due to their internal excitations that also induce 
instability. The limit was demonstrated experimentally for the first time only 
very recently in studies of fission of charged microdroplets [3]. There have 
been a number of theoretical investigations of the charge stability of Ceo- The 
predicted theoretical limit for Ceo'*^ ranges from 17 = 10 [4] to = 13 [5] and 
q=l6 [6]. 

Highly charged ion collisions have been used to produce highly charged 
C6o‘*^ ions from neutral Ceo- Since electron capture occurs at large impact 
parameters where only small amounts of energy are transferred to the 
internal degrees of freedom of the molecule, this is thought to be one of the 
most promising methods to probe the charge stability limit. Charge states up 
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to q = 9 [7] and <7 = 10 [8] have been observed in collisions with and 
Xe^^’^, respectively. The ionisation and fragmentation patterns obtained in 
highly charged ion collisions are very similar to those obtained in 
femtosecond laser ionisation experiments or electron impact experiments. 
This is not too surprising since the interaction times are similar in all three 
cases and the interaction predominantly induces electronic excitation. Fig. 
6.2 shows a comparison between mass spectra produced by collisions with 
Ar*’^ and by 800 nm, 25 fs laser irradiation. By making use of coincidence 
collision experiments, it is possible to select out events that minimise the 
transfer of energy to internal degrees of freedom of the fullerene and thus 
reduce the amount of fragmentation observed in the mass spectra. An 
example of this is shown in Fig. 6.3 for collisions between Ar*"^ and Ceo at a 
collision energy of 16 keV [9]. Here, the spectra have been recorded in 
coincidence with different outgoing projectile charge states Ar^*’**"^. In Fig. 
6.4, we show the evidence for production of in collisions with Xe^^"^ at 
a collision energy of 280 keV [8]. 




Figure 6.2. Comparison between ionisation /fragmentation spectra for highly charged ion 
collisions (upper spectrum, collisions with Ar*"^ at 26 keV [10]) and fs laser excitation (lower 
for 25fs 800 nm pulses with an intensity of 10*^ W/cm^). The large series of peaks in the laser 
spectrum just below is due to Xe"^, used to calibrate the laser intensity. 
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Mass/Charge 

Figure 6.3. Fragmentation spectra measured in coincidence with different outgoing projectile 
charge states (8-.s) with s =1-4,6 for collisions between Ar*’^ and at a laboratory collision 
energy of 16 keV. Adapted from [9]. 
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Mass/Charge 

Figure 6.4. Single stop time-of-flight mass spectrum of highly charged € 60 “*^ produced in 
collisions with Xe^^"^ at a laboratory collision energy of 280 keV. The spectrum provides 
evidence for € 60 **^^, superimposed on a singly charged Cf fragment. Adapted from [8], 

Recently, unambiguous evidence has been obtained for via 

multiphoton absorption of infrared photons (1800 nm) in femtosecond laser 
experiments [11]. This is the highest fullerene charge state observed so far. 

All the experiments reported so far have used an effusive Ceo beam 
source and thus the molecules have a relatively high initial internal 
vibrational energy (see Chap. 2.1.3). It may be possible to observe even 
higher charge states if the fullerene targets are initially in the vibrational 
ground state. Experiments with vibrationally cold Ceo are possible by using a 
cluster aggregation source [12] and collision experiments with such a source 
are expected to become available in the near future [13]. 

Recent measurements of the kinetic energy release in fragmentation 
reactions of highly charged C6o‘*^ have provided information on the stability 
of the multiply charged fullerene ions [14]. These experiments were carried 
out using the coincidence set up shown in Fig. 2.15. The kinetic energy 
release was determined from the flight time and lateral spread of the 
fragment ion signals detected with a position sensitive detector. The 
fragment ions were detected in coincidence with a final projectile charge 
state. The results of the kinetic energy release measurements for the process 
where C6o‘‘^ emits a single Cf' fragment are shown in Fig. 6.5. The 
competing channel where the fullerene emits a neutral C 2 fragment shows a 
significantly lower value for the kinetic energy release [14]. The 
measurements are compared with measurements using different methods 
such as the MIKE (mass-analysed ion kinetic energies, Fig. 2.16) [15,16] 
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Figure 6.5. Kinetic energy release for the fission process where a multiply charged C6o‘’’^ ion 
emits a C 2 * fragment. Adapted from [14], 

and time-of-flight techniques [17,18]. In the latter experiment, the 
coincidences between two fragment ions from the same event are 
determined. A simple model was used to interpret and discuss the kinetic 
energy release measurement results [14]. The two separating fragments are 
treated as conducting spheres The model was first developed by Naher et al. 
in 1997 [19]. This was later extended to take the effect of electron transfer 
during fragmentation into account [20]. This is the basis of the autocharge 
transfer mechanism [15] in which it is assumed that the fragmentation starts 
out as a separation between a charged heavy fragment (Cs8‘*^) and a light 
neutral fragment (C 2 ) which at a rather large critical distance Rc loses an 
electron to the larger charged fragment. If one considers the potential 
energies in such a model. Fig. 6.6, one can see that it is energetically more 
favourable to emit a neutral C 2 than a € 2 ^ from For q = 3, the two 

potential energy curves cross at a distance of approximately 14 ao which is 
close to but just outside the critical distance for electron transfer, 12.7 ao. 
This implies that although the potential energies of the sphere-sphere system 
are the same for asymmetric fission and evaporation at the crossing point, 
the potential barrier which the electron experiences as it attempts to move 
from the lighter to the heavier sphere is higher than its total energy. Under 
these circumstances the autocharge transfer process would require tunneling 
through a barrier and is not expected to be efficient. This is in agreement 
with the presented experimental results [14]. For q = 4, 5 the fission potential 
energy curve lies below the evaporation curve for all distances and therefore 
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Figure 6.6. The interaction energies for two conducting spheres, used to model the + 

Cf and + C 2 potentials as a function of R for q = 2,3,4 and 5. The potential energies for 
fission approach zero as R approaches infinity. Adapted from [14]. 

there is no crossing and direct Cf emission is expected to dominate. The 
kinetic energy release values estimated from the sphere-sphere model are 
shown on Fig. 6.5 as the dotted line. It is assumed in the model that the 
fission products are in their electronic and vibrational ground states. This is 
not the case in the experiments and the final states will be at least 
vibrationally excited, both due to the initial temperature of the neutral 
fullerene target and also to any energy transferred during the collision, and 
may also be electronically excited. This would yield a smaller difference in 
relation to the maximum of the potential energy barrier and therefore smaller 
values of the kinetic energy release. 

The kinetic energy release in asymmetric fission is related to the fission 
barrier, via the expression [14] 

^E,,, + D^ + IP{C,)-IP{C,,) (6.1.1) 

where D‘‘ is the dissociation energy for C 2 emission from This is 

plotted as a funcdon the initial charge of Ceo in Fig. 6.7, using the presently 
accepted value for the C 2 dissociation energy (10 eV, with the additional 
assumption that this is q independent) (Chap. 3.2.1) and ionisation potentials 
estimated in [21]. The sphere-sphere model gives a stability limit for the 
charge on the fullerene cage of qr = 18. The extrapolation of the experimental 
data (which unlike the model does not involve ions with zero internal 
energy) gives a critical charge of 11. The theoretical predictions of the 
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stability limit are shown by arrows and the thick arrow shows the highest 
charge state observed in fs laser experiments. It will be very interesting to 
compare the present set of experimental results with those that will hopefully 
be obtained in the future using initially vibrationally cold fullerenes. 




Charge of fragmenting C0O , q 

Figure 6.7. Semi-empirical and model fission barriers for ^ + C^. The curves 

show a fit to the experimental data and the results of the sphere-sphere model. The thin arrows 
indicate theoretical predictions of the stability limit [4-6]. The thick arrow indicates the q=\l 
charge state observed in laser experiments [11]. Adapted from [14]. 

The branching ratios between singly-charged light fragments and singly and 
doubly charged light fragments emitted in asymmetric fission reactions have 
been studied by Martin et al. in a similar coincidence setup involving 
collisions between Xe^^"^ and Ceo [21]. The branching ratios for the number 
of carbon atoms in a singly-charged light fragment are plotted in Fig. 6.8 for 
different initial charges, q on the parent molecule. The light fragment 
obviously shift to smaller numbers of carbon atoms as the parent molecule 
charge state increases. The experimental results are well reproduced with a 
simple statistical model showing that the asymmetrical fission of for q 
up to 9 can be described by a thermally activated evaporative process. The 
branching ratios of the total singly- and doubly-charged light fragment 
intensities are plotted in Fig. 6.9 as a function of charge on the parent 
molecule. The total branching ratio of these channels increases from less 
than 2% to 25% as the charge on the parent molecule increases from^^' = 5 to 
q = 9. This is because the fission barrier for the emission of the doubly 
charged light fragment decreases faster than that of the singly-charged 
emission channel as the parent molecule charge increases [21]. 




6. Charge Transfer In Collisions With Highly Charged Ions 



119 



U.B 


r-- 

{ 


OJ 


[..• 







r. A I 

' A / \ 1 





0.? 

rw 

0^ t / \ 

If \ / A 
0 I )sC. 

0.8 






q=5 i 

^ ^ >■ p. — 



[ ! 1 r- 



•R 0.4 


: .A 

/ \ 


CD 

^ 0 ^ 


1 , , 


O) 08 . 




1' • 

■P i q=7 

-C 04 ^ ^ \ ^ 

o r/ . 


C : 

CD 0 1 


1 . . i 


^ ^ 1 


1 ' • ■ 1 




Figure 6.8. Comparison of experimental branching ratios (squares) with a statistical model 
(circles) for the fission channels ^ C6o-m*'' + Cn,^ with q ranging from 4 to 9. Adapted 

from [21], 



2. CHARGE TRANSFER MECHANISMS AND 
FINAL PROJECTILE STATES 

Although many electrons can be removed from the Ceo in collisions with 
highly charged ions, not all of the transferred electrons will be stabilised on 
the projectile ion. A collision system that has been applied by a number of 
groups to the study of this effect is Ar*"^ + Ceo [9,22,23]. The reaction 
equation therefore reads 
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^ Ar ^ Ar<*-^>" + + (^ - ^ )e- (6.2.1) 

where q is the number of extracted ("active") electrons and 5 is the number 
of electrons stabilised on the projectile ion. In experiments, only the products 
on the right-hand side of the equation can be observed, along with any 
fragments from the metastable fullerene ion. The first experiment of this 
kind was carried out by Walch et al. [22]. At an ion energy of 80 keV, up to 
eight electrons were seen to be stabilised by the projectile ion. Projectile ions 
with fixed 5 were detected in coincidence with the corresponding fullerene 
ions and fragments. With increasing number 5 of stabilised electrons, 
enhanced fragmentation of the fullerene targets was observed. In 
coincidence with the Ar’"^ (5 = 1), fullerenes up to CtJ"* but no fragments 
were seen (similar to the situation in Fig. 6.3(a)). This indicates that a larger 
number of stabilised electrons is achieved in close collisions, while pure 
electron transfer processes occur at large impact parameters with relatively 
little excitation of vibrational modes in the fullerene. In the latter case, most 
of the active electrons are lost again by the projectile after capture. This was 
also concluded from a coincidence study of Ceo ions and fragments with Ar 
K x-rays from - 1 - Ceo collisions [24]. Most of the captured electrons 

are lost again shortly after capture via Auger transitions. 




q 



Figure 6.9. Total branching ratios of the singly and doubly charged light fragment emission 
channels versus the charge state q of the parent ion. Adapted from [21]. 

Total absolute cross sections <7^ were measured for the formation of 
stabilised projectile ion states [9,22,23] and are shown in the lower part of 
Fig. 6.10. Below 5 = 4 the cross section increases strongly. There is a small 
maximum around 5 = 5,6 and then it drops down towards 5 = 8 





6. Charge Transfer In Collisions With Highly Charged Ions 



121 



(corresponding to projectile neutralisation). From the cross section data in 
the figure in combination with time-of-flight spectra of the fullerene ions and 
fragments as a function of 5 one can obtain absolute cross sections for the 
production of C6o‘*^, (7q- Experimental capture radii are then obtained from 

( 6 - 2 - 2 ) 

with the total cross section 

= = (6.2.3) 

1 

Rq decreases from 23 ao for = 1 to 14 ao for q = 6. The values of the capture 
radii can be understood in terms of an over the barrier model (see Chap. 
3.4.1). The active electron is assumed to be transferred to the projectile at the 
distance where the potential energy of the electron matches its Stark-shifted 




J I I 1 1 1 1 L 

1 2 3 4 5 6 7 8 

number s of stabilised electrons 



Figure 6.10. Top: charge state fraction fg* = cr/Zas of the Ar projectile ions after collision 
with various targets. Bottom: absolute cross section for stabilisation of s electrons on the 
projectile after Ar®’^ + C 60 collisions. Data from three independent measurements are shown 

[9,22,23]. 
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binding energy in the fullerene (ion). Similar approaches in more or less 
sophisticated forms were employed by a number of groups for the analysis 
of electron transfer in ion-fullerene collisions for various collision systems 
[9,25-31]. It is interesting to note that the capture radii and cross sections 
obtained from an over the barrier model in which the positive charges on the 
Ceo are assumed to be point-like and to move freely on the fullerene surface 
are not too different from the results obtained by treating the fullerene as a 
conducting sphere. 




AE [eV] 



Figure 6.11. Energy gain and loss of projectiles having stabilised 1 to 6 electrons after 
collisions with at a laboratory collision energy of 16 keV. In the case of Ar^"^ and Ar®"^ 
ions the full and dashed lines correspond to ^ = 1-4 active electrons. Adapted from [32]. 

The relation between Ar*"^ -i- Ceo collisions and collisions of highly charged 
ions with atoms and surfaces as well as penetration of highly charged ions 
through carbon foils is shown in the upper part of Fig. 6.10. Here, the 





6. Charge Transfer In Collisions With Highly Charged Ions 



123 



projectile charge fraction // is plotted as a function of the number 

of stabilised electrons, 5 [9]. For low 5 the behaviour of collisions with Ceo is 
very similar to those with Ar atoms. For 5 > 5, ff increases, similar to 
observations in surface scattering and passage of the projectile through a 
carbon foil. The hump around 5 = 6 is present in the Ceo as well as in the 
carbon foil data. It is interpreted in terms of equilibration of the charge state 
on the projectile in small impact parameter collisions, which are frequent in 
the carbon foil experiments. In the case of Ceo it is assumed that trajectories 
through the fullerene cage have a similar effect and are therefore responsible 
for the hump. 

Measurements of the projectile kinetic energy gain in the case of single 
electron capture allows the assignment of the atomic ion state into which the 
electron is captured [33]. Since no projectiles with kinetic energy loss or 
negligible energy gain are observed, it is concluded that no electronic 
excitations in the Ceo target are involved in the single electron transfer events 
(a reason for the fragmentation free spectra. Fig. 6.3). The situation changes 
when a larger number of electrons is captured from Ceo- For s = 6, the 
projectile starts to lose kinetic energy as a consequence of strong electronic 
excitations in Ceo which are produced in small impact parameter collisions. 
Fig. 6.11. 
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Chapter 7 

ELECTRON COLLISIONS 



Electron impact ionisation of fullerenes is one of the most common ways 
of producing a beam of fullerenes for a variety of collision experiments but 
is also interesting in its own right. The ionisation/fragmentation mass spectra 
that are produced on electron impact ionisation show many similarities with 
those produced by highly charged ions or by interaction with ultrashort laser 
pulses. The characteristic feature of the mass spectra is the presence of 
multiply charged fullerenes [1-4]. The extraordinary stability of the fullerene 
structure is manifested in the relative ease with which non-fragmented 
multiply-charged molecules can be produced and detected. The detection of 
electron-impact produced C6o‘*^ with q = 5 [3], 6 [5] and 7 [6] was reported 
in rapid succession in 1994. Appearance energies for the different charge 
states have been reported [7] and are in good agreement with photoionisation 
values [8] and calculations [9], see also Chap. 3.2.1. The fragmentation 
behaviour of fullerenes following electron impact excitation is very similar 
to that of atom or ion impact excitation and can basically be described in 
terms of the statistical models and phenomenological fission models we have 
discussed in previous chapters. We will not repeat this here but instead will 
focus on the electron impact cross section data for both 
ionisation/fragmentation and electron attachment. 



1. ELECTRON IMPACT IONISATION CROSS 
SECTIONS 

Absolute partial and total cross sections were reported by a number of 
groups for electron impact ionisation of Ceo [10-13]. The original cross 
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sections reported by Diinser et al. [10] were later corrected to account for the 
occurrence of metastable fragmentation in the field free region of the mass 
spectrometer [14]. Some of the original and corrected data from the Mark 
group are shown in Fig. 7.1. Fig. 7.2 compares the measured (corrected) 
cross sections for single ionisation of neutral Ceo measured by the different 
groups along with two theoretical calculations. Curve 1 is a calculation using 
the semi-empirical Deutch-Mark formalism [15]. The second theoretical 
curve (curve 6) is based on the same approach but uses a modified additivity 
concept [14]. The calculation uses a Mulliken population analysis of the 
fullerene orbitals in which the outermost 120 Ceo orbitals are treated as linear 
combination of the atomic C(2s) and C(2p) orbitals. The calculated curves 
are significantly larger than the experimental cross sections and also are not 
able to reproduce the shape of the experimental curves More sophisticated 
theoretical models are obviously needed to satisfactorily explain the 
experimental measurements. 




Electron energy (eV) 

Figure 7.1. Partial electron impact ionisation cross sections as a function of electron impact 
energy. Circles from [11]. Squares from [14]. From [14]. 
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Electron energy (eV) 

Figure 7.2. Cgo total single ionisation cross sections. 1: calculated, from [15]. 2: experimental 
from [12]. 3: [13]. 4: experimental from [14]. 5: experimental from [11] corrected by [14]. 6: 
theoretical from [14]. From [14]. 

Cross sections have also been obtained for electron impact induced 
fragmentation and ionisation of fullerene ions [16-19]. the apparatus used for 
these studies is shown in Fig. 2.17. Absolute cross sections for electron 
impact induced C 2 fragmentation of C6o‘‘^ ions (q = 1,2,3) are shown in Fig. 
7.3 [17]. The shape of the cross sections is very similar to those for collision 
induced dissociation of positively charged fullerenes in collisions with 
atomic targets (Fig. 5.5), however the time window in the experiment is 
quite different and the interpretation of the significant peak that appears with 
its maximum around 30 eV is also quite different. Here, the high energy part 
of the cross sections was fitted with the Lotz formula [20], normally used to 
calculate electron impact ionisation cross sections. This scales linearly with 
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the charge state of the fullerene ions (Fig. 7.3). The dominant peak is 
interpreted as being due to the excitation of the giant plasmon resonance. 
Evidence for this was also found in gas phase electron diffraction 
experiments [21]. 




Electron energy [eV] 



Figure 7.3. Absolute cross sections for electron impact induced C 2 fragmentation of 
ions. Circles: q= I ', squares: ^ = 2; triangles: q = 3. The solid lines are fits of the Lotz formula 
through the data points at energies higher than 120 eV. From [17]. 

A considerable effort has been made to study absolute cross sections for 
electron impact ionisation and fragmentation of negatively charged fullerene 
targets [16,18,19]. Cross sections for the double, triple and quadruple 
ionisation of Cm, C70' and CsT are shown in Fig. 7.4 [16]. There is very good 
agreement between the experimentally determined ratios of the cross 
sections for the double ionisation and the ratios of the geometrical cross 
sections. The triple ionisation data fit well with the squared ratios of the 
geometric cross sections and the ratios for quadruple ionisation match the 
cube of the ratios of the geometric cross sections. This can be described in 
the form of a scaling law [16] 





Cross section (lO'^o m*) Cross section (10 » m*) Cross section (10 » m^) 



7. Electron Collisions 



129 




Figure 7.4. Absolute cross sections for electron impact double ionisation (upper), triple 
ionisation (middle) and quadruple ionisation (lower). From [16] 
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a^{n,q) = a^{60,q) 
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(7.1.1) 



where cr,, and ageom are the electron impact ionisation and the geometrical 
cross section, respectively. The dependence is shown in Fig. 7.5. Instead of 
the number of electrons that are removed, (<?+l), the charge state of the 
product ion q corresponds with the power multiplying the geometrical ratios. 
This implies that the ejection of one of the electrons does not contribute to 
the magnitude of the cross section in each of the observed ionisation 
processes. It is likely that the detachment of the extra electron is responsible 
for this behaviour. In the case of pure ionisation of anions to positively 
charged ions, the energetic thresholds are about 10 eV higher than expected 
on purely energetic grounds [19]. This indicates that during the approach of 
the projectile electron, the attached electron will be pushed to the opposite 
side of the fullerene by Coulomb repulsion and will therefore not interact 
with the projectile. In addition, its binding energy (2.65 eV) will be reduced. 
Close to the threshold energy of an ionisation cross section all electrons 




Figure 7.5. ratios of the measured cross sections for double (q = 1), triple (q = 2) and 
quadruple (q = 3) ionisation and the ratios of the geometric cross sections raised to the power 
of the respective final charge states q. From [16]. 
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leaving the ion are slow. There will therefore be a relaxation of the formerly 
attached electron rather than its removal. However, at sufficiently large 
energies above the threshold, one of the emitted electrons which is normally 
ascribed to the former projectile electron will have enough energy to push 
the weakly bound electron completely away from the fullerene. On its way 
through the fullerene a fast projectile electron will most likely interact with 
the 7i-orbitals. As the number of electrons is so large, the removal of only a 
few electrons can be assumed to have only a little effect on the size of the 
orbitals. Therefore the cross section for the removal of q electrons from the 
7i-orbitals can be expected to scale with the size of the 7i-orbital cloud to the 
power q and thus also with the size of the fullerene to the same power [16]. 
The scenario is illustrated schematically in Fig. 7.6. 






Figure 7.6. Model for the multiple ionisation of negatively charged fullerene ions hy electron 

impact. From [18]. 
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2. ELECTRON ATTACHMENT 

Fullerenes have high electron attachment cross sections for low-energy 
electrons with a long plateau in the cross section for attachment to Ceo 
existing until energies of around 10 eV. There has been a lot of discussion in 
the literature about the electron attachment cross sections. The results of 
early experiments [22-27] indicated that there was a barrier for attachment at 
zero incident energy and thus no s-wave capture. This was rationalised 
theoretically by Tosatti and Manini [28] (admittely with some rather crude 
approximations) and the apparent activation energy for attachment was 
attributed to a centrifugal barrier to electron capture in the p-wave channel of 
about 0.24 eV, in excellent agreement with experiment [23]. However, this 
has since been shown to be misleading. Other experiments [29,30] gave clear 
indications that s-wave capture can occur. This is in agreement with Rydberg 
electron transfer experiments [22,31,32] that also showed s-wave capture 
with either alkali or argon Rydberg atoms up to very high n. More recent 
experiments have provided additional convincing evidence for a strong 
enhancement of the electron attachment cross section for E — >0 [33]. In this 
work, a beam depletion technique was used to determine the absolute 
integral attachment cross section. The results are shown in Fig. 7.7. The 
points show a strong rise for decreasing electron impact energy. The 
experimental data is compared with the Langevin cross section, calculated 
using the experimental value of the polarisability and convoluted with the 
experimental electron energy spread. There is a correlation between the 
experimental data and the trend in the Langevin cross section although the 
experimental data are consistently lower in intensity. The capture of free 
electrons by fullerenes can be considered as a two-stage process. The first 
stage is governed by the polarisation interaction and the Langevin curve is a 
measure of how efficiently the electrons are attracted by the particle's 
polarisation field. The full attachment cross section, however, also depends 
on how efficiently these electrons "stick" to the fullerenes [33]. One can 
define a sticking probability via 

^ attahment ^ Langevin^ ‘ (7.2.1) 

P, the ratio of the experimental cross section to the Langevin value is clearly 
considerably less than 1 over the experimental energy range and shows some 
interesting structure. The results have consequences for modeling the reverse 
process of thermionic electron emission from fullerene anions. 
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Electron Energy [eV] 

Figure 7. 7. Experimentally determined integral cross sections for free electron attachment to 
Ceo- Dashed line: Langevin polarisation capture cross section. Adapted from [33]. 
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Chapter 8 

ENDOHEDRAL FULLERENE FORMATION 



Endohedral fullerenes (from the Greek "endo" meaning within) are 
fascinating molecules where one or more atoms are "captured" inside the 
carbon cage. An entire volume of the series "Developments in Fullerene 
Science" has been devoted to the production and characterisation of these 
novel materials [1]. The name "endohedral fullerene" was coined by Helmut 
Schwarz, the first person to show that an atom could be captured inside an 
intact carbon cage during a gas phase collision [2]. The idea of filling the 
fullerene cage is almost as old as fullerene research itself. Shortly after the 
discovery of buckminsterfullerene [3], evidence was found by the Rice 
group that lanthanum atoms can be trapped inside these empty carbon cages 
in experiments where a LaCla-doped graphite target was laser vapourised 
[4]. In the corresponding mass spectra, mass peaks according to LaC„, 44 < n 
<76, were observed. The fact that in these compounds only even numbers of 
carbon atoms occurred was interpreted as a consequence of the presence of 
fullerene-like carbon cages surrounding the La atom. This hypothesis was 
later confirmed in "shrink-wrapping" experiments where the carbon cage 
was photofragmented down to the smallest size still capable of containing 
the metal atom [5]. Atom-containing fullerenes were given the appealing 
designation M@C„. A number of different endohedrally metal-doped 
fullerenes could be synthesised by laser vapourisation [6, 7]. In addition to 
the shrink-wrapping technique, other experiments concerning the reactivity 
of the metal-fullerene compounds in the gas phase provided convincing 
evidence for the endohedral nature. Reactions of endohedral compounds 
with molecules such as O2, NH3, N2O or H 2 O do not occur in the gas phase 
whereas externally ("exohedrally") bound metal-fullerene species are readily 
oxidised [5, 8]. 
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Other, more-efficient methods for the production of endohedral 
fullerenes were also developed. Using metal or metal oxide doped graphite 
rods it is possible to produce macroscopic amounts of endohedral fullerenes 
using the Kratschmer-Huffman arc discharge method for fullerene 
production [1, 9]. Although many different endohedral fullerenes can be 
observed in the mass spectra of the raw soot produced by this method, only 
the most stable compounds can be isolated and purified. In particular, it has 
proved very difficult to isolate M@C6o using this method although some 
progress has been made in recent years [10, 1 1]. 

In this Chapter we will focus on the collisional formation of endohedral 
fullerenes, both in the gas phase and on substrates. For details of purification 
and characterisation of the more conventional and readily available arc- 
discharge produced materials, the reader is referred to the literature [1,9]. 



1. GAS PHASE COLLISIONAL CAPTURE OE 
ATOMS AND IONS 



1.1 Capture of Rare Gas Atoms 

The collisional production of endohedral fullerenes in single gas-phase 
collisions was first reported by Weiske et al. [12] shortly after fullerenes 
became generally available in the early nineties. These experiments involved 
high-energy (8 keV laboratory collision energy, corresponding to a centre- 
of-mass energy of 44 eV) collisions between Ceo^ ions and neutral rare gas 
(RG) atoms. The apparatus used for an extensive series of studies is shown 
in Fig. 2.9. The energetic fullerene ions collided with helium atoms in the 
central field-free region of the four-sector mass spectrometer. Apart from the 
usual appearance of fragment ions from C„^, the collision-induced 
fragmentation pattern of Ceo^ also showed a second distribution of mass 
peaks shifted upwards by four mass units compared to the pure carbon 
fragments. The observation was interpreted as being evidence for the 
formation of an endohedral He @€50^ complex since an exohedrally bound 
HeCeo^ ion should lose the loosely bound He atom before it undergoes loss 
of tightly bound C2 molecules. Calculations indicated that the collision 
energy was sufficient to force a helium atom through a six-membered ring of 
the carbon cage [13, 14]. The experimental result was reproduced shortly 
afterwards by two different groups [15, 16] using three- and four- sector mass 
spectrometers. Figure 8.1 displays two mass spectra that were obtained by 
Ross and Callahan for collisions between Ceo^ and He at 44 eV (centre-of- 
mass) [16]. The upper spectrum displays the product distributions from these 
collisions. Apart from the primary ions, Ceo^, endohedral product ions. 
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mass-to-charge ratio 



Figure 8.1. Three-sector mass spectra of collisions between and He at 44 eV (centre-of- 
mass). Adapted from [16]. Top: product distribution after unimolecular loss of C 2 in the field- 
free region of the mass spectrometer. Bottom: product distribution after collisions with Xe 
atoms at 200 eV. The fact that the He complexes survive the second collisions provides 
evidence for their endohedral character. 



He@C6o^ (mass 724u) and fragment ions and He@C 5 g'^ are visible. 
Evidence for the endohedral character of the products can be seen in the 
lower part of Fig. 8.1 that shows the result of an experiment where the 
primary collision products were collided with Xe atoms at 200 eV (centre- 
of-mass). The He-carbon complexes evidently survive these collisions, 
which can only be explained if capture of He inside the fullerene cage is 
assumed. 
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Another clear indication of the endohedral nature was obtained in time- 
of-flight mass spectrometry experiments. In these experiments it was much 
more straightforward to change the collision energy than in the complex 
commercial multi-sector machines. Measurement of the collision energy 
dependence of the formation of He@Ceo^ yielded an energetic threshold of 
approximately 6 eV for He to penetrate the fullerene cage [17], in good 
agreement with predictions for ring penetration [13, 14, 18]. 

In subsequent work by the Schwartz group, it was shown that He can also 
be inserted into doubly and triply charged fullerene ions [19] and that 
sequential insertion of ^He and “^He can occur in separate collisions [20]. One 
of the nicest experiments of this time was the illustration of the stability of 
the endohedral complex by carrying out neutralisation-reionisation 
collisions. Fig. 8.2 [21]. Interestingly, no direct capture products He@C 5 g/ 68 ^ 
but only their post-capture fragments He@Cn'^, n < 58,68 were observed in 
collisions between and He at 44 eV [22], reflecting the lower stability 
of Cs 8/68’" ions relative to Ceo/To""- 

A deeper understanding of the capture process has been obtained from 
measurements of the capture cross sections as a function of the collision 
energy. Two different methods have been used. Time-of-flight mass 
spectrometry was used by the author's group in order to investigate the 
capture of He and Ne atoms by positively and negatively charged fullerene 
ions [17, 23-25]. The Anderson group used a triple-sector guided ion beam 
mass spectrometer to study capture of Ne''' ions in neutral Ceo [26, 27]. In 
contrast to the majority of studies discussed so far in this chapter, these 
experiments allowed the clear determination of the energetic thresholds for 
capture and of absolute capture cross sections. It was shown that the internal 
energy of the fullerene ions has a considerable influence on the capture 
threshold. Different ion sources were used to study this effect. Using Ceo^ 
ions produced by laser desorption from a fullerene fdm yielded thresholds of 
6±2 eV and 9±1 eV for the capture of He and Ne, respectively [17]. These 
values increase by approximately 1 .5 eV if negatively charged Ceo projectile 
ions, produced by electron attachment in an electron impact ionisation 
source are used [24], Fig. 8.3. Negative fullerene ions produced in this way 
have lower internal energies than do hot positive ions from a laser desorption 
or electron impact ionisation source because electron attachment occurs 
efficiently only for electron energies below 10 eV [28]. This internal energy 
effect on the capture threshold has also been confirmed by classical MD 
simulations [29]. The simulations implied that the internal energy of the 
laser desorbed positively charged ions should be approximately 30 eV higher 
than that for the negative ions produced by electron attachment in order to 
explain the shift in capture threshold observed in the experiments. The 
excitation energy of Ceo can be estimated to be 7.5 eV (thermal energy from 
the oven plus the electron affinity) which would then imply an internal 
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Figure 8.2. Neurtralisation - reionisation experiment [30] A: Mass -resolved € 60 "^ ions. The 
multiple peak structure reflects the 1.1% natural abundance of B: mass-selected 
endohedral He@C6o products of the first collision. C: Selection of "'He@'^C6o. D: Ion signal 
after neutralisation of He@C6o followed by collisional reionisation. Adapted from [30]. 



energy of approximately 35 eV for Ceo^- This can be compared with the 
estimate of 30-40 eV based on the amount of metastable fragmentation 
observed from the Ceo^ parent ions from the source in these experiments and 
the fragmentation rate constant, discussed in Chap. 3. 

The decreasing capture threshold with increasing internal energy reflects 
the higher flexibility of an excited fullerene molecule. The effective size of 
the hexagons in the fullerene cage increases and therefore leads to a decrease 
of the penetration barrier [29]. The measured He capture threshold decreases 
even down to as low as 3 eV if Ceo projectiles from a laser desorption source 
with high fluence (and thus high internal energy) are used in the collision 
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Figure 8.3. Cross sections for the production of endohedral He@Cjo- a) negatively charged 
Ceo (b) positively charged Ceo^- The different threshold can be attributed to the initial internal 
energy of the fullerene ions. From [24]. 



experiment [25]. This value is extremely low, however, a 3 eV escape 
threshold has also been estimated from high pressure experiments (discussed 
in section 8.1.2, below). A possible explanation for these observations was 
provided by Murry and Scuseria [31]. In their semi-empirical and ah initio 
calculations they observed the formation of thermodynamically stable 
windows on the potential energy surface of the triplet manifold of Ceo 
corresponding to nine- and ten-membered rings in the fullerene cage. Fig. 
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8.4. Such "windows" have also been seen in MD simulations of Ceo with an 
internal energy of approximately 40 eV, due to the high propensity for defect 
formation in the cage at such high excitation energies [32]. This, combined 
with calculations that show that annealing processes are much less efficient 
with anions than with cations [33], could explain the reduced capture 
threshold for fullerene anions with very high initial internal excitation 
energies in the collision experiments [25]. It also provides a plausible 
explanation for the rather low value determined for capture of Ne using 
laser-desorbed fullerenes [17]. 




Figure 8.4. Window on the potential energy surface of the triplet manifold of Ceo suggested as 
an explanation for the high pressure encapsulation of rare gas atoms in Ceo- Adapted from 

[31]. 

The collision energy dependence of the reaction products in collisions 
between Ne"^ and Ceo was reported by Christian et al. [26]. Fig. 8.5 shows the 
relative cross sections for collision-induced fragmentation (top) and capture 
followed by fragmentation (bottom) as a function of collision energy. The 
formation of these ions is accompanied by charge transfer from the neutral 
Ceo target to the Ne"^ ions (note that this process is exoergic by 13.95 eV). 
The charge transfer cross section, represented by the Ceo^ ion signal in the 
upper part of the figure, is approximately a factor of 100 times larger than 
the cross section for the reactive collisions. The shapes of the and 
Ne@C„'^ ion signals are very similar. With increasing collision energy, the 
intensity goes from a threshold region through a maximum after which it 
decreases. The positions of threshold and maxima shift to higher energies 
with decreasing fragment size, n. The threshold shift, however, is not present 
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Figure 8.5. Relative cross sections for charge transfer and collision induced fragmentation 
(top) and capture followed by fragmentation (bottom) in Ne"^ + Cgo collisions. Adapted from 

[26], 

for the first few members of the Ne-containing fragments with n = 60, 58 
and 56 for which the value is close to 25 eV. The NelgiCeo^ signal (not 
shown on plot) mostly originates from capture of Ne"^ by C 70 and subsequent 
fragmentation (15% of the target gas consisted of C70) since most of the 
primary Ne@C 7 o'^ fragment on their way to the detector. Note that the flight 
times in these experiments are roughly three orders of magnitude longer than 
in the time-of-flight or three- or four-sector mass spectrometer experiments 
discussed above. The threshold energy of 25 eV for the detection of 
Ne@C 6 o^ is in rough agreement (although somewhat high) with capture 
thresholds reported from other experiments and theory (see below), however, 
the electron transfer which occurs before the actual atom capture in this 
study may limit direct comparability. The similarity of the fragment patterns 
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of pure carbon and Ne-containing ions demonstrates that evaporative 
cooling, essentially via loss of C 2 molecules, is basically the same for 
fullerenes and rare gas endohedrally doped fullerenes. This is also in 
agreement with measurements of kinetic energy releases in metastable 
fragmentation of endohedral Ne@C6o^ [34]. 

There are actually two mechanisms for capturing He in collisions with 
fullerene ions. Figure 8.6 depicts relative cross sections for capture of He 
and Ne atoms by Ceo^ as a function of inverse collision energy [23]. An 
expression describing the reaction cross section in terms of a simple 
absorbing sphere or line-of-centres model [35] is fitted to the experimental 
data 



< 7 „ 







( 8 . 1 . 1 ) 



where E,hres is the energetic threshold for capture. The total energy during the 
collision can be expressed as the sum of kinetic energy T and potential 
energy V(R), E = T + V(R). The kinetic energy is the sum of radial and 
centrifugal contributions 
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where R is the intermolecular distance, / the rotational quantum number, // 
the reduced mass and b the impact parameter. It follows that 
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(8.1.3) 



The capture cross section can be obtained from the assumption that the 
intermolecular potential V(R) with the barrier Ethres at the contact distance R 12 
of the collison partners, together with the centrifugal energy, has to be 
overcome by the radial kinetic energy 



E-E., 



EE 



> 0 . 



(8.1.4) 



The capture cross section (eq. 8.1) follows from the maximum impact 
parameter that satisfies this inequality. The same simple model is applied to 
fullerene-fullerene collisions in Chap. 9. 
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Two separate linear regimes are seen for He capture in Fig. 8.6, 
corresponding to two different capture mechanisms. The low energy capture 
for collisions with He is characterised by a penetration barrier of E,hres= 
6±0.5 eV, in excellent agreement with earlier studies [17]. The slowly rising 
cross section as the collision energy increases indicates a very low capture 




Figure 8.6. Capture cross sections for collisions between and He (left) or Ne (right) 
atoms [23]. Squares: experimental capture cross sections. Circles: results of MD simulations. 
Lines are fits of the line-of-centres model (eq. 8.1.1) to the data. 



probability with Rj 2 «0.5 A. This indicates that only central collisions 
through a six-membered ring contribute to He capture at these energies. The 
threshold energy determined for the process that onsets at higher collision 
energies is E,hre.i = 17±0.5 eV. The steeper slope indicates a larger probability 
for capture via this second process. 

In order to understand the mechanism for the high energy capture, we 
consider additional information provided from measurements of inelastic 
scattering in Ceo + He collisions. The upper part of Fig. 8.7 shows the 
experimental kinetic energy distribution of product ions from + He 
collisions at a centre-of-mass collision energy of 27.6 eV [18]. Similar data 
were reported in [36] and [23]. The sharp feature in the broad scan 
corresponds to the endohedral collision product as is evident from the 
narrow scan (below) which was taken with an additional mass filter set to 
724 u before ion detection. The lower part of Fig. 8.7 shows the result of a 
corresponding classical MD simulation that agrees very well with the 
experimental result [23]. The low energy foot of the kinetic energy 
distribution does not correspond to fully elastic collisions with a kinetic 
energy loss in the laboratory system of approximately 4Ecm as would be 
expected for a very light target like He (Chap. 3.1). The energy loss in the 
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Figure 8.7. Top: kinetic energy scans of products from Cso* + He collisions at 4.95 keV (lab). 
Adapted from [18]. In the broad scan both and He@C6o'^ were transmitted to the detector 
while the narrow scan was obtained with an additional mass filter set to 724 u, corresponding 
to the endohedral helium complex. Bottom: kinetic energy distributions of projectiles 
after collisions with He at 5 keV (lab) calculated with classical MD simulations [23]. The 
white bar indicates the position of the kinetic energy of the endohedral complex He@C6o^- 

laboratory frame of reference is smaller and can be understood in terms of 
the spectator model discussed in Chap. 3.1.1. The He atom first collides fully 
elastically with only a part of the fullerene cage e.g. a C atom or a dimer, 
leading to the observed laboratory kinetic energy loss. Subsequently, this 
sub-unit on the fullerene cage that has been scattered in the collision, 
collides with the remaining fullerene cage, finally leading to a statistical 
distribution of vibrational excitation in the fullerene. The effective mass m^jf 
that the atom "sees" on the cage depends on the collision energy and can be 
calculated from the maximum of the kinetic energy loss spectrum [23]. This 
has been plotted in Fig. 3.4 and it can be seen that decreases rapidly for 
collision energies between 10 and 17 eV from 720 u to 24 u and remains 
constant close to the mass of a Cj unit in the energy range where the most 
effective He capture mechanism is active. This indicates that for collision 
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energies above 17 eV, direct knock-on collisions with C 2 sub-units on the 
cage are efficient. During such events, carbon-carbon bond breaking can 
occur, temporarily opening a window in the fullerene network and therefore 
facilitating atom capture. This picture is in qualitative agreement with 
classical MD simulations that predict a similar decrease of m^jf, however at 
lower collision energies [23]. Other classical MD simulations [15, 29, 37] 
have explicitly shown the direct ring penetration mechanism and also what 
are referred to as "window" and "hole" mechanisms. The "window" 
mechanism refers to trajectories where one C-C bond is broken on impact 
and rapidly re-closes leaving an intact fullerene cage with a captured atom 
(as in the higher energy He-capture mechanism discussed above). The "hole" 
mechanism refers to trajectories where multiple bonds are broken on atom 
impact with a large hole remaining on the surface of the cage for up to a few 
hundred fs. These trajectories have a high probability for the captured atom 
to escape from the cage. In these simulations. He capture is dominated at low 
collision energies by the ring penetration mechanism with a significant 
contribution from the window mechanism for collision energies beyond 
roughly 20 eV, in good qualitative agreement with the experiments (Fig. 
8.6). The experiments do not cover the energy range (> 30 eV) where the 
hole mechanism is predicted to contribute (but see section 8.2.1 below), 
however, earlier experiments noticed a strong decrease in the capture cross 
section for collision energies beyond this value [17, 18]. This is also 
observed for collisions between Cgo and Lf^ ions (see section 8.1.2, below). 

The situation for Ne capture is rather different. Experimentally, there are 
also two mechanisms observable with highly excited laser desorbed fullerene 
ions (Fig. 8.6). The low-energy threshold at about 10 eV is attributed to 
penetration of the atom through defects, already present on the cage surface 
prior to collision. The higher energy process is due to a bond-breaking 
mechanism [23]. In the Anderson experiments. Fig. 8.5, the fullerenes are 
not so highly excited prior to collision and the low-energy capture 
mechanism is not observed. Taking into account the experimental errors and 
the differences in temperature in the two kinds of collision experiments, the 
threshold found by the Anderson group [27] is consistent with the high- 
energy threshold determined by the author's group [23]. Interesting 
behaviour is observed when one considers the scattering dynamics by 
analysing the effective mass seen by Ne, as discussed above for He. In 
contrast to He, the effective mass seen in the Ne experiments does not drop 
to 24 u but reaches a minimum of only about 60-100 u for collision energies 
beyond 30 eV [23]. This is consistent with the larger Ne atom colliding with 
a Ce-ring on the cage, inducing multiple bond-breaking and thus 
predominantly being captured by the hole mechanism as discussed by Cui et 
al. [37] and inferred by them to be the main capture mechanism for Ne. 
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The capture cross section drops with increasing collision energy (Fig. 
8.5). This is due, partly, to fragmentation of the endohedral complex on the 
ps or ms time scale of the experiments. For small captives such as Fie or Li 
(discussed below) the loss of the guest atom may compete with the 
evaporation of C 2 molecules. This is particularly likely at high collision 
energies where there is a greater probability for the production of a 
metastable hole in the fullerene cage on impact. Secondly, at high collision 
energies, direct "shoot-through" trajectories are possible during which an 
atom enters the fullerene but is not trapped and leaves the cage some fs after 
its capture. Such trajectories have been observed in MD simulations [38, 39]. 

Ar-doped endohedral fullerenes have also been reported in collision 
experiments with the subsequent formation of Ar@C6o^ fragments with an 
odd number of carbon atoms [15]. This observation rules out the standard C 2 
loss fragmentation mechanism and indicates a direct collisional induced 
fragmentation rather than a statistical mechanism, although to the author's 
knowledge it has not been confirmed in other studies. Collision-induced 
fragmentation of pre-prepared Ar@C6o^^ leads to the typical even-numbered 
fragment distribution [40]. The larger Kr and Xe noble gas atoms cannot 
form metastable endohedral complexes in collisions under single-collision 
conditions. Flowever insertion of these atoms into Ceo is possible with the 
high-pressure technique of Saunders [41] or by using the nuclear recoil 
method [39, 42]. Both these methods are discussed below. 

1.2 High Pressure Encapsulation of Rare Gas Atoms 

The discovery that Fie atoms could be incorporated into fullerenes during 
the arc discharge production process was made by M. Saunders [43], 
however the yield was extremely low and the discovery was only made 
possible by using highly specialised mass spectrometry methods. 
Subsequently, the Yale group developed a high pressure - high temperature 
method of incorporating rare gas atoms into the fullerene cage [44]. In this 
method the fullerenes are heated to a temperature of 600 °C for a number of 
hours under a high pressure (2500 atm) of the rare gas atom. Analysis of the 
fullerene material after this treatment shows that approximately 0.03-0.1% of 
the fullerene molecules contain the rare gas atom. The presence of Fie inside 
the fullerene cage does not significantly affect its chemical or other 
properties. NMR experiments with endohedral fullerenes containing ^Fle 
have therefore been used to monitor the magnetic field at the centre of the 
cage during chemical reactions [45]. The presence of the larger rare gas 
atoms, Ar and Kr, inside the cage does affect the cage properties slightly. 
Both Ar@C6o [46] and Kr@C6o [47] have been partially purified using 
FIPLC (high pressure liquid chromatography). Although the UV-vis 
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absorption spectrum of Kr@C6o is very similar to that of Ceo, it does show 
some significant differences such as a red shift of peaks in the 'Tig — > 'Aig 
dipole forbidden band and a blue shift of the electronically allowed band at 
400 nm. 



1.3 Gas Phase Collisional Capture of Alkali Ions 

The capture of alkali ions in neutral fullerenes is similar in many respects 
to the capture of rare gas atoms in charged fullerenes, discussed in section 
8.1.1 above. The radii of the two smallest alkali ions (Li"^ and Na"^) are 
comparable to those of He and Ne and we may expect similarities in the 
dynamical behaviour on capture. Capture of alkali metal ions by fullerenes 
has been investigated in detail by the Anderson group [48, 49]. In their 
experiments, metal ion projectiles collide with neutral fullerene target 
molecules. Mass-separated collision products are detected as a function of 
the collision energy. The experiments are similar to the Ne"^ capture studies 
from the same group (discussed above in section 8.1.1) with the important 
difference that in the case of alkali metal ions no charge transfer is involved 
in the collision process due to the lower ionisation potentials of the alkali 
atoms compared to the fullerenes. At the low collision energies of these 
experiments endoergic charge transfer is of no significance. Figure 8.8 
shows the relative cross sections of the different reaction channels in 
collisions between both Lf^ and Na"^ with Ceo (the target gas cell contained 
15% C70). The intensity of the metal -doped parent ions Li@C6o^ and 
Na@C6o^ is plotted in Fig. 8.9 as a function of 1/Ecm, as was done for He 
and Ne collisions, above. The parent ion signals behave in a very similar 
way to the rare gas capture collision products. Both sets of data show 
evidence of two collision mechanisms. The Li product has a low energy 
threshold at approximately 5.5 eV with a more efficient process onsetting at 
approximately 15 eV. As for He collisions, this can be related to the 
threshold for low impact parameter penetration of a six-membered ring and 
for capture via a C 2 bond-breaking mechanism. The Ce penetration barrier 
for Li"^ is in good agreement with calculations [50, 51]. There does not 
appear to be a temperature dependence in the capture probability for Lf^, as 
shown in a detailed series of quasiclassical trajectory calculations [50]. This 
is in contrast to the situation with He [25, 29] and may reflect the smaller 
radius of Li"^ (0.2 A compared with 0.3 A for He). No experimental results 
on the temperature dependence are available for the alkali ion projectiles. 
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Figure 8.8. Relative cross sections for the formation of endohedral metallofullerenes in alkali 
ion - C 60 collisions as a function of collision energy. Adapted from [48]. Left: Li’^; right: Na"^ 

projectile ions. 

The Na@C6o^ product has a similar collision energy dependence as the 
Ne@C6o^- The first threshold lies at approximately 15 eV (corresponding to 
the C 2 bond-breaking also seen for Li"^). Na"^ is too large to penetrate an 
intact six-membered ring in the carbon cage. A second capture mechanism is 
seen to onset for a collision energy of approximately 28 eV. This can be 
related to the "hole" mechanism discussed in section 8.1.1. A more recent set 
of experimental results for Na"^ - 1 - Ceo [52] is also shown on Fig. 8.9(b) for 
comparison. The newer data was obtained using time-of-flight mass 
spectrometry in contrast to the ion guide method used for the earlier 
experiments. The newer measurements have a much shorter and better 
defined time window for detection of the products and also a better collision 
energy resolution [52] . The original reports gave threshold energies of ca. 1 8 
eV [49] and 29.6 eV for Na"^ capture [52]. The discrepancy was attributed to 
the collision energy uncertainty in the earlier experiments. The plot of the 
two sets of data versus 1/Ecm show that in fact the two sets of data do not 
give conflicting threshold energies. Both data sets show a threshold at 
approximately 30 eV. The earlier data also show a lower energy capture 
threshold. The linear extrapolation places this threshold at approximately 15 
eV, in agreement with the Lf^ measurements and the C 2 bond-breaking 
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Figure 8.9. Yield of (a) Li@C6o^ and (b) Na@C6o'^ plotted versus 1/E. Circles: data from [48]. 
Squares: data from [52]. Lines: Fits to the data using eq. (8.1.1). Thresholds for Li capture are 
found to be at 5.5 eV and 15 eV. The Na data shows a high energy threshold for both data sets 
of 29+1 eV and a low energy threshold from the Anderson data of 15+1 eV. 

picture. It is not clear why the low energy mechanism is not observed in the 
later experiments perhaps it is simply a matter of sensitivity. The two data 
sets have been normalised to the same intensity maximum in Fig. 8.9. One 
would expect a larger fragmentation of the product ion in the guided ion 
beam experiment since there is a much longer time before detection (0.5-5 
ms compared with 10 ps for the time-of-flight experiments). Thus the 
maximum in the parent endohedral ion intensity should be decreased relative 
to the low energy tail compared to the time-of-flight data. (For an internal 
energy of 40 eV, a rough estimate, taking account of radiative cooling, yields 
that ca. 90% of the parent endohedral ion should survive to be detected in the 
time-of-flight setup compared to ca. 70-80% in the guided ion beam 
experiment.) Of course, one cannot rule out the effect of the poorer collision 
energy resolution in the early experiments, and the different widths of the 
data in Fig. 8.9 show that this is indeed significant, but it is unlikely to be 
able to explain the full extent of the low energy tail in the guided ion beam 
data. 

The initial internal energy of the collision partners and thus the total 
excitation energy after formation of the endohedral compound is much better 
defined in the alkali ion implantation experiments compared to experiments 
using fullerene ion beams. The fullerene target has an internal energy 
defined by the oven temperature with a corresponding thermal energy 
spread. The ion source has an energy spread of approximately 0.3 eV and the 
amount of energy transferred to internal degrees of freedom is just the 
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centre-of-mass collision energy. Deng and Echt made use of this to estimate 
the binding energy of C 2 in Na@C6o^ by modelling the collision energy 
dependence of the parent ion and the Na@Cs8^ fragment ion intensities [52]. 
By adopting a pre-exponential factor of A = 5 x 10*^ s’^ (eq. 3.2.1) they 
obtained a value of D = 10.17±0.02 eV (only statistical error). This is 
consistent with the value of 10.3 eV for C 2 loss from Ceo^ derived by 
considering the competition between thermionic electron emission and 
fragmentation in the neutral channel [53], however this work fitted a pre- 
exponential factor of 3.4x10^^ s '. If the pre-exponential factor was increased 
for Na@C6o^ this would require a corresponding increase of D in order to 
still fit the experimental results (an increase of A by one order of magnitude 
would increase the derived binding energy by 0.62 eV [52]). A surprisingly 
large enhancement of the C 2 binding energy by ca. 23% for rare gas 
endohedrals Rg@C6o^ (Rg = Ne, Ar, Kr), formed by the high pressure 
method discussed above, compared to was reported by Laskin et al. 
[54]. However, there are many difficulties associated with determining the 
binding energy from kinetic energy release measurements, as was done in 
this work, (see Chap 3.2.1) and one would need an independent confirmation 
to be sure that this is a genuine effect. 



2. ION IMPLANTATION 



2.1 Low Energy Ion Implantation 

The gas phase experiments, discussed above, are related to a low energy 
ion implantation method developed as an alternative way to obtain 
macroscopic amounts of endohedral Ceo compounds [55]. The experimental 
set up used to produce macroscopic amounts of alkali endohedral 
metallofullerenes is shown in Fig. 8.10. Alkali ions are produced from a 
thermal source, accelerated to the desired energy and deposited on a rotating 
metal cylinder with an ion beam cross-sectional area of approximately 2 cm^ 
and a current of typically a few pA. The ion beam energy is controlled to 
within ± 1 eV. Fullerenes are deposited simultaneously from an oven at a 
temperature of approximately 450 °C. The oven temperature, fullerene 
deposition rate and speed of rotation of the metal cylinder are controlled to 
ensure that one monolayer of fullerenes is deposited on each full rotation. 
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Figure 8.10. Schematic diagram of apparatus used to produce macroscopic amounts of 
Li@C6o via low energy ion implantation. 

The deposition and irradiation is continued until films of typically up to a 
few hundred-nm thickness are produced. For ease of extraction of the 
endohedral fullerene material, the fullerenes are deposited on aluminium foil 
wrapped around the cylinder. After a production run the aluminium foil is 
removed from the vacuum apparatus, cut into small pieces and put into CS 2 
to dissolve the fullerenes. 

The films produced by irradiation of fullerenes with alkali ions are 
routinely characterised by laser desorption time-of-flight mass spectrometry 
(LDMS) [56]. An example of an LDMS mass spectrum from a fdm prepared 
at 30 eV Li ion energy and a high ion current corresponding to 
approximately 17 Li"^ deposited per fullerene molecule is shown in Fig. 8.1 1. 
Under such conditions it is possible to produce also Li 2 and Lis endohedral 
species. Flowever, films produced with such high ion currents are rather 
difficult to dissolve and for the routine material production it is more 
efficient to use a lower current corresponding to Lf^iCeo of 1:1. The LDMS 
method gives a lower limit to the efficiency of production of endohedral 
fullerenes. Since the molecules are desorbed with very high internal 
energies, corresponding to temperatures of a few thousand Kelvin, both 
thermal ionisation and fragmentation takes place on desorption [57]. The 
magnitude of the endohedral mass peak compared to empty fullerenes and 
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fragments thus generally depends on the time scale of the ion extraction and 
detection [58]. 




Cftn I Li@C6n 
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Figure 8.11. Time-of-flight mass spectrum from fullerene film implanted with a high dose of 

Ld ions at an energy of 30 eV. 

The yield of endohedral species, as a function of ion implantation energy, 
as determined by measuring the mass peaks in LDMS mass spectra, is 
plotted in Fig. 8.12 for the different alkali metals [56]. The absolute yields 
are not comparable since different ion currents were used for the different 
alkali metals. The endohedral yield has been compared with the sum of 
M@C6 o^ and all endohedral fragment ions measured in the gas phase 
experiments of Wan et al. (Fig. 8.8 and [49]). The distributions coincide very 
nicely. In particular, the energetic thresholds are in good agreement as well 
as the widths of the energetic window within which endohedral fullerene 
formation is possible. The fact that all endohedral fragments of the gas phase 
studies have to be taken into account for the comparison reflects the situation 
that, in the fullerene film, fragmentation is suppressed due to the efficient 
energy transfer from the collisionally excited fullerenes into the film. No 
evidence for exohedral MCeo production during laser desorption was found 
for the small alkalis (Li and Na). Flowever, there is a significant exohedral 
MCeo production for K and Rb [56]. This shows up as an ion energy 
independent background signal on the plots in Fig. 8.12 for energies below 
100 eV (beyond this the fullerenes in the film are destroyed on ion impact 
and the signal decreases). The magnitude of this “exohedral” signal depends 
on the amount of alkalis present in the fdm (i.e. on the ion current used for 
production) and on desorption laser fluence but not on the ion energy. 
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Superimposed on this background signal is the genuine endohedral signal. 
As the size of the alkali increases, the energetic threshold for the endohedral 
production increases and the maximum in the cross section moves to higher 
energy values just as one would expect from the gas phase results. 




Figure 8.12. Comparison of the endohedral yield of material obtained in alkali ion - fullerene 
film collisions (squares) [56] and the sum of all endohedral product ions from the gas phase 
experiments of Anderson and co-workers (circles) [48]. The endohedral fraction of the ion- 
implanted material was determined by laser desorption mass spectrometry. 

Other critical parameters for optimising the efficiency of the production 
process are the ion current (i.e. the ratio of metal ions to fullerene molecules 
deposited) and the film thickness. The alkalis show quite different 
dependencies on the ion current [59]. The differences can be explained in 
terms of the amount of energy transferred to the fullerene film per ion - the 
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larger the alkali ion, the higher is the ion energy required for endohedral 
fullerene production and the higher is the likelihood of destroying the 
fullerenes. However, as discussed briefly above, the films produced with 
high currents of Li and Na ions are not very soluble and there appears to be a 
strong degree of polymerisation occurring for high ion currents which does 
not show up in LDMS. 

Studies of the dependency of the endohedral production efficiency on the 
thickness of the fullerene film deposited on each rotation of the metal 
cylinder showed that the optimum conditions are for the deposition of one 
fullerene monolayer per rotation. The capture rate drops by more than a 
factor of two when going from monolayer to multiple layer deposition per 
rotation. Finally, there is also an influence of the total thickness of the 
deposited fdm. The efficiency decreases practically linearly by almost a 
factor of two as the total film thickness is increased from a few monolayers 
to between 50 and 100 nm. After this value the behaviour reaches a plateau, 
staying constant up to very large film thicknesses [59]. This is a consequence 
of the efficiency of energy transport within the films and will possibly also 
depend on the conductivity of the underlying substrate, although this has not 
been investigated systematically. 

It is clear from the above that a good control of the deposition parameters 
is essential for optimisation of the production method. It is also important to 
avoid neutral metal vapour from reaching the substrate. If large amounts of 
neutral metal are present the dominant extractable reaction products are 
small oligomers (dimers and trimers) of empty fullerenes. This presumably 
happens in a reaction process closely related to the ball-milling production 
method for (€ 60)2 [60]. The ion irradiation serves to heat up the fullerene 
substrate thus facilitating the alkali-metal catalysed oligomerisation 
reactions. Such reactions do not occur to a noticeable degree when only 
metal ions are present (at least at the currents routinely used for production). 

Reports have also been made of a high efficiency production of 
endohedral alkali fullerenes using a plasma method [61]. However, the ion 
energy dependence of the production efficiency is quite different from that 
discussed here (as determined by LDMS), indicating a different formation 
mechanism, and purification and further analysis of the purified material is 
needed to ensure that it is the same product. 

Low energy ion implantation can also be used to produce endohedrally 
doped fullerene derivatives. For example, hydrogenated C 60 H 36 , prepared by 
the transfer hydrogenation reaction [62] was deposited at low temperatures 
and irradiated with 30 eV Lf. The endohedral fragment ions containing even 
numbers of hydrogen atoms are more dominant in LDMS spectra and thus 
more stable than those containing odd numbers [59]. This is the opposite 
tendency from the empty hydrogenated fullerenes [62] and reflects the 
transfer of charge from the Li to the fullerene cage. No attempts have yet 
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been made to extract the endohedral hydrogenated fullerenes using high 
performance liquid chromatography (HPLC). 

Recent spectroscopical studies of the purified Li@C6o materials obtained 
from the low energy ion implantation production method have shown that 
the two fractions that are purified with HPLC are in fact dimers and trimers 
[63]. However, when the purified samples are evaporated at temperatures of 
a few hundred °C the weak inter-cage bonds are broken and the properties of 
monomer molecules can be investigated in the gas phase [64] or they can be 
deposited on substrates [65]. 

A similar method was developed by Weidinger and co-workers to 
implant fullerenes with and P"^ ions [66, 67]. It has proved to be very 
difficult to purify and isolate N@C6o, however, the molecule shows a very 
sharp line in EPR spectra that can be used as a fingerprint. Due to the long 
spin lifetime, it has also been discussed as a possible candidate for quantum 
computation applications [68]. 

The studies described in this section provide a nice example of a situation 
where detailed, fundamental gas phase studies provided the know-how and 
inspiration for the development of a new macroscopic production method for 
interesting new nanomaterials. 

2.2 Nuclear Recoil 



Braun and Rausch [69] showed evidence for radionuclide endohedral 
fullerene formation from neutron irradiation of polycrystalline Ceo as a result 
of Ar impurities trapped in the solid. They proposed that the endohedral 
nucleotide formation was the result of penetration of radioactive "^'Ar atoms 
into the fullerene cages on account of the high nuclear recoil energy 
produced from prompt gamma emission following the nuclear reaction 
‘*°Ar(n,y)"^^Ar. These results were later confirmed by Gadd and co-workers 
[42]. '^’Xe and ^^Kr endohedral compounds were also produced using a 
nuclear recoil method by mixing the fullerene material with or K^®Br 
respectively and irradiating the mixed material with 16 MeV deuterons [39]. 
Similarly the incorporation of tritium atoms (via neutron activation of ^Li) 
[70] and ^Be via neutron activation of ’Li and [71] has been shown 
among many other reactions [72]. 

The recoil energy of a radionuclide resulting from prompt gamma 
emission following neutron activation is given by 



^recoil 



2m, 



( 8 . 2 . 1 ) 
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where Frecou is the recoil energy of a radionuclide with mass Mrecon when a y- 
ray of energy £^is emitted. The maximum energy available for y-emission is 
determined by the Q-value for the neutron activation reaction. This 
parameter is defined as the difference between the nuclidic masses of the 
reactants and those of the products and is easily obtained from nuclear data 
tables [73]. The maximum recoil energies that are obtained for various rare 
gas nucleotides are on the order of 200-500 eV in the centre-of-mass frame 
of reference [74]. These energies are rather high compared to those 
discussed earlier in this chapter but the formation of endohedral fullerenes 
under these conditions is certainly possible via mechanisms of the kind 
discussed for the gas phase reactions. The first volume of this book series 
"Nuclear and Radiation Chemical Approaches to Fullerene Science" [75] 
contains a number of chapters related to the nuclear recoil approach and the 
properties of radioactive endohedral fullerenes. 
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Chapter 9 

FULLERENE-FULLERENE COLLISIONS 



Collisions between fullerenes provide a means to study cluster-cluster 
interactions both theoretically and experimentally. This makes it possible to 
explore the similarities and differences in the dynamical processes in 
collision systems covering many orders of magnitude in size and collision 
energy; from heavy ion collisions in nuclear physics to macroscopic liquid 
droplets or even colliding galaxies. Atomic clusters provide an interesting 
intermediate-size system that can in principle be explored in considerable 
detail, both theoretically and experimentally, with respect to collision partner 
size, binding energy and collision energy. Theoretical investigations of 
cluster-cluster collisions were pioneered by Schmidt and co-workers for the 
case of collisions between two small alkali clusters. Using molecular 
dynamics simulations combined with the density functional formalism 
(Chap. 3.3.2) they were able to show that a number of reaction channels, 
analogous to the dynamics of heavy ion collisions and colliding liquid 
droplets were likely to appear in atomic cluster-cluster collisions, such as the 
fusion of two clusters [1]. Unfortunately, controlled collisions between 
mass-selected atomic clusters under single-collision conditions are extremely 
difficult to study experimentally and so far such studies have only been 
possible for colliding carbon clusters, or fullerenes. Fullerene-fullerene 
collisions can be carried out with relatively little experimental effort due to 
the fact that molecular beams of mass-selected carbon clusters can be 
produced by simply evaporating purified fullerene material. In this Chapter 
we discuss the experimental and theoretical work involving the study of 
molecular fusion and subsequent fragmentation in fullerene-fullerene 
collisions. We also discuss the other reaction channels that can be observed 
from these collisions such as scattering, fragmentation without fusion and 
charge transfer. Fig. 9.1. 
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Figure 9.1. Snapshots from MD-DFT simulations of CC + C^o collisions illustrating the 
different observable reactions channels in fullerene-fullerene collisions. 



Apart from the very fundamental collision dynamics aspects, fullerene- 
fullerene collisions are related to other areas of fullerene science and 
technology where aggregation of fullerene molecules plays a role, e.g. 
aggregations in the laser plume during laser desorption [2,3], in fullerene- 
fullerite surface collisions (Chap. 10) and during polymerisation of 
fullerenes [4]. In addition, insight into the growth of larger fullerenes may be 
obtained from fullerene-fullerene fusion studies. The experiments and 
theoretical studies described in detail in this chapter are carried out under 
well-defined (single collision) conditions in order to provide information on 
the relevant dynamical processes in terms of reaction cross sections. This 
information can then be applied in order to improve the understanding of the 
above-mentioned processes which generally occur under less well-defined 
conditions. 

Before going into the detailed discussion of the topic, it is useful to get a 
general picture of possible reaction channels in fullerene-fullerene collisions 
and of the collision energies of interest. Figure 9.1 shows some snapshots 
from MD-DFT calculations (Chap. 3.3.2) illustrating the reaction products 
that can be expected in collisions between a singly charged and a 
neutral Ceo at different centre-of-mass energies. For low collision energies 
the main reaction channels are charge transfer and inelastic scattering. For 
low impact parameter collisions, the two fullerenes are compressed and 
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deformed on collision and "bounce" apart again with high vibrational 
excitation. Molecular fusion is a significant reaction channel for collision 
energies of around 100-200 eV where a metastable fullerene-like product is 
formed. As the collision energy increases beyond this, the increasing amount 
of vibrational excitation in the fusion product leads to increasing sequential 
evaporation of particles, even on the ps time scale of the simulations. At 
energies > 200 eV the reaction product fragments into multiple fragments on 
the 100 fs to 1 ps time scale. 

By combining experimental results with simple phenomenological 
models, statistical modelling, classical MD and detailed but time-consuming 
MD-DFT simulations it is possible to obtain considerable insight into the 
dynamics of these collisions. This is the topic of the following section. All 
collision energies given will be in the centre-of-mass reference frame, unless 
otherwise stated. 



1. FUSION AND FRAGMENTATION 



1.1 Fusion 

The first observation of molecular fusion between fullerenes was 
reported in 1993 [5] and was the starting point for a series of experiments on 
the reaction dynamics in low energy (50eV - 6keV) fullerene-fullerene 
collisions [6-11]. The experimental work was accompanied by molecular 
dynamics simulation studies by the Schmidt group and others [12-20]. The 
experiments were performed in the apparatus shown in Fig. 2.8. The 
positively charged products of the collision reaction were detected by a time- 
of-flight reflectron mass spectrometer. The reflectron could be rotated 
around the scattering cell allowing the determination of the angular 
distribution of ions. This was not considered in the early experiments [5-7] 
and in these papers the cross sections for fusion were underestimated at high 
energies where substantial fragmentation of the fusion product occurs. 
Newer measurements determined the angular distributions of the fusion 
products and integrated over these to obtain corrected values of the absolute 
cross sections [8,21]. A typical series of "raw" mass spectra is shown in Fig. 

9.2 for a scattering angle of 0±0.5°. The energy window in which a clear 
signal of a fusion product can be seen is rather small. It is possible to 
distinguish between fused reaction products and inelastically scattered 
projectile ions due to the kinematics of the collisions. When two fullerenes 
fuse together the velocity of the products and any fragments from a fused 
product are around half of the velocity of the parent ion. Fig. 9.3. This 
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time of flight/ gs 



Figure 9.2. Reflectron time-of-flight spectra of Ceo^ + C^o collisions at different centre-of- 
mass collision energies. At the lowest and highest collision energy shown there is no clearly 
discernible signal arising from fused fullerenes. Beyond the fusion threshold energy the 
product ion distribution broadens and shifts to lower masses as the collision energy is 
increased. Expected arrival times for the fusion products and some of their fragment ions are 

marked. Adapted from [7] 

allows them to be separated in a reflectron time-of-flight mass spectrometer. 
However, as can be seen in Fig. 9.2, the arrival times of small fragments 
from a fusion product overlap with large fragments from scattered 
projectiles. In order to separate out the different contributions to the time-of- 
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flight spectrum an additional pulsed electric field was used to accelerate and 
detect the products of the collision. This allowed all fusion products to be 
clearly separated from non-fused projectile ion fragments [9]. 



V, Vj = m^/(mi + mj) v., 




release 

•r fragmentation 



Figure 9.3. The velocity of the products of a fusion reaction between Qo^ and C^o is half that 
of the projectile ion plus the influence from kinetic energy release on fragmentation. 

The fusion cross sections are obtained from integrating over the intensities 
of all product ions with a laboratory velocity corresponding to that of fusion 
products, after corrections to the intensities have been made due to the 
detection probability at the micro-channel plate ion detector [7] (see Chap. 
2.1.3). The absolute fusion cross section for Ceo^ + Ceo collisions is plotted in 
Figure 9.4. Note that the error bars shown are relative errors in the 
measurement. There is a large error with respect to the absolute cross section 
due mainly to uncertainty in the fullerene vapour pressure. 

It is tempting to compare the form of the cluster fusion cross sections 
with those obtained in nuclear collisions. Fig. 9.4(b) gives a schematic 
picture of the behaviour in such collisions. The rise in the cross section at 
energies close to threshold can be well described using a simple line-of- 
centres or absorbing sphere model [22]. The fall in the cross section towards 
higher energies is due to the centrifugal instability of the fused compound. 
Fusion is an activated process and the simplest expectation is that the 
magnitude and energy dependence of the fusion cross section will be limited 
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Figure 9.4. (a) Measured absolute fusion cross section for + C(^ collisions. The 
experimental results are compared with a simple line-of-centres (or absorbing sphere) model 
which is often used to describe the fusion reaction in nuclear heavy ion collisions. Note that 
the calculated cross sections have been scaled by a factor of 0.07. 

(b) schematic cross section behaviour for nuclear collisions. 

by the condition that the energy available along the line of centres of the two 
fullerenes (or colliding nuclei) exceeds this threshold. At an initial relative 
translational energy E and impact parameter b the energy along the line of 
centres is £ = {\-b^/R^) where R is the centre-to-centre distance. We locate 
the barrier to fusion at a distance R equal to the sum of the radii of the 
collision partners R = cl ~ 7A for two Ceo- Then the maximum impact 
parameter that allows the crossing of the barrier is determined by 
E = {l-b„J/R^) = Ebfus or 
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(9.1.1) 



By plotting the fusion cross section as a function of 1/E it is possible to 
determine the energetic threshold for the reaction, assuming the collision can 
be well described in this way. However, the measured fullerene fusion cross 
section actually differs from the simple model in several key aspects. The 
rise in threshold is not quite linear in 1/E, Fig. 9.5, and the slope is quite 
wrong. The geometrical forefactor is much smaller than mcf. It is well known 
that the line-of-centres model typically overestimates the magnitude of the 
cross sections for chemical reactions. This is often ascribed to steric 
requirements of the reactions and the simplest correction is to lower the 
probability of reaction by a "steric factor" p, p<\ [6,7,23]. The forefactor is 
then Tofp and fitting values of p to the experimental results (after correcting 
for the scattering losses, not considered in the original articles) yields values 
from 0.07 for C 6 o "t"C 6 o 1^ 6.4 for C 70 +C 60 C 60 ^tnd 0.45 for C/o +C 70 

collisions. A more refined approach makes the steric factor also depend upon 
the impact parameter but to compute it requires knowledge of the barrier as a 
function of the approach geometry. The experimental threshold energies for 
the various fullerene fusion reactions were determined by assuming that the 
steric factor p increases linearly with collision energy. This gives a 
functional form 

[E-E,J (9.1.2) 



which provides a very good fit to the experimental data [7]. The 
experimental fusion threshold energies, determined by fitting eq. (4.2) to the 
data, are tabulated in Table 9.1. Such a dependence was rationalised by 
assuming that different orientations of the fullerenes on impact have 
different energetic barriers to fusion. 

Various molecular dynamics studies have been applied to determining 
the threshold energy for fullerene fusion [14,20,23,24]. The results of these 
studies are also given in Table 9.1. Note that there are also rather large errors 
in the MD determinations of the threshold energy due to the limited number 
of trajectories calculated but these have not been given explicitly. The 
simulations all give rather good agreement with experiment as long as the 
initial internal energy of the projectile ions is properly accounted for. 
However, extensive MD-DFT simulations [11,23] showed no evidence for a 
dependence of the fusion probability on the relative orientation of the 
colliding cages. The dependence on the initial orientation was tested in two 
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Figure 9.5. Absolute cross section for fusion in Qo"^ + C(,q collisions plotted as 1/E. The 
threshold energy for the reaction is determined by fitting eq (9.1.2) to the experimental data, 
dotted line. The simple line-of-centers model (with cut-off due to centrifugal instability), 
dashed line (linear dependence on 1/E close to threshold) does not give a good fit to the data. 

From [21]. 



different ways. First, two fullerenes were oriented hexagon-to-hexagon (6x6) 
or pentagon to pentagon ( 5 x 5 ). Then one of the fullerenes was rotated 
around the collision axis by the angle Sd. The collision trajectory was 
calculated for each value of dOaaA for several collision energies [ 8 , 11 ]. The 
probability for fusion was found to be more dependent on the collision 
energy close to threshold than on the initial orientation. Results for + 
Ceo collisions calculated with MD-DFT are given in Table 2 . No correlation 
with orientation can be found but the number of different orientations that 
lead to fusion does increase with increasing collision energy. Secondly the 
principal orientation of C70 was chosen prior to C70 + Ceo collisions [ 23 ]. The 
simulations showed no dependence on the orientation of the symmetry axis 
of C70. However, the absolute cross sections calculated with MD-DFT are in 
excellent agreement with the experimental data in contrast to those 
calculated with classical MD [8]. The absolute cross sections for the 
different collision systems studied are shown in Fig. 9.6 along with the 
results of the MD simulations. Unfortunately, not enough trajectories have 
been calculated for Ceo ^ + Ceo using the MD-DFT approach to give a reliable 
estimate of the fusion cross section. An early report gave an average 
calculated fusion probability close to threshold of 0.06 which is in agreement 
with the experimental values [6]. Although classical MD gives a reasonable 
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value for the threshold energy it vastly overestimates the absolute magnitude 
of the cross sections and also the width of the collision energy window for 
which fusion is an observable reaction channel. 



Table 9-1. Threshold energies for molecular fusion between fullerenes in units of eV 



System 


Exp. 

Threshold [7] 


Classical MD 
[8] 

T = 2000 K 


Classical MD 
[14] 

T = 2000 K 


MD-TB [20] 
T = 2000K 


MD-DFT 
[8,23] 
T=2000 K 


C 60 +C 60 


60+1 


60 


50 


60 


60 


Cvo+Ceo 


70+7 


70 






70 


C70+C70 


76+4 


80 






75 



Table 9-2. Some results of 6x6 orientation studies for C^o* + Cgo collisions using MD-DFT. 
One of the fullerenes was rotated about the collision axis by an angle Jff. Impact parameter 
b=0. Here "f" stands for fusion, "s" for scattering without fusion. From [11], 



66 


E=110eV 


E=120 eV 


E=130eV 


E=140eV 


E=150eV 


0.0 
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2.0 
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2.5 
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3.0 
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3.5 
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f 


4.0 
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f 


f 



Beyond the threshold, the cross section rises rather sharply to a maximum 
and then falls steeply to higher energies. In nuclear physics this decrease in 
the probability of fusion is associated with the compound nucleus not being 
formed. The impact parameter becomes so high that the effective potential 
V(R) + eFir^ no longer supports a well. MD simulations show that the 
product of the collisions close to the fusion threshold takes the form of a 
peanut, Fig. 9.7. The estimated stability of the peanut shaped isomer of Cno^ 
is about that of the reactants [25]. So Cno has a barrier for dissociation 
which about the same as the fusion barrier. To fill-in the well to the left of 
the barrier the centrifugal energy must balance the dissociation barrier which 
gives the condition for the impact parameter 

^2 ^ (9.1.3) 

E 



The energy dependence of the cross section due to such a cut-off is then 

(9.1.4) 
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Here too, the forefactor is much too large and the decrease as 1/energy is too 
moderate to fit the steep drop of the experimental cross section (Fig. 9.5). 
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Figure 9 . 6 . Cross sections for fusion between two fullerenes (a) + C^o', (b) + C^o and 

C^o*^ + C6o; (c) €70^ + C70. Triangles: experimental data; filled circles: MD-DFT simulations; 
open circles (right hand axis): classical MD simulations. Adapted from [8]. 
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0.12 0.48 0.72 0.96 1.2 



1 1 ps 

Figure 9. 7. MD-DFT simulation of a collision between Cjo* + Cso for a collision energy of 
130 eV and b=0. A metastable peanut isomer is formed on a time scale of 1 ps. From [26] 

In order to at least partially explain the discrepancies with the simple 
model one has to consider the possibility of competing reaction channels 
[21]. In the threshold regime the available phase space for the inelastic 
scattering channel is larger than that for the formation of a very hot 
compound fullerene, say Such a competing channel can be observed in 
molecular dynamics simulations, Fig. 9.8, and is also seen experimentally as 
will be discussed in more detail in the following section. To estimate the 
branching ratio for fullerene-fullerene collisions we require the density of 
vibrational states of a polyatomic molecule. For our purpose a simple RRK 
estimate is sufficient because we only need ratios and not absolute numbers. 
For an n-atomic species with 5 vibrational modes, s = 3N - 6, and a 
(geometric) mean vibrational frequency v (2.7x1 Hz = 900 cm ' for C^j,) 
the number of vibrational states below the energy E is 

1 f £ Y 

N{E)^~ (9.1.5) 

The essential difference between the two channels is that C 120 has many 
more vibrational modes than does Ceo- Once the energy is high enough to 
reach the vibrational quasi-continuum of C 120 , it will have overwhelmingly 
more states than Ceo- But it takes less energy to reach the quasi-continuum of 
Ceo- At the same internal energy, the ratio of the number of states is given by 




A^(Qo) 

A^(Q2o) 



2 (^+ 3 ) ’ 



5 = 174 



(9.1.6) 
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Figure 9.8. MD-DFT trajectory showing inelastic scattering where one of the fullerenes 
emerges from the collision with a high vibrational excitation. The time scale of the simulation 

is 1 ps. From [11] 

The mean energy per vibrational mode needs to exceed the mean 
vibrational spacing for C 120 to offer more vibrational states to accommodate 
the internal energy. In the immediate post threshold region we thus expect 
that many collisions do not end up in fusion but in inelastic collisions. At 
higher energies the situation is reversed and, because of the steep rise of the 
density of states with energy, once formation of C 120 is the dominant 
channel, it is the overwhelmingly dominant one. This very simple 
consideration leads to a much steeper rise of the cross section close to 
threshold [21], assuming an energy transfer of 0.3£ in the collisions leading 
to scattering rather than fragmentation (see 9. 1 .2 for a justification of this 
value). The increasing phase space that becomes available for fusion at 
energies past the barrier means that the fusion cross section rises more 
rapidly than as , as seen in figures 9.4, 9.5 and 9.6. However, the 

experimental situation is rather more complicated with a broad but unknown 
internal energy distribution for the projectile ions and also a distribution in 
the energy transferred in the scattering collisions. It is only possible to give 
an upper limit to the internal energy of the projectiles by considering the 
amount of metastable fragmentation in the beam. If one assumes that the 
fusion barrier shifts with the initial internal energy of the projectile and 
target, as indicated by the results of MD simulations, then it is also possible 
to obtain an overall cross section with positive curvature but for reasonable 
assumptions the rise in the cross section at threshold is not as steep as that 
observed in the experiment. 

Fusion is operationally defined as formation of a compound species 
accompanied by an extensive dissipation of the internal energy. It is however 
possible for the hot compound species to fragment promptly. This process is 
very clearly seen in molecular dynamics simulations where many trajectories 
leading to C 120 formation also show it breaking apart on a time scale of one 
or two vibrations. For such short times energy is not fully dissipated and 
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these events are not detected in the experiments as 'fusion'. The competition 
with the prompt fragmentation (with incomplete energy equilibration) means 
that the fusion cross section will decrease strongly with energy, past the 
threshold for three-particle production. What we therefore require is an 
estimate for the cross section for direct (i.e. prompt) collision induced 
dissociation (CID). The CID cross section can be computed if all the excess 
energy is disposed of as kinetic energy of the three outgoing particles, The 
result is [21] 



= A 






(9.1.7) 



where A is a geometrical factor and Ecm is the threshold for CID via fusion. 
The branching ratio for fusion is (cr -cr„„,)/(j where cr is the cross 
section for complex (=compound species) formation, given by the capture 
cross section of equation (9.1.1). The fusion cross section is 



*^fus *^2-body ( ( *^com *^CID ) ^ ^com ) 



(9.1.8) 



where a2-body incorporates the centrifugal cut-off, equation (9.1.4) and 
accounts for the competition with scattering close to threshold, as discussed 
above. The rapid decline of the fusion cross section beyond Ecm as given by 
equation (9.1.7) is in good agreement with the experimental results. This is 
shown in Fig. 9.9 for Cvo"^ -i- C70 collisions (the calculations do not take the 
range of initial internal energies into account which would presumably give 
a more "gentle" behaviour at the wings). The model works equally well for 
Cbo^ + C70 and €70^ -1- Cm collisions but still overestimates the + Ceo 
cross section by a factor of 6 although the form of the cross section is 
correctly reproduced. The reason for the exceptionally small cross section in 
Ceo^ + Cbo collisions is still not fully understood. 
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Figure 9.9. Comparison of experimental cross section for Cyo* + C^o collisions with fusion 
models. Triangles: experimental data. Dashed line: line-of-centres model with centrfugal cut- 
off. Thin full line: line-of-centres model scaled with p = 0.45. Thick full line: fusion cross 
section with competing scattering and direct collision induced dissociation channels. Adapted 

from [21]. 

The experimental and theoretical studies discussed in this sub-chapter 
were carried out under single-collision conditions. Coalescence and fusion of 
fullerenes also occurs under less well-defined conditions. Irradiation of 
fullerene films with laser radiation leads to the formation of fullerene 
aggregates up to the pentamer [3] with maxima in the corresponding mass 
spectra at 118, 178, 236 and 300 carbon mass units. These aggregates are 
very stable, surviving surface reflection at 160 eV collision energy. It can 
therefore be speculated that they have structures similar to those molecules 
formed in the collision experiments. It was also shown that in laser 
desorption, fusion occurs predominantly between Ceo and a Csg fragment. Cjg 
is chemically more reactive [27] but also less stable than Ceo- Thus, the 
fusion barrier in Csg-Ceo collisions is expected to be lower than in Ceo-Ceo 
collisions. The relative collision energy in the laser desorption experiments 
is much less than in the single-collision experiments. However, the 
interaction time is longer and the total available energy (kinetic plus 
internal) is similar in both cases. Molecular fusion of fullerenes was also 
observed in collisions between Ceo^ ions and fullerene films at 275 eV 
impact energy [28]. It is very likely that the reaction mechanism in these 
experiments is similar to the gas phase collisions. High mass products of a 
fusion reaction at this relatively high collision energy are possible since the 
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substrate can absorb energy, thus stabilising the fused compound and 
reducing fragmentation. 

1.2 Fragmentation 

At collision energies slightly above threshold, the fusion product can 
fragment on the microsecond time scale of the experiment [9]. The average 
fragment ion mass and the range of fragment ions produced from the fusion 
product of Ceo^+Ceo collisions is plotted in Fig. 9.10 as a function of internal 
energy (centre of mass collision energy plus initial maximum internal 
energy). The maximum entropy procedure (Chap. 3.2.2) has been used to 
characterise the outcome of the collision at these higher energies [21]. These 
products are classified as coming from fusion because they result from a 
process where the initial centre of mass energy has been converted to 
internal energy. The application of the maximal entropy model to the 
products of collisional fusion is much more straightforward than for laser or 
collisional excitation of Ceo [29] since we have a much better knowledge of 
the energy transferred to internal degrees of freedom in the collision. The 
same values of binding energies and ionisation potentials for all possible 
fragments up to and including are taken to be the same as those used 
earlier for applications of the model to Ceo and fragmentation [29]. 
Extrapolation to larger clusters was done by using the following expressions 
for the binding energy (BE(n)) and ionisation potential (IP(n)) of cluster C„ 
in units of eV 



5£(n) = -7.30n + 20.241 (9 P 9 ) 

IP(n) = 5.6 + 13.8/ yfn 

These expressions are based on [30] but the binding energy of the larger 
clusters has been reduced slightly to obtain a smoother dependence of 
binding energy on n on going from n < 60 to n > 60. 

The comparison of the experimentally determined average fragment mass 
as a function of total internal energy (assuming a total initial internal energy 
of 20 eV for the projectile ion and 5 ev for the neutral target) [9] with the 
maximal entropy model is shown in Fig. 9.10 and is rather good. The range 
of detected fragment ions is shown in Fig. 9.11. The detection limit for 
fragments predicted by the calculations was taken to be 0.001 (i.e. if the 
predicted fragment intensity was less than 10'^ of the initial parent ion 
intensity it was not counted). The agreement here is not quite so good. The 
behaviour at low collision energies follows the experiments very well and 
the onset of the bimodal fragment distribution is in good agreement. The 
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Figure 9.10. Average fragment size from the fusion product in + Cgg collisions 
versus total internal excitation energy (centre-of-mass collision energy plus estimated initial 
internal energy of projectile and target). Squares: experimental data. Full lines: result of 
maximal entropy calculation. Adapted from [21]. 

main discrepancy is that the model shows the presence of Cuo up to high 
collision energies beyond the onset of the bimodal fragment distribution. 
This is not seen in the experiments where there is a steady decrease in the 
size of the fullerene-like fragments up to a total internal energy of ca. 200 
eV. This is not too surprising since the binding energies used in the 
calculations are those of the most stable fullerene isomers. The agreement 
can be improved slightly by making some additional assumptions concerning 
the binding energies of the isomers formed in the collision experiments [21]. 

Overall the agreement between experiment and the simple statistical 
model is very good. This provides us with a confirmation that there is 
equilibration of the centre-of-mass collision energy among the vibrational 
degrees of freedom of the product of the fusion reactions and that the 
subsequent fragmentation behaviour is completely statistical in nature. MD- 
DFT calculations have also provided evidence for the efficient energy 
equilibration on a time scale of ps [23]. 
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Figure 9.11. Fragment size distribution as a function of total internal energy. Squares: 
experimental data, thin lines indicate range of detected fragment sizes (not error bars). Thick 
lines: maximal entropy calculations. Adapted from [21]. 

The picture of a purely statistical fragmentation mechanism is borne out by 
studies of the angular dependence of the fragmentation products [9]. Fig. 
9.12 shows angular distributions measured for the products of a fusion 
reaction at three different collision energies. The signal has been integrated 
over all fragment masses. The distributions are concentrated at small 
laboratory scattering angles. Assuming that the fused compound formed in 
the collision undergoes a very rapid energy equilibration and then undergoes 
a completely statistical fragmentation it is possible to estimate the expected 
cut-off angle as a function of collision energy. A Newton diagram for the 
collisions is shown in Fig. 9.13. Immediately after collision, the velocity of 
the fused compound in the centre-of-mass reference frame will be zero. The 
final centre-of-mass velocity of the fragment ions (w'ker) will be given by 
the accumulated kinetic energy released in the fragmentation process. This 
will be randomly distributed in space. The maximum detectable laboratory 
scattering angle (d„ax) is thus determined by the magnitude of the kinetic 
energy release which, in turn, depends on the temperature of the fused 
compound and the mass of the fragment. It was assumed that the small 
fragment ions (n < 30) are produced by a multifragmentation process in 
which the highly excited fused compound disintegrates into a number of ring 
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and chain fragments, one of which carries the positive charge and is detected 
in the experiment. 
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Figure 9.12. Integrated fusion signal intensity in Co* + Co collisions, (I( 0)sin( 0}), plotted as 
a function of scattering angle. Collision energies of (a) 122 eV (b) 131 eV (c) 165 eV. The 
energies shown on the plot are the total internal energy (centre of mass collision energy plus 
initial internal energy of projectile ion and target) Intensities have been normalised to the 
same maximum value. The arrows indicate the estimated maximum scattering angle for the 
maximum and minimum sized fullerene-like fragments detected in the experiment. The 
horizontal bar indicates the estimated maximum scattering angle range for the small fragment 
ions detected at this collision energy. From [9] 
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Figure 9.13. Newton diagram illustrating the relationship between laboratory and centre-of- 
mass frames of reference for fragmentation following fusion of Qo* + C(^. 



This may undergo a further evaporative step by emitting a C 3 molecule (the 
most energetically favourable decay pathway for small carbon clusters). The 
large fullerene fragment ions are assumed to be produced by successive 
evaporative loss of C 2 molecules. Obviously this is an over simplification- 
the small fragments may undergo further evaporative cooling on the time 
scale of the experiments and the large fragment ions may well be produced 
via the loss of neutral fragments larger than C 2 , especially at high collision 
energies. However, to a first approximation these seem to be reasonable 
assumptions. In order to estimate the temperature of the fragments one can 
make use of the maximal entropy calculations, discussed above. Fig. 9.14 
gives the relationship between total internal energy of the fused compound, 
<E> and the fragment temperature T, as obtained from the maximal entropy 
model. By treating successive evaporation processes as a random walk, 
whereby the net velocity after a chain of fragmentation steps is given by 
Wker = ^/iv iv with N the number of steps and vv the average velocity gained 
by a fragment ion in a single fragmentation, it is possible to estimate the 
maximum possible laboratory scattering angle from Fig. 9.13. Considering 
just the maximum and minimum sized fullerene fragment ions detected in 
the experiments at each collision energy, one can thus estimate the 
maximum expected laboratory scattering angles. These positions are marked 
with arrows in Fig. 9.12. The simple estimate gives rather good agreement 
with the experimental data for the two lowest collision energies. At the 
highest energy it is likely that the assumption of purely successive C 2 
fragmentation breaks down and the maximum laboratory scattering angle is 
therefore over-estimated. Similar estimates can be carried out for the small 
fragment ions although they are of low intensity for these collision energies 
and do not contribute significantly 
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Figure 9.14. Relationship between average internal energy of fused compound, C 120 *, and the 
temperature of its fragments, calculated with the maximal entropy model (Chap. 3.2.2). The 
onset of the "phase transition" corresponds to the threshold for observation of small (n > 30) 

fragment ions. 

to the overall scattering intensity. The horizontal bar on Fig. 9.12b indicates 
the range of maximum scattering angles at this collision energy for the 
observed small fragment ion mass range of n = 15-23. 

The angular distributions of the integrated fragment ion intensities from 
scattered projectiles that have not undergone a fusion reaction are shown in 
Figure 9.15 and are quite different from the distributions from the fusion 
products in Fig. 9.12. The corresponding average fragment ion sizes are 
plotted in Fig. 9.16. The angular distribution is centred around a laboratory 
scattering angle of about 8° for the lowest collision energy investigated of 
130 eV. It then broadens and moves towards smaller angles, reaching its 
maximum width for a centre-of-mass collision energy of 860 eV. As energy 
is increased beyond this the angular distribution becomes narrower again. 
Measurements of the ion kinetic energies, using a retarding field energy 
analyser, show that the fragment ions velocity is about 85% of the initial 
projectile velocity [9]. With the help of the Newton diagram in Fig. 9.17 it is 
possible to estimate the amount of energy converted to internal excitation of 
the projectile and target for projectiles scattered at an angle corresponding to 
the maximum in the measured distributions. For all collision energies 
investigated this amounts to close to 60% of the centre-of-mass energy. This 
is in good agreement with molecular dynamics statistical trajectory ensemble 
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Figure 9.15. Integrated fragment intensities from scattered projectile ions at different collision 
energies as a function of laboratory scattering angle. Adapted from [9] . 

simulations at a collision energy of 500 eV that indicated the transfer of ca. 
50% of the centre-of-mass energy to internal degrees of freedom [31]. If we 
further assume that on average this energy is equally divided among the 
projectile and target and we know the average fragment size from Fig. 9.16 
then it is also possible to estimate the width of the angular distribution. We 
follow the procedure discussed above for the fragments from the fused 
product. We assume the projectile ions are first scattered inelastically and 
then subsequently undergo statistical fragmentation. The fragment 
temperature is obtained from Fig. 3.10. The estimated widths are in very 
good agreement with the measurements. They increase from ±2° to ±5° on 
increasing the collision energy from 132 to 500 eV. For collision energies of 
750 eV and 860 eV where one still has some fullerene-like fragments as well 
as small fragment ions (average fragment ion size of around n = 30, Fig. 
9.16) the angular distribution is estimated to be ca. ±1°. On further 
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increasing the collision energy beyond this range, the fragment ion mass 
spectrum consists of only small fragment ions (where we consider only two 
steps, as discussed above, to estimate the angular distribution - a 
multifragmentation followed by one evaporative cooling step) and the 
estimated angular distribution decreases again to about ±5°. 




Figure 9.16. Mean fragment ion size from scattered projectile ions as a function of centre of 
mass collision energy. The bars show the range of detected masses (not error bars). Adapted 

from [9]. 




Figure 9.17. Newton diagram for Cai* + C^o collision with inelastic scattering, v: velocity in 
laboratory system; w: velocity in cente-of-mass system; p: projectile; t: target; w'ker: velocity 
due to kinetic energy release on fragmentation. 

Both the fragmentation of the fused compound and the fragmentation of 
the inelastically scattered projectile ion are consistent with a completely 
statistical fragmentation at the low collision energies relevant to these 
experiments. At the high energy range of the investigated collision energies 
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the fused compound has an extremely short lifetime and breaks up into 
multiple fragments before the centre-of-mass collision energy can be 
statistically distributed among all vibrational degrees of freedom. This is the 
reason for the sharp fall in the fusion cross sections. 

The experiments discussed in this section provide convincing evidence 
for the occurrence of a phase transition and the onset of multifragmentation, 
even at the relatively low collision energies investigated experimentally. 
However, the experiments are restricted since only the positively charged 
fragments can be detected at the low collision energies involved and this 
makes an unambiguous detection of multifragmentation impossible. The 
situation is different for very high collision energies in the MeV range where 
there is also a high probability for detection of neutral fragments. 
Coincidence experiments have been used to investigate cluster-cluster 
collisions under such high collision energies, using Ceo as the target for H2s'^ 
projectile ions [32,33]. The observed power law in the fragment distributions 
with increasing collision energy along with the increasing multplicity is a 
strong indication of a liquid-gas phase transition such as discussed by 
Bonasera and co-workers [34]. 



2. CHARGE TRANSFER 

The study of charge transfer in collisions where fullerenes are both the 
projectile ion and the target has been carried out for collision energies 
ranging from thermal to 100 keV and for projectile charges ranging from -i-l 
to -1-6 [35-40]. These collision systems have been of interest for exploring the 
application of simple phenomenological charge transfer models, borrowed 
from atomic physics, such as the over-the-barrier model (Chap. 3.4). The 
results also help to provide a detailed understanding of the electronic 
properties of the isolated fullerenes and the conditions for electrical contact 
that may be highly relevant in future for producing nanoscale electrical 
devices based on fullerenes and carbon nanotubes. 

The measured cross sections for collisions between fullerene ions and 
neutral fullerene targets at a laboratory collision of 100 keV are given in 
Table 9.3 [36]. Here cr^^. refers to collisions with a q+ charged projectile ion 
that result in a q'+ charged product. The single capture cross sections (027 
and (T 32 ) are very similar in magnitude. This can be explained with the 
simple classical over-the-barrier argument, in this case using the "movable 
hole" picture (Chap.3.4.1). In this case, the high polarisability of the neutral 
fullerene target will tend to orient the initially delocalised holes on the 
projectile ion so that one positive charge will be at the minimum distance to 
the target. At the same time, the target will itself become polarised so that if 
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an electron is transferred, it will start out from the point on the target that is 
closest to the projectile. In a first approximation, the electron transfer from 
the target to the projectile can be modeled as a singly-charged ion-atom 
collision, with an ion-atom separation equal to the minimum distance 
between the fullerenes. The classical barrier argument then yields a critical 
separation between the fullerenes, Rc = llao- The cross-section for single- 
electron capture then becomes: 



a{a-\) 



q + \ 



-n 






(9.2.1) 



where is the fullerene radius. This gives values of (J21 = 280 or 240 
calculated for = 8.2 ao and 6.7 Uq, corresponding to the geometrical cage 
radius and to the value obtained from the polarisability, respectively. The 
corresponding values for (732 are 310 A^ and 270 A^. The agreement with 
experiment is really very good. 

To explain the double-electron capture by a triply-charged projectile one can 
consider the relation between the geometry of the two positive charges left 
on the projectile after capture of the first electron and a positive charge now 
on the target. Fig. 9.18. Since the two ions are close to one another there are 
two different possible configurations. The first is that the positive charge is 
repelled by the projectile and moves to the far side of the fullerene cage. 
This is the most likely situation for large impact parameters and it does not 
lead to any further electron capture. The second possibility, which is only 
stable for small impact parameters, is where the three positive charges all fall 
on a straight line passing through the fullerene centres. Fig. 9.18. This is 
possible because there is a polarisation attraction on one of the positive 
charges of the projectile by the target. The competition between this 
attraction and the Coulomb repulsion from the target positive charge can be 
estimated to give a critical distance for further charge transfer of Rc ~ 7ao. 
The classical over-the-barrier cross section is then 



— — K 
2 






(9.2.2) 



which gives 100 A^ or 120 A^ depending on the value used for r^, again in 
reasonably good agreement with the experimental values [36]. 
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Table 9-3. Experimental cross sections (Tq,- for charge transfer in units of for collisions 
between fullerenes at a laboratory collision energy of 100 keV. q and q' refer to the projectile 
ion charge before and after collision. Data from [36]. 
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Figure 9.18. Schematic picture of double electron transfer to a triply charged fullerene form a 
neutral fullerene. Adapted from [36]. 

The cross section for double electron transfer to a singly-charged 
projectile to form the anion is orders of magnitude smaller than the other 
charge transfer reactions and cannot be estimated with the simple classical 
over-the-barrier arguments. 

Studies of collisions between a neutral Ceo target and charged projectile 
ions, € 60 “*^ (q = 2-5) at collision energies of 2q keV confirmed the basic 
picture for the charge transfer mechanism [40]. An extension of the simple 
model of Fig. 9.18 to higher projectile charges predicts that the maximum 
observable target charge after collision with projectiles of q = 1,2, 3, 5 should 
be 1,1, 2,2. This is in perfect agreement with the experimental results. The 
results also show that electrical contact between two isolated, charged Ceo is 
already established at a centre-to-centre distance in excess of 20 ao and that 
the charge re-distribution must occur on a time scale shorter than lO '"^ to 
10 *^ s. 

The values shown in Table 9.3 are for a relatively high collision energy. 
However, the absolute values of the cross sections do not change very much 
with collision energy. Cross sections are plotted as a function of collision 
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energy for much lower energy collisions in Fig. 9.19 and 9.20. Fig. 9.19 
shows the measured cross sections for C6o‘*^ + Ceo charge transfer collisions 
for q = 2,3. The value of the cross sections for single and double electron 
transfer to a projectile with q = 3 are ca. 300-400 and 100-200 A^, 
respectively in the collision energy range of a few keV [38], in very good 
agreement with the simple model discussed above. Similarly good 
agreement is obtained from single electron transfer collisions involving 
doubly charged projectiles where the measured cross section is scattered 
around 100-300 A^, to be compared with the model estimate of 240-280 A^. 
The low energy collisions have also been compared with the two-state Rapp 
and Francis model (Chap. 3.4.2), modified to take the size of the fullerene 
cage into account [37]. The results of these calculations are also shown in 
Fig. 9.19 and again there is good agreement, especially if one assumes a 
small off-resonance in the charge transfer reaction. 




collision energy [eV] 

Figure 9.19. Absolute charge transfer cross sections for Qo^'^+Ceo collisions as a function of 
centre of mass collision energy. Top: q = 'i [38] single and double electron transfer. Bottom: q 

= 2 [37] single electron transfer. 
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In all of these studies it is assumed that charge transfer occurs 
predominantly for grazing collisions and the experiments generally 
discriminate against low impact parameter collisions that would lead to 
scattering as well as transfer of excitation energy to vibrational degrees of 
freedom, leading ultimately to fragmentation. This is a reasonable 
assumption for high energy collisions and collisions with multiply-charged 
projectile ions. One has to be more careful with low energy collisions using 
a singly-charged projectile ion. Results of the measurements of charge 
transfer cross sections for single-electron transfer in collisions between Ceo^ 
+ Ceo are shown in Fig. 9.20. The original publication only presented results 
for products detected within a laboratory scattering angle range of 0 ± 0.6° 
[39]. Later investigations of differential cross sections showed that scattering 
outside this range was in fact significant for fragmentationless charge 
transfer reactions with singly charged projectiles and collision energies in 
the range of a few hundred eV to a few keV [10]. The original data is 
shown in Fig. 9.20 along with the cross sections obtained from integrating 
over the scattering angle dependence of the charge transfer signal as a 
function of collision energy. Clearly, the earlier cross sections were 
underestimated by a factor slightly larger than 2 for the investigated collision 
energy range. The correction brings the cross section value into very good 
agreement with the other single-electron transfer cross sections and the 




Figure 9.20. Total resonant charge transfer cross section for Qo *+C(^ collisions. Squares: 
original data not accounting for effects of scattering [39]. Circles: data corrected for scattering 
[10]. Dotted lines: upper and lower limits estimated from the Rapp and Francis model [39]. 
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prediction of the over-the-barrier model (210-240 A^). It is also in good 
agreement both in terms of magnitude and in terms of energy dependence for 
low energies with the prediction of the Rapp and Francis model for resonant 
charge transfer [39]. 



REFERENCES 

1 R. Schmidt and H. O. Lutz, Phys. Rev. A 45, (1992) 7981. 

2 R. Mitzner, B. Winter, C. Kusch, E. E. B. Campbell, and I. V. Hertel, Z. Phys. D 37, 
(1996) 89. 

3 C. Yeretzian, K. Hansen, F. Diederich, and R. L. Whetten, Nature 359 , (1992) 44. 

4 A. M. Rao, P. Zhou, K.-A. Wang, G. T. Hager, J. M. Holden, Y. Wang, W.-T. Lee, X.-X. 
Bi, P. C. Eklund, D. S. Cornett, M. A. Duncan, and I. J. Amster, Science 259 , (1993) 955. 

^ E. E. B. Campbell, V. Schyja, R. Ehlich, and I. V. Hertel, Phys. Rev. Lett. 70, (1993) 263. 

^ F. Rohmund, E. E. B. Campbell, O. Knospe, G. Seifert, and R. Schmidt, Phys. Rev. Lett. 
76,(1996) 3289. 

^ F. Rohmund, A. V. Glotov, K. Hansen, and E. E. B. Campbell, J. Phys. B. 29 , (1996) 
5143. 

^ A. V. Glotov, O. Knospe, R. Schmidt, and E. E. B. Campbell, Eur. J. Phys. D 16 , (2001) 
333. 

^ A. V. Glotov and E. E. B. Campbell, Phys. Rev. A 62 , (2000) 033202. 

A. V. Glotov and E. E. B. Campbell, Chem. Phys. Lett. 327 , (2000) 61. 

^ ^ A. V. Glotov, PhD Thesis, Gdteborg University, Sweden, (2000). 

O. Knospe and R. Schmidt, in Theory of Atomic and Molecular Clusters with a Glimpse at 
Experiments, edited by J. Jellinek (Springer Verlag, 1999), p. 111. 

L. Liu, K. Chen, and Y. Li, Chin. Phys. Lett. 10 , (1993) 720. 

1^ D. H. Robertson, D. W. Brenner, and C. T. White, J. Phys. Chem. 99, (1995) 15721. 

13 D. H. Robertson, D. W. Brenner, and C. T. White, Mat. Res. Soc. Symp. Proc. 359, (1995) 
199. 

ll’ R. Schmidt, J. Schulte, O. Knospe, and G. Seifert, Physics Letters A 194 , (1994) 101. 

1^ J. Schulte, O. Knospe, G. Seifert, and R. Schmidt, Physics Letters A 198 , (1995) 51. 

1^ Y. Xia, Y. Xing, C. Tan, L. Mei, and H. Yang, Nucl. Instr. Meth. Phys. Res. B 111 , 
(1996)41. 

19 Y. Xia, Y. Xing, C. Tan, and L. Mei, Phys. Rev. B 53 , (1996) 13871. 

20 B. L. Zhang, C. Z. Wang, C. T. Chan, and K. M. Ho, J. Phys. Chem. 97 , (1993) 3134. 

21 E. E. B. Campbell, A. Glotov, A. Lassesson, and R. D. Levine, C.R. Physique 3, (2002) 
341. 

99 

R. D. Levine and R. B. Bernstein, Molecular Reaction Dynamics and Chemical Reactivity 
(Oxford University Press, New York, Oxford, 1987). 

23 O. Knospe, A. V. Glotov, G. Seifert, and R. Schmidt, J. Phys. B. 29, (1996) 5163. 

24 R. Ehlich, PhD Thesis (Universitat Freiburg), (1994). 

23 D. L. Strout, R. L. Murry, C. Xu, W. C. Eckhoff, G. K. Odom, and G. E. Scuseria, Chem. 
Phys. Lett. 214,(1993)576. 

26 F. Rohmund, PhD Thesis, FU Berlin, (1997). 

22 S. Petrie and D. K. Bohme, Nature 365 , (1993) 426. 

2^ T. Lill, F. Lacher, H.-G. Busmann, and I. V. Hertel, Phys. Rev. Lett. 71, (1993) 3383. 

29 E. E. B. Campbell, T. Raz, and R. D. Levine, Chem. Phys. Lett. 253, (1996) 261. 

30 G. Seifert, R. Gutierrez, and R. Schmidt, Phys. Lett. A 211 , (1996) 357. 




9 . Fullerene-Fullerene Collisions 



189 



J. Schulte, Int. J. Mass Spec. Ion Proc. 145, (1995) 203. 

32 B. Farizon, M. Farizon, and M. J. Gaillard, Int. J. Mass Spectr. 192, (1999) 259. 

33 B. Farizon, M. Farizon, M. J. Gaillard, F. Gobet, C. Guillermier, M. Carre, J. P. Buchet, P. 
Scheier, and T. D. Mark, Eur. Phys. J. D 5, (1999) 5. 

34 M. Belkacem, V. Latora, and A. Bonasera, Phys. Rev. C 52, (1995) 271. 

33 D. K. Bohme, Int. Rev. Phys. Chem. 13, (1994) 163. 

36 H. Shen, P. Hvelplund, D. Mathur, A. Barany, H. Cederquist, and N. Selberg, Phys. Rev. 
A 52, (1995) 3847. 

32 F. Rohmund and E. E. B. Campbell, Chem. Phys. Lett. 245, (1995) 237. 

38 F. Rohmund and E. E. B. Campbell, Z. Phys. D 40, (1997) 399. 

39 F. Rohmund and E. E. B. Campbell, J. Phys. B 29, (1997) 5143. 

40 H. Cederquist, P. Hvelplund, H. Lebius, H. T. Schmidt, S. Tomita, and B. A. Huber, Phys. 
Rev. A 63, (2001)025201. 




Chapter 10 

SURFACE COLLISIONS 



"I hurried to complete the paper sphere... When I dropped it on the floor, 
it actually bounced. Surely, this had to be the correct shape." This quotation 
from Rick Smalley [1] comes from a description of the search for the 
structure of Ceo and can perhaps be regarded as the first report on a surface 
collision experiment with the molecule! Fullerene surface collisions received 
a considerable amount of interest after initial reports of the extremely high 
resilience to fragmentation of at relatively high impact energies [2]. The 
observation was also in agreement with the predictions of classical MD 
simulations for collisions between a fullerene and a diamond surface [3]. 
Many of the phenomena observed in gas phase collisions have their 
counterpart in surface collisions. A schematic summary of the processes that 
have been investigated is shown in Fig. 10.1. In this chapter we will review 
the field of fullerene-surface collisions. The emphasis will be on the basic 
scattering dynamics but we will also touch upon possible applications that 
exploit the unique nature of the fullerene collisions. 



1. FRAGMENTATION AND SCATTERING 



1.1 Eragmentation 

Most of the experiments that will be described in this chapter have been 
carried out in apparatus that is similar to that shown in Fig, 2.17. A beam of 
fullerene ions is produced (normally, laser desorption/ionisation has been 
used to produce the fullerene cations for the surface scattering experiments), 
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Metastable 
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Thermionic 
Electron Emission 



Backscattering of the 
Intact Cluster 
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Secondary Ions 



Deformation and Heating of 
and Target 



Figure 10.1. Schematic summary of the processes observed in collisions of Co+ ions with 

surfaces. 

mass-selected, accelerated or decelerated if necessary and collided with a 
surface (silicon, graphite and diamond have been popular choices). The 
scattered fullerene ions and their fragments are then detected by time-of- 
flight mass spectrometry. Early experiments reported that the surface 
collisions in the impact energy range of 50 - 500 eV were highly inelastic, 
resulting in a nearly constant kinetic energy of the scattered ions of 5-10 eV 
[2, 4-6]. A value of 25-35 % energy transfer to internal degrees of freedom 
was assumed based on classical MD simulations for collisions with diamond 
[3]. It was later shown [7] that the kinetic energy of the scattered fullerenes 
as well as their internal energy after the surface collision was strongly 
dependent on the cleanliness of the surface. The graphite targets used in the 
early experiments had a residual hydrocarbon coating that affected the 
collision dynamics. Time-of-flight mass spectra of the products of collisions 



10. Surface Collisions 



193 




Mass [u] 



157 




Mass [u] 

Figure 10.2. Time-of-flight mass spectra for collisions of Cgo* with a graphite target. 
Upper: with residual hydrocarbon adsorbates present on target. Lower: clean graphite. 

Adapted from [7] 

between Ceo^ and a clean graphite surface are compared in Fig. 10.2 with the 
products obtained for the same impact energies but a graphite target that has 
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hydrocarbon adsorbates present. The small masses present at short flight 
times from the target with adsorbates are secondary hydrocarbon ions. It is 
clear that more energy is transferred to internal degrees of freedom in the 
collisions with a clean surface. The fragmentation behaviour is very similar 
to that observed in the gas phase. For low excitation energies one sees the 
typical pattern due to successive evaporation of C 2 molecules. At impact 
energies between 256 and 353 eV the small ring and chain geometry carbon 
fragment ions appear in the spectrum. This is the same as for collisional (or 
laser) excitation in the gas phase and is the signature of a phase transition in 
the excited fullerene leading to the rapid break-up into small fragments as 
modelled with the maximal entropy formalism, discussed in Chap. 3.2.2. 
There is a rather large uncertainty in the initial internal energy of the 
fullerene ions in these experiments but if we consider the threshold energy 
for the appearance of the small fragment ions, as predicted by the maximal 
entropy model (ca. 65 eV [8]) we can easily estimate that approximately 20- 
25% of the collision energy must have been transferred to internal degrees of 
freedom in the collisions. The kinetic energy of the scattered fullerene ions is 
also affected by the presence or absence of the adsorbate layer. Fig. 10.3 
shows the velocity distributions for the scattered parent ions. The post- 
impact velocity is approximately twice as high when the graphite target has 
no hydrocarbon adsorbates [7]. A simple impulsive excitation model can 
qualitatively explain the observations. The model predicts that the energy 
transfer efficiency will increase with an increase in the effective surface 
mass that is determined by the mass of the functional group hit by the 
incident molecular ion. The observation that collisions with the densely 
packed carbon atoms of the clean HOPG (highly oriented pyrolytic graphite) 
surface result in a more efficient energy transfer than collisions with a 
disordered hydrocarbon layer having much lower average atomic mass and 
density is consistent with the model. 

In addition to the scattered parent and fragment ions and the sputtered 
hydrocarbon ions it is also possible to observe ions with masses larger than 
when scattering from the hydrocarbon-coated surface [7] (barely visible 
in Fig. 10.2). This can be attributed to "pick-up" reactions at the surface. 
The high internal energy of the scattered ion makes it possible for the 
collision complex to rearrange ("anneal") to form stable higher fullerenes. 
The addition of up to 5 C 2 units has been resolved in such experiments. 

It is also possible to detect small negatively charged carbon clusters as 
the products of collisions between positively charged Ceo^ and various 
targets (e.g. graphite, silicon, tantalum) [9]. Some results are shown in Fig. 
10.4 as a function of collision energy. The positive and negative integrated 
fragment ion yields are shown in Fig. 10.5 as a function of collision energy. 
The small negatively charged carbon clusters are neutral fragments from the 
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Figure 10.3. Measured velocity distributions for Ceo^ scattered from graphite at an impact 
energy of 200 eV. Left: target with residual hydrocarbon adsorbates. Right: clean target. 



Adapted from [7] . 

Anions 




Number of Carbon Atoms 

Figure 10.4. Small carbon anions observed from collisions between and a graphite target 
for different impact energies. The signal is due to small neutral fragments "picking-up" an 
electron from the surface. Adapted from [9] . 
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Figure 10.5. Integrated intensities of detected fragment species as a function of C(C impact 
energy on graphite. Adapted from [9] . 

prompt break-up of the fullerene cage into ring and chain species that pick 
up an electron from the surface. The high electron affinities of the small 
carbon clusters [10] makes them ideal candidates for dynamic electron 
transfer/pick-up processes at or near the surface (negative surface 
ionisation). Positive surface ionisation processes will be much less efficient 
due to the large ionisation potentials of the small clusters [11]. The electron 
transfer process can be via a surface harpooning mechanism [12]. An 
electron transfer is favourable for a molecule closer than a distance 
R^=e^ /4{(p-Ej from the surface, where (f) is the work function of the 
surface and Ea is the electron affinity of the molecule, due to the Coulombic 
stabilisation energy of the generated anion and the image charge on the 
surface. If the kinetic energy of the particle leaving the surface is greater 
than (f) - Ea then there is a non-zero probability that the negative ion can 
escape from the surface without being re-neutralised. The velocity 
distributions measured for the small carbon anions in these experiments give 
average velocities in the range 2000-8000 ms ' which, together with 
maximum distances of = 2 - 6 A for electron transfer, gives a sub- 
picosecond time scale for the break-up of the excited fullerene cage. This 
agrees very well with MD simulations of highly excited fullerenes [13-15]. 
The formation of negative parent ions can be seen in experiments involving 
hyperthermal collisions of neutral fullerene beams on surfaces [16]. 
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Tight-binding MD simulations [15] of Ceo collisions on a diamond 
surface indicated the presence of three different collision regimes. For low 
collision energies, £ < 120 eV, the molecule bounces off the surface without 
forming any bonds with the surface atoms and recovers its original shape 
after severe distortions. In the intermediate regime, 120 < £ < 240 eV, Ceo 
forms covalent bonds with the surface. The molecule can either be adsorbed 
or can leave the substrate with defects in its original structure, leading to 
metastable evaporation of C 2 units. Similar results have been found for 
simulations using a semi -empirical density functional approach [17], Fig. 
10.6. Finally, for £ > 240 eV, the molecule breaks after impact with the 
substrate, forming large ring and chain fragments. Pieces of the broken cage 
can form stable bonds with the substrate which can induce a local surface 




Figure 10.6. MD simulations of a C^o colliding with a diamond surface at an impact 
energy of 150 eV and an angle of incidence of 20°. Adapted from [18]. See also [17] for 

details of the calculation. 
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reconstruction. The calculations also showed that the kinetic energy of the 
fullerene centre of mass, when it bounces off the surface shows a clear 
dependence on the incident energy in the low collision energy regime. 
However, in the intermediate regime, the centre-of-mass energy after impact 
depends only weakly on the impact energy [15]. The calculations are in good 
qualitative agreement with measurements of Ceo^ scattered on diamond [18] 
in the collision energy range of 100-1500 eV. Velocity distributions for 
collision energies in the intermediate range are shown in Fig. 10.7. All 
distributions measured in this range (220 - 455 eV) show a velocity 
distribution that peaks at approximately 2000 ms ' and is thus independent of 
the impact energy. The intensity of scattered ions decreases strongly for 
impact energies greater than 300 eV and for energies beyond 500 eV, no 
scattered ions could be observed at all. This can be explained by complete 
destruction and reaction of the fullerene with the surface atoms. Such a 
process has been used to form a stable, resistant amorphous carbon coating 



Collisions with Diamond 
C0O fragments 




0 2 4 0 

velocity [km/s] 



Figure 10.7. Velocity distributions of scattered C^o* (left) and fragment ions (right) for 
collisions with a diamond (111) surface at different impact energies and an angle of incidence 

of 20°. Adapted from [18]. 






10. Surface Collisions 



199 



on air- sensitive fulled te films [19]. 

Lower energy collisions were investigated by decreasing the angle of 
incidence of the fullerene beam on the surface thus decreasing the normal 
component of the initial velocity. The experimental geometry is shown in 
Fig. 10.8. For these experiments it was not possible to rotate the time-of- 
flight mass spectrometer, as in Fig. 2.17, instead two fixed detector positions 
were used and the target could be rotated to vary the angle of incidence. The 
dependence of the final normal velocity to the initial normal velocity is 
plotted in Fig. 10.9 and shows an initial increase in the low collision velocity 
regime, reaching a plateau for velocities corresponding to kinetic energies 
beyond 120 - 140 eV, in good agreement with the MD simulations [15]. 

y= 140° y=45° 



Figure 10.8. Experimental geometry for the experiments described in [18, 20]. The positions 
of two time-of-flight detectors are fixed corresponding to the shown y angles. The angle of 
incidence can be varied by rotating the target. 

A similar dependence of the final kinetic energy with respect to the initial 
kinetic energy was shown in collisions between a hyperthermal neutral Ceo 
beam and a carbon-passivated nickel surface for impact energies in the range 
10-50 eV [21]. The kinetic energy of the scattered molecules is plotted as a 
function of the initial kinetic energy for two angles of incidence in Fig. 
10.10. A linear dependence of the scattered energy on the initial energy is 
clearly seen from the low impact energy data. 

1.2 Scattering 

In the above section we concentrated on the component of the initial 
kinetic energy of the projectile normal to the surface. Here we focus on the 
tangential component. The behaviour of the tangential velocity component is 
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Figure 10.9. Final velocity component normal to the surface as a function of initial velocity 
component normal to the surface (maxima of the velocity distributions) for impact. 
Squares 7=45°. Circles: "^140°. Adapted from [20]. 




Initial kinetic energy [eV] 



Figure 10.10. Relationship between initial and final kinetic energy of colliding on a 
carbon-passivated Ni surface. Squares: incident angle of 61°. Circles: incident angle of 15.7°. 

Data from [21]. 
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closely related to the rotational angular momentum transferred to the 
molecule. The fullerene behaves very much like a rubber ball picking up top- 
spin in a collision with a solid surface [20-23]. 

The collision of a Ceo molecule with a solid surface is illustrated 
schematically in Fig. 10.11. We assume that the Ceo hollow sphere touches 
the surface at any given time at just one point P only and in such a manner 
that no torque is applied with respect to that contact point. Assuming the 
absence of initial rotation, conservation of angular momentum gives a simple 
kinematic constraint 

amUi—amv^ + L ( 10 . 1 ) 

where u, and v, are the initial and final tangential velocity, respectively, a is 
the distance of closest approach, m the fullerene mass and L = I£2 the 
angular moment of the ball after the collision. I = qma^ is the moment of 
inertia at the time of closest approach with the numerical factor q = 2/3 for 
an undeformed hollow sphere. This can be taken as a good approximation 
for low energies. As the energy increases, the spherical molecule becomes 
deformed under the influence of normal forces into an ellipsoid as shown in 
Fig. 10.11. Equation (10.1) can be written as 

w^-u^-qaD. ( 10 . 2 ) 

which gives the relationship between the tangential velocities and the 
velocity of rotation ai2 at the surface of the fullerene. Equation (10.2) 
predicts that the final tangential velocity, Vt, will be smaller than the initial 
velocity component parallel to the surface, Uf Also, the loss of tangential 
velocity during the collision is predicted to increase for larger factors, q, i.e. 
with increasing deformation of the sphere. An upper bound for a£2 is given 
for complete sticking of the contact point to the surface, a£2 = Uf Then one 
obtains Vt = Ut {l-q), giving Vt = 1/3 for an undeformed sphere and a 
decreasing v/ut ratio at higher impact energies where q becomes larger. 
Eigure 10.12. shows the experimentally determined v/ut ratio for different 
angles of incidence and impact energies [20]. The experimental values are 
very close to 1/3 for low impact energies and start to deviate significantly to 
lower values for energies greater than 800 eV (corresponding to a normal 
impact energy of 120 eV), in very good agreement with the simple rolling- 
ball model. A similar model was applied to the experiments using a 
hyperthermal neutral fullerene beam [21]. A deviation from the simple 
model was observed which could be attributed to the existence of 
translational slip for the low energy neutral collisions [21]. 
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FigurelO.il. Schematic picture of the scattering geometry illustrating the transfer of 
tangential velocity into rotational angular momentum of the during the surface collision. 

Adapted from [22]. 
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Figure 10.12. Ratio of final and initial tangential velocities as a function of the angle of 
incidence for different impact energies. The data are clusteres around V( = u,/3 as expected 
from the simple rolling ball picture. Adapted from [20]. 





10. Surface Collisions 



203 



1.3 Collisions with Fullerite Films 

Fullerite films (epitaxial films of fullerene molecules) have three 
properties that are important for the Ceo-surface scattering process. Firstly, 
the projectile and the target molecules in the surface have the same mass. 
Secondly the surface corrugation is large and its period matches the diameter 
of the impinging fullerene. Thirdly, the surface consists of only weakly 
bound molecules. These three characteristics are responsible for the different 
results obtained from fullerene-fullerite collisions [18] compared to 
scattering of fullerenes from more conventional surfaces, discussed above. 

Fig. 10.13 shows a typical mass spectrum of ions after impact on a 
fullerite film [24]. In this particular example the initial kinetic energy of the 
projectile beam was 275 eV and the incident angle was 25 ° with respect to 
the surface normal. The products from colliding with a graphite surface 
under identical experimental conditions are shown in the bottom part of the 
figure. Practically no fragmentation is observed in the fullerite scattering 
case but there is a signal corresponding to masses much larger than the 
projectile ion. This is due to a molecular fusion reaction between the 
projectile fullerene ion and a target fullerene on the surface. The reaction is 
very similar to the gas phase fusion discussed in Chap. 9. One can clearly 
see a typical fullerene fragmentation pattern for the larger molecules 
corresponding to the loss of multiple C 2 units indicating that there is 
sufficient internal energy present for a rapid rearrangement of the carbon 
bonds to produce a metastable Ci 2 o^ isomer rather than a Ceo dimer. The 
collision energy range for which these fusion products are detected is also 
consistent with the gas phase measurements. The velocity distributions of the 
product ions in these experiments also support this interpretation. The 
projectile ion beam consisted of mass-selected Ceo^ so that the detected Cvo^ 
ions must have been sputtered from the substrate (the fullerite film contained 
10% C70) [24]. There is a low velocity component of the Ceo^ ions that 
matches the € 70 ^ distribution and that is also attributed to sputtering of the 
substrate material. The fusion products are significantly faster than the 
sputtered molecules and the velocity distribution is consistent with the 
formation of large fullerenes via fusion of the projectile ion with a 
surface fullerene in a direct reaction mechanism of the Eley-Rideal type 
[25]. 

Another interesting observation in fullerene-fullerite experiments is the 
occurrence of rainbow scattering. The rainbow effect, which can be observed 
in many kinds of particle scattering studies, causes enhanced scattering 
intensities at certain scattering angles, the so-called rainbow angles. It is 
named after the atmospheric phenomenon of rainbow arcs in the sky. In 
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general, all rainbows are the result of two counter-acting effects. For 
atmospheric rainbows, it is the interaction of the reflection and refraction of 




Figure 10.13. Time-of-flight mass spectra of the products of collisions between and (a) a 
fullerene film (b) graphite. Collision energy 275 eV and angle of incidence of 25”. Adapted 

from [24]. 

sunlight in water droplets that causes a minimum in the deflection angle. In 
atom-atom scattering it is the combination of the attractive and repulsive 
parts of the interaction potential. The angle of deflection changes from 
negative to positive as the impact parameter decreases. The minimum value 
of the angle of deflection is called the rainbow angle. According to the 
expression for the classical differential cross section (eq. (3.1.28)), at the 
rainbow angle dxidb = 0 and the classical deflection function goes to 
infinity. This shows up as a clear maximum in the experimental differential 
cross section. In surface scattering, the periodic surface potential gives rise 
to two extrema in the deflection function. The simplest case is scattering 
from a "zig-zag" potential where the reflection angle depends on which side 
of the "zig-zag" the projectile hits, see Fig. 10.14(a). The situation for 
fullerenes is shown schematically in Fig. 10.14(b). The primary angle of 
incidence, a, can be varied in the experiments while the angle yis fixed at 
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140°. Therefore the detected reflection angle is determined by /i = 180° - 
(y+a). The scattering plane is orthogonal to the fullerite surface and 
intersects it in the <110> direction [24]. The nearest neighbour distance in 
the Ceo lattice is 10.02 A in this direction. If the surface Ceo molecules 

a) b) a = Tl0° 





Figure 10.14. a) Illustration of specular reflection from a "zig-zag" potential. The particles are 
scattered at two angles: (11 and (12. (b) Scattering geometry for interaction of with a 
<1 1 1> oriented C^o crystal. The dot-dashed line shows an effective interaction potential. 
Scattering from this point under the condition a + (1 = 40° results in enhanced scattering 

intensity for a = -10°. 

interact with the projectile fullerene ion as hard spheres then the interaction 
process can be viewed as reflection of a point-like particle from the potential 
barrier indicated by the dash-dot line in Fig. 10.14(b). The dependence of the 
scattering signal on the angle of incidence a for a fixed y= 140° is shown in 
Fig. 10.15 for different final velocities (the impact energy was 275 eV) [24]. 
Three partially overlapping maxima are found at approximately a= -10°, 20° 
and 50° which shift slighdy as the detected velocity increases. The effect is 
even more pronounced for lower impact energies [18]. The calculations 
shown on the same figure come from an extremely oversimplified classical 
model and consider elastic reflection of a point mass from a rigid wall 
consisting of overlapping spheres (Fig. 10.14(b)). The calculations using a 
cosine form of the potential do indeed reproduce the pronounced rainbow 
structure. The central maximum can be attributed to lattice dynamics effects, 
as indicated by classical MD simulations [18]. 



2. SPUTTERING 

In the previous section it was pointed out that there is a strong signal 
from sputtered ions when a fullerene projectile ion impinges on a fullerite 
target [24]. This has also been observed for very high energy collisions of 
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fullerene ions (20 MeV) on bio-molecular samples [26]. For the 
desorption/ionisation of trypsin molecules, the secondary ion emission yield 
was found to be about two orders of magnitude higher than with MeV 
atomic ions. This is believed to be a consequence of the very high energy 
density deposited via electronic stopping in the solid by the highly energetic 




Angle of Incidence / degrees 



Figure 10.15. Angular distributions for C^o* scattered on a crystalline Cgo films for different 
final velocities and an impact energy of 275 eV. The two lowest plots are the results of 
numerical 2D simulations with simplified potentials. From [24]. 

fullerene ions leading to the ejection of a large amount of matter after each 
ion impact. At lower impact energies fullerene ions are a very surface- 
sensitive probe and could potentially be used for SIMS (secondary ion mass 
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spectrometry) applications [27]. The yield of ions per unit mass from organic 
surfaces was compared for small gold cluster ions and Ceo^ and found to be 
4-5 times higher for the fullerenes than for gold projectiles of similar mass 
[28]. A TRIM calculations suggests that 20 keV deposits most of its 
energy within 30 A of the surface whereas gold clusters of similar mass 
deposit their energy very much deeper in the solid (> 100 A). A comparison 
of the positive secondary ion yields for Ga"^, Au"^, Aus"^ and Ceo^ at 10 keV 
impact energy on poly(ethylene-terephthalate) is shown in Fig. 10.16 [27]. 
The advantage of using a fullerene ion beam is apparent. The damage cross 
sections for the organic materials under fullerene ion bombardment depend 
on their chemical and physical properties. High molecular weight polymers 
and monolayer films of molecular solids are damaged at approximately the 
same rate as under Ga"^ irradiation (an ion commonly used for SIMS 
applications). However, thick films of molecular solids show extremely low 
levels of damage and it appears that the molecules are being peeled off the 
surface with very little sub-surface damage. Generally, less fragmentation is 
seen with the fullerene projectile ions and there is a very large increase in the 
molecular ion yield compared to conventional projectile beams [27]. 




Figure 10.16. Comparison of positive secondary ion yield during analysis of PET using Au"^, 
Au 3 '^ and at 10 keV impact energy. The intensities are normalised to the signal from Ga"^. 

Adapted from [27]. 
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